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Chapter 1 
Introduction 



1.1 General remarks 

Up to now, no convincing solution for the problem of describing baryons and 
mesons in terms of a few underlying fields has been found, yet it is widely be- 
lieved that quarks and gluons are the basic entities which build the hadrons. The 
associated theory is called Quantum Chromodynamics (QCD). The complexity 
of this theory in the low energy sector is prohibitive for a straightforward deduc- 
tion of hadronic properties from it. Therefore, models with a simplified dynamics 
have been employed over the last decades which aim at describing hadrons with 
degrees of freedom borrowed from QCD (like quarks) or from the observable par- 
ticle spectrum (like pions) that might be better suited in the treatment of hadrons. 

The aim of the present thesis is to study a relativistic framework for the de- 
scription of baryons employing quark and diquark degrees of freedom. Actual 
calculations are performed with rather simple assumptions about these, but it is 
hoped that further progress in knowledge about QCD's quark propagator and the 
structure of 2-quark correlations will clarify whether the reduction to quarks and 
diquarks in the treatment of baryons is compatible with QCD. 

We refrain from discussing QCD here. The relevance of QCD and its (approx- 
imate) symmetries for hadronic physics is discussed in many good textbooks on 
quantum field theory, see e.g. refs. [1, 2]. 

In the remaining part of this Introduction we will give a short description of 
the basic ideas underlying three classes of hadronic models which have become 
popular (and have remained so) over the last years. This brief survey is by no 
means meant to be exhausting, and serves only as a motivation to consider a fully 
covariant approach like the diquark-quark model. A short overview about the 
topics which are covered in the subsequent chapters is given in the last subsection 
of this chapter. 
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1.2 Modelling Baryons 
1.2.1 Constituent quark models 

In the early sixties, it has become recognized that the eight spin- 1/2 baryons today 
referred to as "octet baryons" belong to an eight-dimensional tensor representation 
of the unitary group SU{3). Obviously this unitary symmetry is broken, as can 
be seen from the mass differences within the baryon octet, nevertheless a couple 
of useful relations among e.g. the octet masses and magnetic moments could be 
deduced from unitary symmetry [3] that appeared to be in fair agreement with 
experiment. 

Quarks with their denominations up, down and strange as building blocks 
of hadrons have been introduced somewhat later as the fundamental representa- 
tion of the flavor group SU{3). According to this hypothesis, the baryon octet is 
an irreducible representation which arises from the decomposition of the tensor 
representation formed by three quarks. Another irreducible representation, the 
decuplet of spin-3/2 baryons, has been identified in experiments in subsequent 
years. 

To match the baryon quantum numbers, the quarks had to be assigned spin 1/2. 
Spin and flavor may be combined to the larger group SU (6) and then octet and 

decuplet baryons can be interpreted as a 56-dimensional representation of this 
group, with a symmetric spin-flavor wave function. However, the overall wave 
function needs to be antisymmetric, and if one is not willing to sacrifice the spatial 
symmetry of a baryon ground state wave function, one is led to the hypothesis 
that all quarks appear in three colors. Baryons are then fully antisymmetric color 
singlets. 

Still, these group-theoretical considerations do not encompass a dynamical 
model for baryons in terms of quarks. To explain the plethora of baryon reso- 
nances, one of the very first models [4] aims at describing baryons by solutions of 
a 3 -particle Schrodinger equation using the non-relativistic Hamiltonian 



Massive quarks with ~ 0.3 GeV move in a harmonic oscillator potential which 
(i) confines the quarks, (ii) allows for the separation of variables easily and (in) 
leads to a rich spectrum of excitations, even far more than observed. However, 
the splitting between octet and decuplet has to be provided by some sort of spin- 
dependent interaction, from the theory of atomic spectra usually known as hyper- 
fine interaction. As meanwhile QCD became favored by many theorists to be the 



H — i?kin + -f^conf 



(1.1) 




(1.2) 
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theory accounting for the strong interactions, a tensorial and a spin-spin interac- 
tion term in the Hamiltonian could be motivated by a non-relativistic truncation 
of one-gluon exchange between quarks [5, 6], 



Here, ag is identified with the QCD fine structure constant and Si is the spin vector 
for quark i. Spin-orbit interactions which also result from the non-relativistic 
expansion of the one-gluon exchange are usually neglected using more or less 
convincing arguments as they lead to practical complications {e.g. three-body 
interactions) and spoil the calculated hadronic spectrum [6]. 

We have introduced the basic elements (i^kin, ^^conf and i^spin) of nearly all 
quark potential models that have been devised since the first models have been 
developed some thirty years ago. The conceptual problems associated with this 
approach are manifest: lack of covariance and not even an approximate chiral 
symmetry as the quarks appear with masses of around 0.3 GeV in the Hamiltonian. 
We shortly discuss approaches that have been taken to remedy this. 

In order to include relativistic effects, one certainly has to modify the kinetic 
energy operator. 



Fortunately, it is in this case still possible to separate overall bound state variables 
(total momentum and spin) from relative variables such that the following holds: 
the hereby constructed mass operator is only a function of the relative variables 
and the overall and relative variables (considered as operators) fulfill the commu- 
tation relations for the Poincare algebra separately [7]. However, Poincare invari- 
ance does not guarantee a covariant description within Hamiltonian dynamics. 
Rather, demanding covariance leads to fairly complicated spin-dependent con- 
straints on the Hamiltonian. They have been explicitly formulated for light-cone 
variables in ref. [8] but are also present for the above discussed conventional use 
of Galilean time as the Minkowski time variable in relativity. As these constraints 
are usually overridden in phenomenological applications, covariance is definitely 
lost in semi-relativistic potential models. 

Chiral symmetry and its dynamical breaking lead to massless Goldstone bo- 
sons. They are identified with the pseudoscalar mesons which do have mass. 





3 




(1.4) 



i=l 
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but this being comparatively small. Thus, chiral symmetry appears to be ex- 
plicitly broken only by small parameters which in QCD are the current quark 
masses. Within the last years, it has become popular to replace potentials mo- 
tivated through one-gluon exchange by non-relativistic approximations to one- 
meson exchange between quarks [9], to pay credit to the importance of pions. 
This does not render the Hamiltonian chirally symmetric, though, but solves some 
problems of the level ordering in the baryon spectrum. 

To conclude this brief summary of quark potential models, we remark that 
the practical impossibility of a covariant description limits these models to an 
effective parametrization of the spectrum and of some static observables. In a non- 
covariant description, the calculation of dynamical observables, with the simplest 
being form factors, introduce too much of a frame dependence into the results as 
that they could be trusted as model predictions. 

1.2.2 Soliton models 

We now turn to a short discussion about another kind of baryon model which 
in its original form also predates QCD and even the first quark models. In our 
exposition, we follow closely ref. [10]. The central idea is to regard baryons as 
solitons of an effective, non-linear meson theory. To illustrate this, let us start with 
a chirally symmetric theory describing massless pions and its chiral partner, the a 
meson. The so-called a model is defined by the Lagrangian density 

£ = i {d^TT-d^TZ + d^ad^'a) -t- V"(7r-7r + a^) . (1.5) 

The a field may be regarded as the fourth component of the pion field, and we 
see that this Lagrangian density is invariant under rotations in the (tt, a) space 
leaving tt^ -|- ct^ constant. Thus, the Lagrangian is invariant under transformations 
of the group 0(4) which in turn shares the same Lie algebra with the chiral group 
SU{2) X SU{2). Depending on the functional form for V, the model may exhibit 
also the feature of spontaneous symmetry breaking. The non-linear version of this 
model is defined by demanding 

C7^ TT^ = = const . (1.6) 

Introducing U — exp(iT-0), the Lagrangian simplifies to 

£(2) ^ ^rj^ df.Ud'^U^ , (1.7) 

since the potential V reduces to a constant and may be omitted. Now we look for 
static, finite energy configurations with a finite extension in space which minimize 
the energy - solitons. Due to the constraint in eq. ( 1 .6), tt is a variable which labels 
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the points of the three-dimensional unit sphere. A finite energy configuration only 

exists if U approaches a constant value for |a;| — > oo, with x being the position 
variable. Evidently, all points x on the sphere with radius r may be identified in 
the limit r — > cxo, and a finite energy configuration is determined by the map 

x^-nix) : S^^S\ (1.8) 

which is labelled by its topological invariants, the winding numbers. These con- 
siderations can be easily extended to time-dependent fields 7v{x,t) since the same 
boundary conditions as for the static configurations must hold at infinity. Thus a 
topological current may be defined, 

S'* = ^i^e'^^^^Tr [{U^dM){U^dpU){U^dxU)] , (1.9) 

which conserves the winding number as the corresponding charge. It has been 
the conjecture of Skyrme [11] to identify the topological current with the 
baryon current and the winding number with the baryon number. Unfortunately, 
one finds by scaling U (r, t) ^ U (Ar, t) that a configuration of minimal energy 
is only obtained for A — > oo, i.e. the soliton is not stable and collapses to size 
zero. This already indicates the way to resolve the problem, by adding a term to 
the Lagrangian (1.7) with more derivatives (but no more than two time-derivatives 
in view of later quantization) which stabilizes the soliton, 

C^'^ = y^Tr {[{U^d.U), {U^dM)] [{U^d^U), {U^d-'U)]) . (1.10) 

Then a static soliton configuration is given by the ingenious hedgehog ansatz, 

Uo^exp(iT-xF(r)) . (1.11) 

We see that the isovector pion field is simply proportional to the position vec- 
tor. The strength of the pion field is determined by F{r), the chiral angle. The 
boundary conditions for finite energy configuration may be translated into the re- 
quirements F{0) = rni and F{oo) = for a configuration with baryon number 
n. Thus we have found a classical description of baryons within the context of a 
simple non-linear meson theory. 

There arises the immediate question of how the baryons obtain their half- 
integer spin and isospin, since the fundamental pion field possesses spin zero and 
isospin one. To answer this question, the soliton needs to be quantized and to this 
end, time-dependent soliton solutions are needed. We do not want to go into detail 
here and just line out the basic ideas. Full time dependent solutions to the theory 
described by the Lagrangian C — C^'^^ + C^^^ are hard to come by. Nevertheless 
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one observes that global rotations in isospace, caused by a time-dependent matrix 



do not change the "potential" energy of the hedgehog, i.e. the part in the La- 
grangian depending only on spatial derivatives. Therefore they might be a suit- 
able ansatz for a time-dependent soliton solution.^ The time derivatives in the 
Lagrangian create terms which represent the rotational energy. Upon using the 
hedgehog properties and a suitable definition of angular velocities as canonical 
variables one finds [10]: 

• The absolute values of spin and isospin are equal, T^. 

• The Hamiltonian of the problem reduces to the one of the spherical top, 
H = ■^j'^ + M. Here, 6' is a moment of inertia and M is the static energy 
(the mass) of the soliton. 

Thus, upon quantizing the top, J and T as the quantum numbers for spin and 
isospin may be well assigned half-integer values, though there is no compelling 
reason to do this. 

Alas, one would like to extend the theory to three flavors and to chiral SU{3) x 
SU (3) symmetry in order to describe the baryon octet and decuplet. A detailed 
review on this subject is provided by ref. [12]. Let the Lagrangian C = C^'^^ +C^^^ 
remain unchanged but define the chiral field as an SU (3) field using the mesons 
0", a = 1, . . . , 8, built by pions, kaons and the octet component of the rj. 



One may wish to add some flavor symmetry breaking terms to the Lagrangian and 
can then proceed as described above, i.e. by introducing the collective rotations 

and quantizing the soliton. At some point one would like to relate the effective 
meson theory to QCD, and one could start by comparing symmetries. It turns out 
that this model possesses an extra discrete symmetry that is not shared by QCD. 
Under parity transformations P the pseudoscalar meson fields as described by 
QCD should obey 



A{t) e SU{2), 



U{x,t) = A{t)Uo{x) A^{t) , 



(1.12) 




(1.13) 



P: 4>%x,t) 



(1.14) 



The static soliton violates isospin or, equivalently, rotational symmetry. Therefore it needs to 
be projected onto good spin/isospin states. The introduction of the collective rotation in isospace 
exactly does this job. 
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In the meson theory, one therefore defines P as 



P : X 



(1.15) 



The Lagrangian is however invariant under cc — > — cc and U ^ separately. To 
break this unwanted symmetry, one needs to add some extra term to the meson 
action, the famous Wess-Zumino-Witten (WZW) term [13]. We state the conse- 
quences which come about by adding this term: 

• If photons are added in a gauge invariant manner to the meson theory, the 
WZW term generates a vertex for the anomalous decay tt*' — 77. This 
decay is determined in QCD by the triangle graph which is proportional 
to the number of colors, Nc- Thus, the strength of the WZW term can be 
related to Nc. 

• If one considers an adiabatic 2n rotation of the soliton during an infinite 
time interval, the WZW term gives a contribution Nqt^ to the action while 
all other contributions vanish. Thus the soliton acquires a phase (—1)^'^ for 
such a rotation. For Nc odd, the soliton must therefore be quantized as a 
fermion, in accordance with Nc = 3 in QCD. 

We find that incorporating the QCD symmetries (for three flavors) already gener- 
ates a good argument to identify baryons with solitons. 

There is one more connection of the soliton picture to QCD. Witten has shown 
in ref. [14] that in the limit Nc — ^ 00 QCD transforms into a theory of weakly 
interacting mesons and baryons emerge as solitons in this theory. However, an 
exact mapping has not been found in this limit. 

Instead, one can consider a Nambu-Jona-Lasinio (NIL) model involving ex- 
plicit quark degrees of freedom, defined by the Lagrangian [15] 



Here, quarks (g) are assumed to come in 3 flavors and Nc colors, rrr denotes a 
current quark mass matrix which is assumed to contain small diagonal terms and 
are the generators of the color group SU (N). The model is an approximation 
to QCD if one assumes the gluon propagator to be diagonal in color space and 
~ 5*{x). Additionally, in the limit mP = it conserves the chiral SU{3) x 
SU (3) symmetry of QCD (with massless current quarks) which is of importance 
for meson and baryon properties. Upon Fierz transforming the NIL Lagrangian 
into attractive channels, one can show that in the limit Nc 00 only the meson 
channels survive indeed. 




(1.16) 
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Furthermore, this specific model and other quark models of the NJL type ad- 
mit solitonic solutions which are reviewed in refs. [16, 17]. In contrast to the 
non-linear a models discussed before, these solitons are non-topological. In mini- 
mizing the energy of a field configuration, one determines the quark energy levels 
in a background meson field of the hedgehog type, and a unit baryon number 
is still provided by a configuration of Nc quarks. It is an advantageous feature 
to retain quark degrees of freedom in the solitonic picture, if one considers phe- 
nomenological applications such as calculations of structure functions and parton 
distributions which have been performed in the last years. 

We will stop at this point and summarize some important features of soliton 
models. They possess a deep conceptual appeal as with the requirement of chiral 
symmetry imposed onto some effective meson theory the baryons arise as lumps 
of these without further ado. The possible transformation of QCD in the limit 
Nc oo into such a nonlinear meson theory creates the link to the fundamental 
quarks. However, in practical calculations one has to abstract from quarks entirely 
or must resort to NJL models as crude approximations to QCD. For both strategies 
it turns out that in leading order of Nc (quasiclassical soliton configuration with 
quantized collective variables) static baryon observables are only accurate on the 
20 . . . 30 % level, as known from nearly all hadronic models. Corrections due 
to the finite number of colors (or, phrased otherwise, due to mesonic fluctuations) 
yield substantial contributions, and their technically involved calculation tends 
to hide the original conceptual beauty. Another point of weakness is the non- 
covariant formulation which will hamper the possible application of this model to 
hadronic processes in the intermediate energy regime. 

1.2.3 Bag models 

Soon after QCD surfaced in theory circles, it became popular to associate hadrons 
with a different phase of the theory as compared to the vacuum phase. This picture 
motivates the original M.I.T. bag model, where the vacuum phase is assumed to 
prohibit the propagation of quarks and gluons but creates bubbles of hadronic size 
in which quarks and gluons may propagate ordinarily. The model is reviewed in 
ref. [18]. It is formalized by the Lagrangian 

>Cbag = (>Cqcd -B)ev- ^qqAs . (1.17) 

Inside the bag, to be specified, Oy = 1, and outside, 9v = 0. The surface function 
As is defined by |^ = n^As and n'^ — (nP, n) is a normalized unit vector 
perpendicular to the surface. According to this setup, no quark or gluon current 
can escape the bag, due to the boundary conditions (for a static bag) 

nxBc^O, n-Ec^O, qi^q^O, (1.18) 
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where Be and Ec denote the color magnetic and electric field, respectively. The 
bag constant B has been introduced to balance the surface pressure of the confined 
color fields. 

The usual procedure in finding the ground state is to assume a spherical, static 
bag and evaluate the quark wave functions and energy levels perturbatively. To 
this end, one solves the free Dirac equation with the spherical boundary condi- 
tions, and the lowest-lying energy level for a massless quark is approximately 400 
MeV for a bag radius of 1 fm, slightly above the constituent quark masses used in 
the potential models. The splitting between nucleon (spin 1/2) and A (spin 3/2) is 
evaluated by a one-loop gluon correction to the "bare" quark levels. In the end one 
finds that the bag radius usually comes out too large, above 1 fm, and for the QCD 
fine structure constant > 2 to account for the — A mass difference. Hence 
the validity of a perturbative treatment is highly questionable. Furthermore, pionic 
bags give just the wrong pion mass. 

This led to the idea to couple a mesons and pions to the bag in order to regain 
chiral invariance. These models come under different guises as 

• a "hybrid" bag where quarks and a — n couple on the bag surface and the 
latter are also allowed to propagate inside the bag, 

• a "little" bag where in contrast to the hybrid bag pions are excluded from 
the interior of the bag, 

• a "cloudy" bag where the mesons are constrained to the chiral circle and are 
allowed inside the bag, and 

• a "chiral" bag where the constrained mesons are described outside the bag 
by the Skyrme Lagrangian C = C^"^^ + C^'^'' which was discussed in the 
previous subsection, cf. eqs. (1.7,1.10). 

All these facets of a softened bag are described on an introductory level in ref. [19]. 
Usually all these models improve on the shortcomings mentioned above. The bag 
radius can be chosen considerably smaller and as needs not be as large since the 
— A splitting arises partially from the pions. Additionally, they support the 
popular picture of a nucleon as consisting of a valence quark core surrounded by 
a pionic cloud. As it will be of interest later on, the non-zero electric form fac- 
tor of the neutron is frequently attributed to "the" pion cloud around the neutron 
but such a statement has to be interpreted with care as it is only valid in a model 
context such as a bag model, where valence quarks alone yield zero for the form 
factor indeed. 

Especially the chiral bag mentioned last in the above list aims at combining 
the soliton picture of baryons with the valence quark picture. At first sight, it is ap- 
pealing to think of the bag surface as a borderline which indicates where to switch 
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from a description with quark degrees of freedom to a description with mesons. 
Then observables should depend very little on the size of the bag ( Cheshire Cat 
Principle) and this has been found indeed. Unfortunately the static properties 
of baryons are not improved as compared to the pure soliton, and come out even 
worse than for solitons which are formulated with pseudoscalar and vector mesons 
[12, 20]. Non-covariance of the formulation and possible large modifications of 
observables by quantum fluctuations also remain present as problems to be solved. 

1.2.4 The diquark-quark picture: baryons as relativistic bound 
states 

Hopefully, the short review of the basic principles of three established classes 
of baryon models has illustrated the conceptual variance that exists among the 
models. Yet, all of them describe static and a few dynamic observables fairly 
well on the proverbial twenty-percent- level of accuracy, and their degrees of free- 
dom may be motivated from underlying QCD and its symmetries. Nevertheless 
we find that none of these is formulated in a manifestly covariant manner, and a 
model which deliberately cares about covariance might bridge the gap between 
soft and high-energy physics. Experimental progress in this intermediate energy 
regime has been considerable within the last years, and therefore understanding 
the corresponding results within a hadronic model with substructure constitutes a 
good motivation. 

In the next chapter, the framework of a diquark-quark model will be described 
in detail. In this model, we identify baryons with poles in the 3-quark correlation 
function. These poles may arise by summing over infinitely many interaction 
graphs, and it is well known that this summation leads to bound state equations 
such as the Bethe-Salpeter equation which is usually formulated for two particles 
but can be extended to n particles as well. We will reduce the 3-quark problem 
to an effective two-particle problem by introducing quasi-particles (diquarks) as 
separable 2-quark correlations. The corresponding Bethe-Salpeter equation can be 
solved without any non- or semi-relativistic reduction. For the model calculations, 
we employ in a first attempt massive constituent quarks with the corresponding 
free-fermion propagator. It is certainly reasonable to assume that a constituent 
quark mass is generated in QCD, and this mass generation is reflected in the quark 
propagator [21], but we acknowledge that the full structure of the quark propagator 
in QCD is certainly described inadequately by just a free, massive propagator. 

In chapter 3, we calculate the nucleon electromagnetic form factors maintain- 
ing covariance and gauge invariance which demonstrates the feasibility of a fully 
relativistic treatment also for dynamic observables. The electric nucleon proper- 
ties can be described very well, we emphasize here in anticipation of the results 
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that the behavior of the neutron electric form factor can be understood without any 
explicit reference to pions. The magnetic moments turn out to be underestimated, 

at this stage. 

The results for the pion-nucleon form factor and the isovector axial form fac- 
tor, presented in chapter 4, are, for the latter and away from the soft point, in 
accordance with the available data. For the pion-nucleon form factor, our predic- 
tions converge with the predictions of other baryon models. However, the results 
reveal an overestimate for the pion-nucleon and the weak coupling constant which 
are obtained in the soft limit. Chiral symmetry demands here a further extension 
of the model to include another diquark correlation, but this proved to be beyond 
the scope of this work. 

In contrast to the baryon models described above, confinement is not included 
by construction in the diquark-quark model. In chapter 5 we investigate the pos- 
sibility of incorporating confinement by a suitable modification of quark and di- 
quark propagators. Thereby we are in the position to calculate the mass spectrum 
of octet and decuplet baryons, and also the nucleon magnetic moments receive 
some improvement. The justification of these modifications is discussed thor- 
oughly. 

Chapter 6 is concerned with a critical comparison to a semi-relativistic solu- 
tion of the model which is shown to be inadequate. Finally, we summarize the 
obtained results in chapter 7 and comment upon further perspectives in the de- 
scription of baryons within the context of the diquark-quark model. 

Most of the results presented here have been published previously, cf. refs. [22, 
23, 24, 25, 26]. Also, the author of this thesis has benefited a lot from the collab- 
oration with the authors of ref. [27] and ref. [28]. 

Throughout the text we work in Euclidean space, if not indicated otherwise. 
The conventions relevant in this connection are summarized in appendix A.l. 
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Chapter 2 



The Covariant Diquark-Quark 
Model 



2.1 Reduction of the relativistic three-quark prob- 
lem 

The three-body problem has attracted a great deal of attention over the last decades 
and the formal framework presents itself to us as being quite established, for an 
introduction see refs. [29, 30]. Most work in the field has been restricted to the 
treatment of quantum mechanical systems, though, where certain theorems about 
existence and uniqueness of solutions could be proved, see e.g. ref. [31] for the 
case of three-particle scattering with just two-body potentials. For relativistic sys- 
tems the few-body problem seems to be somehow ill-posed as formal require- 
ments like covariance and locality have led to the formulation of quantum field 
theories which employ infinitely many degrees of freedom (that come about by 
the possible creation and annihilation of particles). Nevertheless, the spectrum 
of mesons and baryons, their static properties and transitions indicate that only 
the degrees of freedom of a fixed number of (constituent) quarks are relevant. 
Therefore let us introduce the notion of baryon wave functions as matrix elements 
of three quark operators that interpolate between the vacuum and a bound state 
with observable quantum numbers (see below eq. (2.7)). Within QCD, the non- 
perturbative vacuum is a non-trivial condensate and therefore these wave functions 
would contain per definitionem sea quark and gluonic parts. However, it seems 
to be justified phenomenologically to regard wave functions which are matrix ele- 
ments of three effective quark operators between the perturbative vacuum and a 
bound state as the dominant ones compared to other operator matrix elements that 
involve an arbitrary number of particle creation and annihilation operators. Then 
the quantum mechanical formulation of the three-body problem can be taken over 
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to field theory in the formulation of Green's functions. 

2.1.1 Distinguishable quarks 

First, let us consider three distinguishable quarks. Later we will generalize the 
obtained results to indistinguishable particles. 

To save on notation, we will drop Dirac, flavor and color indices. They are 
assumed to be contained in the single particle labels. In our definitions of Green's 
functions or bound state matrix elements we always take out one 5-function rep- 
resenting conservation of the total momentum, e.g. we will denote a (dressed) 
single quark propagator by Si with 

x{0\Tqi{xkMypM (2.1) 
=: F.T. {0\Tqi{xkMypM ■ (2-2) 

We will use the abbreviation 'F.T. (x)' for the Fourier transform of the subsequent 
expression (x) in the following. The annihilation operator for a quark of kind i is 
denoted by qi. In this way, the disconnected (free) three-quark propagator is found 
to be 

GoikuPi) = S\h-p,)S\h-p2)l[S,{kf,Pi) (2.3) 

Sum and product run over i = 1, 2, 3. The full quark 6-point function (or the 
3-quark correlation function) G is accordingly 

{2n)^d'(j2ih-Pi)) G{h;pi) = F.T (0|r J]?,(xfeJg,(y^J|0) . (2.4) 

It obeys Dyson's equation, symbolically written as 

G^Go + GqKG . (2.5) 

K is the three-quark scattering kernel that contains all two- and three-particle 
irreducible graphs and the whole equation involving the symbolic products reads 
explicitly 

G{h;pi) = Go{h;pi)+ (2.6) 

/d^l f dH f d^m f d^m 

Besides the four-dimensional integrations, sums over Dirac indices, color and fla- 
vor labels are implicitly understood in the symbolic multiplication. 
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We define a three-particle wave function ijj to be the transition matrix element 
between the vacuum and a state with an on-shell four-momentum P and appro- 
priate discrete quantum numbers (spin, isospin, . . . ), 

(27r)V (J2p^ - P) HpuP2,P3) = F.T. {0\qi{x,)q2{x2)qs{xs)\P) . (2.7) 

A bound state of mass M with wave function ip will show up as a pole in the 
6-point function, i.e. in the vicinity of the pole G becomes 

^ ^ p^TJp ' ^^-^^ 

where P — pi +P2 +P3- Inserting this into Dyson's equation (2.5) and comparing 
residues, one finds the homogeneous bound state equation 

il) = GoK'4) ^ G-^'iIj = 0. (2.9) 

This is our starting point. As it stands, the equation is far too complicated as nei- 
ther the fully dressed quark propagator nor an expression for all two- and three- 
particle irreducible graphs are known. One has to resort to some sort of approx- 
imation in which the problem is tractable and judge the quality of the chosen 
approximation by the results for the bound state properties. 

The problem is greatly simplified when one discards three-particle irreducible 
graphs from the interaction kernel K. We will call this approximation the Faddeev 
(bound state) problem in the following. The kernel then consists of three terms, 

K^Ki + K2 + Ks. (2.10) 

The Ki, i — 1,2, 3, describe the interactions of quark pairs (jk), i.e. with quark (i) 
as a spectator, (ijk) is here a cyclic permutation of (123). In figure 2.1, examples 

for admitted and excluded graphs in the kernel K for the Faddeev problem are 
shown. The two-quark propagators gi fulfill their own Dyson's equation with 
respect to the kernel K^, 

g.^Go + GoK^g,. (2.11) 

The objects g^ and Ki are defined in three-quark space. The former contain a 
factor Si, the propagator of the spectator quark, and the latter contain a factor of 
(although the spectator quark is not involved in the interactions described by 
Ki). The matrix ii is defined by amputating all incoming and outgoing quark legs 
from the connected part of gi, 

gi^Go + GoiiGo . (2.12) 
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Figure 2.1: Examples for admitted and excluded graphs in the 3-quark interaction kernel 
K for the Faddeev problem. 

Combining the two previous equations yields an integral equation for ii, 

U = Ki + K,GoU . (2.13) 

Still this matrix contains remnants of the spectator quark, this time S~^. Therefore 
the two-quark t matrix of the interacting quarks (jk) is obtained by 

ti = kS,, (2.14) 

and contains just a trivial unit matrix in the subspace of the spectator quark. 

We need one more definition before we can tackle the bound state problem. 
The so-called Faddeev components tpi are introduced by 

^/j, = GoK,ij, (2.15) 

and by virtue of eqs. (2.9,2.10) = ^ipi. Solving eq. (2.11) for GoKi and 
inserting the result into eq. (2.9) yields 

i; = g, {K, + Kk)^. (2.16) 

Using the definition of ii and U in eqs. (2.12,2.14) the equation can be rewritten 
as 

ipi = Go ii {ipj + ipk) = (SjSk) U {ipj + ipk) ■ (2.17) 

These are the famous Faddeev bound state equations relating the Faddeev com- 
ponent ipi to i^jii^k using the full two-quark correlation function ti (instead of 
the kernel K). In turn the tj would have to be determined in a full solution to 
the Faddeev problem by solving eq. (2.13) which involves the kernel components 
Ki. The relativistic Faddeev equations, depicted in figure 2.2, are a set of cou- 
pled 4-dimensional integral equations and are thus a considerable simplification 
of the original 8-dimensional integral equation problem defined in eq. (2.9). Still, 
the wave function components ifji depend on the two relative momenta between 
the three quarks and expanding these components in Dirac space [32] leads to an 
enormous number of coupled integral equations. 
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Figure 2.2: The Faddeev bound state equation for the component ^i. The equations for 
■02 and ^3 follow by cycUc permutation of the particle indices. 

Therefore, we aim at further simplification of the bound state problem. This 
will be achieved by approximating the two-quark correlation function ti in terms 
of a sum over separable correlations, 

^^(^l,^2;Pl,P2) =^Xn^l,^2) Dl^\h + k2) mPl,P2) , (2.18) 

a 

pictorially shown in figure 2.3. Separability refers to the assumed property of 
ti that it does not depend on any of the scalar products ki -pj = 1, 2). We 
call these separable correlations "diquarks". The function x" is the vertex func- 
tion of two quarks with a diquark (labelled with a), x"' correspondingly denotes 
the conjugate vertex function.^ Different a's refer to different representations in 
flavor and Dirac space, and we name D*^"-' the "diquark propagator" for the re- 
spective diquark. The assumption of separability in the two-quark correlations is 
the starting point in the description of baryons within the diquark-quark model. 
In the following section we will introduce scalar and axialvector diquarks as the 
supposedly most important correlations and explore the consequences of quark 
antisymmetrization for their parametrization. 

A suitable ansatz for the Faddeev components ipi, see figure 2.4, takes the 
form 

MPi^Pj^Pk) = ^Go XiiPj^Pk) Df\pj + pk) (j)i{pi,Pj + Pk) , (2.19) 

a 

where we introduced an effective vertex function of the baryon with quark and 
diquark, 0". It only depends on the relative momentum between the momentum 

^We refer here to the conjugate of a Bethe-Salpeter vertex function. For the scalar and axi- 
alvector diquarks which will be of interest, it is introduced in eqs. (2.48,2.49). 




Figure 2.3: The separable matrix ti. 
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Figure 2.4: The ansatz for ij^i using effective baryon-quark-diquark vertex functions 

of the spectator quark, pi, and the momentum of the diquark quasiparticle, pj +Pk- 
A detailed derivation of this property of 4>f can be found in refs. [29, 33]. 

Inserting this ansatz into the Faddeev equations (2.17) yields coupled integral 
equations for the effective vertex functions, 

C = X- {S,S,) Y.^x]Dfct>] + xlD^Ul) . (2.20) 
h 

We will rewrite these equations into a form that will be familiar from the theory 
of Bethe-Salpeter equations [34, 35]. To this end we define a relative momentum 
between quark i and the diquark made of the quarks with label {jk), 

pf^{l-n)p,-rj{pj+pk), (2.21) 

and the total momentum 

P = Pi+P2+P3. (2.22) 

In contrast to the separation of variables in a non-relativistic two-body problem, 
we have no unique definition of the relative momentum between quark and di- 
quark. Therefore we have introduced the parameter rj which distributes relative 
and total momentum between quark and diquark.^ Physical observables, like the 
mass of the bound states, form factors etc. should of course not depend on this 
parameter. 

Using the last two definitions for total and relative momentum and taking into 
account all momentum conservation conditions in eq. (2.20), the bound state equa- 

^Let Xg and Xd be the coordinates of quark and diquark in configuration space. Translational 
invariance tells us that 4>f{xq + h,Xd + h) = exp(— iP-/i)0"(a;g, Xd). Therefore one introduces 
a relative coordinate x'^'^^ = Xq — Xd and an overall configuration variable X = r]Xq + [1 — 'q)xd. 
The momenta p™' and P are the conjugate variables to .t''''' and X. Then (f>°;{xq + h, Xd + h) = 
(j)'l{x^^^, X + a) and, after a Fourier transformation, translational invariance is seen to be fulfilled 
with this choice of coordinates, leaving the freedom to choose r] arbitrarily. 
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tions read 

(2.23) 

Gq^jI^'^ describes the disconnected ("free") quark-diquark propagator in the di- 
quark channel b, 

^0 ' P) = SjivP + ) - V)P - kf) , (2.24) 

and K^^^jj is the quark-diquark interaction kernel given by 

3 

KS%iPT\kf,P) = J2~^iJkmkJ,q)Sk(q)x%q,P^), (2.25) 

k=l 

kj = vP + kf, (2.26) 

Pi = vP + pf, (2.27) 

q = -pY' - kf + {1 - 2r])P . (2.28) 

The symbol Sijk is equal to 1 if i, j, A; are pairwise distinct and otherwise. 

Clearly this Bethe-Salpeter interaction kernel A'^^t describes the quark ex- 
change between the two-particle configurations of quark i/diquark a and quark j/ 
diquark b. To arrive at an equation for physical baryons, we have to project these 
configurations onto the respective baryon quantum numbers. This will be done 
for the nucleons and the A resonances in section 2.3. As it will turn out, after 
the projection of the vertex functions the quark exchange generates the attractive 
interaction that binds quarks and diquarks to baryons. 



2.1.2 Identical quarks 

For identical quarks antisymmetrization is required. Nevertheless, the Bethe- 
Salpeter equation (2.23) will take the same algebraic form regardless of the parti- 
cle index i. Single particle indices i on the quark propagators Si and the diquark 
propagators dI"'^ can be omitted now as well as on the vertices xt their func- 
tional form does not depend on them. However, some care is needed since the par- 
ticle indices i,j, k in eq. (2.23) specified the summation order over color, flavor 
and Dirac indices. Therefore we will introduce for each term in the Bethe-Salpeter 
equation for identical quarks small Greek multi-indices a,i3,... for quarks that con- 
tain all of the indices mentioned above. Consequently, antisymmetry demands 
that the diquark-quark vertices x%-y be antisymmetric in their indices p^^. 
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Figure 2.5: The coupled set of Bethe-Salpeter equations for the effective vertex func- 
tions 

The Bethe-Salpeter equation for 0" = 0" in the case of identical quarks takes 
the form 

/jA / rcl 
j^,Kl^'{k^^\p-\P)x (2.29) 

GlJ^{k'''\P) = Sy.ik,) 2D^'\kd) , (2.30) 
K^^,{k'^\p'''\P) = xW{k„q)S^,^{q)xU<l^P,). (2.31) 

Repeated indices are summed over. Please note that, to save on notation, we 
neglected the indices for the involved diquark channels in the symbols Gg^'^ and 
K^^. Due to the summation over 6ijk an overall factor of two emerged that we 
have to absorb in the definition of the diquark propagator in the case of identical 
quarks [33]. Spectator quark momenta kq\pq\, diquark momenta kd\pa] and the 
exchange quark momentum q are given by 

kM = vP + k'V], (2.32) 
kM = {1 - rj)P - k'^'lp'^'] , (2.33) 
q = -k'^^ - p^^^ + (^i - 2rj)P . (2.34) 

The set of Bethe-Salpeter equations is pictorially shown in figure 2.5. 

To conclude this section, we will write down the completely antisymmetric 
Faddeev amplitude F, obtained by cutting off the three external quark legs from 
the wave function ^ (cf. eq. (2.19)), 

even 
perms 

a {ijk) 

(2.35) 
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The relative momentum pp' is defined as in eq. (2.21). Due to the antisymmetry 
of the diquark-quark vertices the summation over even permutations of (ijk) 
already generates the fully antisymmetric Faddeev amplitude. 

Summarizing the results of this section, we have shown that the determination 
of a 3-quark wave function defined by eq. (2.7) can be simplified enormously by 
the following approximations: 

1. Neglect all three-particle irreducible contributions to the full three-quark 
correlation function G. 

2. Assume that the truncated, connected two-quark correlations (the t matri- 
ces) can be represented as a sum over separable terms that are identified 
with diquark-quark vertex functions, see eq. (2.18). 

The bound state equation (2.9) reduces to a coupled set of Bethe-Salpeter equa- 
tions for effective baryon-quark-diquark vertex functions 0", where a labels the 
diquark types that are summed over in the expression for the separable t matrix, 
eq. (2.18). These Bethe-Salpeter equations (2.29) are the main result of this sec- 
tion and its only ingredients are 

• the quark propagator S, 

• the diquark-quark vertices x" and 

• the diquark propagators D^"-\ 

The aim of this work is to explore the consequences of choosing rather simple 
forms for S, x" and D'^"-^ that in a more elaborate approach should be determined 
by a suitable effective theory of QCD. In the diquark sectors of Nambu-Jona- 
Lasinio (NJL) type of models [36], the corresponding bound state equations have 
been formulated first in ref. [37] and solutions of the full Bethe-Salpeter equation 
have been obtained in ref. [33]. However, the framework of an NJL model is too 
tight: quark and diquark are not confined and additionally the diquark is pointlike 
(x" X°'{Pi:P2) where pi[2] are the quark momenta). Furthermore, the framework 
presented here is well-suited for the calculation of dynamical observables. In 
the case of electromagnetic observables, gauge invariance is also maintained in 
numerical solutions to a high degree of precision. This is somewhat harder to 
achieve in NJL model calculations employing sharp Euclidean cut-offs. Partly due 
to this limitation, no dynamical observables of baryons have been calculated in the 
diquark sector of the NJL model so far, using fully four-dimensional solutions of 
the Bethe-Salpeter equation. 
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2.2 Diquark correlations 

As explained in the previous section, "diquarks" or "diquark correlations" refer 
to the use of a separable two-quark correlation function, the t matrix. As the 
diquark-quark model presented here is inspired by the NJL model, we will de- 
scribe shortly how separable t matrices arise in it and discuss previous approaches 
to the diquark-quark picture of baryons. The following main part of this section 
introduces our parametrization of the t matrix including the explicit form for the 
diquark-quark vertices and the diquark propagators. 

Before we embark upon these points, let us discuss heuristically some prop- 
erties of diquarks. Any two quarks within a baryon have to be in a color anti- 
symmetric state as required by the Pauli principle and the assumed color singlet 
nature of baryons. If we consider for the moment diquarks as asymptotic states (in 
an unphysical subspace of the whole Hilbert space), their spatial wave functions 
(assumed to be separated from their Dirac structure) will be dominated by a sym- 
metric ground state. This requires the combination of flavor and Dirac structure 
to be symmetric. Their mass spectrum and therefore their importance as virtual 
states in baryons can be guessed by looking at the meson spectrum: pseudoscalar 
mesons correspond to scalar diquarks, vector mesons to axialvector diquarks since 
the intrinsic parity of (qq) is opposite to the intrinsic parity of (qq). Although the 
pseudoscalar mesons are especially light due to their suspected nature as Gold- 
stone bosons, the scalar diquarks corresponding to them should be still lighter 
than axialvector diquarks. This point gets support from lattice calculations [38]. 
In quenched QCD employing Landau gauge the scalar diquark correlation func- 
tion corresponds to a massive state of around 700 MeV whereas the axialvectors 
are heavier by around ICQ MeV. Vector and pseudoscalar diquark masses are then 
expected to lie around 1 GeV and therefore assumed to be less important for nu- 
cleons. This is indeed the case in a calculation within the Global Color Model 
[39]. 

Having motivated the importance of scalar and axialvector diquarks, let us 
return to their description within the 2-quark correlations. As is well-known, sep- 
arable t matrices arise from separable 4-quark scattering kernels K = KiSi in the 
following way^. Let 

h;puP2) = X^ih) ^^^\h + k2) riPi) , (2.36) 

be the scattering kernel in the particular channel a. Then the corresponding t 
matrix in this channel can be inferred from eq. (2.13) which for identical particles 

^The kernels Ki introduced in eq. (2.10) contain a factor S~^, the inverse propagator of the 
spectator quark. 
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is modified to 

t = K+^K G^o^t . (2.37) 

(2) 

We denote as the free two-quark propagator, as for distinguishable particles it 
corresponds to just the product of the single propagators. Due to the antisymmetry 
of K in the labels (jA), t is also antisymmetric in these labels. Its solution for the 
separable kernel i^^'^P " reads 

t''{h\v^) = x^{K) D^'^^h + h) r{Pi) ■ (2.38) 
The inverse of the propagator D'^"'^ is determined by 

(^W)-i=(AW)-i_l^aG(V. (2.39) 

We use the symbolic product notation as in the previous chapter. The last term 
in eq. (2.39) describes a quark loop with insertions of the vertex functions x" and 

r- 

In QCD the simplest term in the perturbative kernel, the gluon exchange, has 
already a non-separable form. However, when approximating the gluon propaga- 
tor by a 5-function in configuration space, the resulting model is a Nambu-Jona- 
Lasinio (NIL) model with pointlike four-quark interactions that can be viewed as 
color octet-flavor singlet quark current interactions. This local current-current in- 
teraction can be rewritten into attractive interactions in color singlet meson and 
color triplet diquark channels [15]. Let us pick the quark-quark scattering kernel 
in the scalar channel s which we write in the form 

{K%p,,s = ^Gs{r)ap{x%5 (2.40) 
= AGs il'C \%p (C^7^ A^),^ , (2.41) 

where we consider just an isospin doublet of two quarks that are in a flavor an- 
tisymmetric (r^) and a color antitriplet state (A''', k = 2, 5, 7). The r* are Pauli 
matrices and the A'^ are Gell-Mann matrices [40]. The coupling constant Gg regu- 
lates the strength in this interaction channel and C denotes the charge conjugation 
matrix. The scalar diquark propagator is by virtue of eq. (2.39) 

{D^-r' ik') = 4^ - 2 / (04(^V) S{q + k/2) (7^C) S^{k/2 - q) . 

(2.42) 

After suitable regularization of the divergent integral it turns out that the inverse 
propagator has zeros for certain values of k'^, therefore bound scalar diquarks exist 
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in the NJL model and the propagator can be viewed as a scalar propagator which 
contains the effect of the dressing quark loop. 

In analogy to the meson spectrum where pseudoscalar mesons are the lightest 
ones, followed by vector mesons, axialvector diquarks are expected to be impor- 
tant in the diquark channel. In fact, the axialvector diquark propagator acquires 
poles [41], but their appearance depends on the ratio of the coupling constant in 
the diquark channel to the coupling constant in the pseudoscalar meson channel. 

Having obtained separable t matrices for the scalar and the axialvector cor- 
relations, the spectrum of spin- 1/2 and spin-3/2 baryons is in principle calculable 
from the Bethe-Salpeter equation (2.29). This was done in refs. [42, 43], however, 
for the quark exchange kernel a static approximation was used (which amounts to 
neglecting all momentum dependence in the propagator of the exchange quark). 
In doing this, covariance is lost, and the rich relativistic structure of the effec- 
tive nucleon-quark-diquark vertex functions, to be described in section 2.4, has 
been furthermore approximated by the dominating non-relativistic components, 
the s waves. The full solution of the Bethe-Salpeter equations, i.e. with scalar and 
axialvector diquarks, for nucleon and A was obtained in ref. [33]. Subsequent cal- 
culations of static nucleon observables employed unfortunately only the solutions 
in the scalar diquark sector [44, 45]. Especially the magnetic moments turned out 
to be deficient, thus calling for including the axialvector diquark channel. 

The approximation of the two-quark t matrix by a finite sum of isolated poles 
at timelike momenta k'^ seems to be supported by these calculations. Solutions 
for on-shell diquarks were also obtained in the ladder- approximated Global Color 
model [39]. The Global Color Model employs a non-local current-current interac- 
tion, mediated by an effective gluon propagator and therefore avoids the regular- 
ization problems present in the NJL model. Parametrizing the t matrix with only 
the hereby obtained scalar diquark pole the nucleon Bethe-Salpeter equation was 
solved in ref. [46] . 

However, when going beyond ladder approximation, i.e. including higher or- 
der perturbation graphs in the quark-quark scattering kernel, the diquark poles 
disappear in the NJL model [47] as well as in the Munczek-Nemirowsky model 
[48]. The latter employs for the gluon propagator a (5-function in momentum 
space and can be regarded as quite complementary to the NJL model, indicating 
that in higher-order calculations with a more realistic gluon propagator diquark 
poles may disappear from the quark-quark scattering amplitude too. 

Nevertheless, the phenomenological efficacy of using a diquark correlator to 
parametrize the 2-quark correlations does not rely on the existence of asymptotic 
diquark states. Rather, the diquark correlator may be devoid of singularities for 
timelike momenta, which may be interpreted as one possible realization of diquark 
confinement. In principle, one may appeal to models employing a general, separa- 
ble diquark correlator which need not have any simple analytic structure, in which 
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case no particle interpretation for the diquark would be possible. Indeed, the fact 
that the interaction mechanism in the quark-diquark picture (the quark exchange) 
is attractive, turns out to be independent of the details of the employed diquark 
correlations. It can be traced back to the antisymmetry of the diquark and baryon 
vertices and arises thus from the color and flavor factors in the Bethe-Salpeter 
equation, see the next section. Last but not least, the usefulness of this picture and 
of employing specific diquark correlators has to be judged by comparison to the 
experiment. Calculations of observables have to go beyond nucleon and A mass 
or the octet and decuplet spectrum as these are rather insensitive to the level of 
sophistication in the treatment. 

As the results for the t matrix have such a strong model dependence, we will 
not resort to a specific effective quark model for calculating it in the following. 
Rather, we will employ an admittedly simple ansatz for the t matrix, the pole 
approximation in the scalar and axial vector channel. We will assume an isospin- 
doublet of u and d quarks as we are mainly interested in nucleon properties. Later 
on, in chapter 5, we will generalize our treatment to the case of broken SU (3) fla- 
vor symmetry. We will use vertex functions with a finite extension in momentum 
space, first because if strong diquark correlations exist their size is most likely of 
the order of the proton radius, and secondly, to avoid artificial regularizations in 
solving the bound state equations. The ansatz for the 2-quark t matrix reads, 

t{ka, h;Pa:Pp) = t{k,p, P) = xlp{k. P) D{P) X^p, P) + 

x'apik,P)D>^^{P)r,5{P.P)- (2.43) 

The relative momenta are defined as 



k\p\ = a ka\pc\ - (1 - (t) kp\pi3i a e [0, 1] . (2.44) 

The diquark propagators used in the scalar and the axialvector channel are taken 
to be 

D{P) = -^^^C{P\m^msc), (2.45) 
D^^iP) = --^5^^ + (1 - C{P\ m = m„.) . (2.46) 

"T '"'ax \ '"'ax / 

The choice C{P'^, m) = 1 and ^ = corresponds to free propagators of a spin-0 
and a spin-1 particle, to be employed here. Most generally, the dressing function 
C is not necessarily the same for scalar and axialvector diquarks. However, when 
investigating nontrivial forms for C (which mimic confinement) in chapter 5 we 
will assume the same C in both channels for simplicity. 
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For the trivial form, the diquark-quark vertices x ^iid correspond on-shell 

(P^ = —rrvlcyax]} diquark Bethe-Salpeter vertex functions. With quark legs at- 
tached, they are Fourier transforms of the transition matrix element 

S{aP + p)xZ\p^P)S^{{l-<j)P-p) = (2.47) 

y^^X y^^^ exp(iP-X + ip-x)(0|Tgcj(a;Q)q'/3(x^)|P, sc [ax]) , 

with the relative coordinate x = Xa — xg and X = (1 — a)xa + crxp. The Dirac 
part of the conjugate vertex functions is obtained by charge conjugation, 

f(p,P) = C {x'(-p,-P)f , (2.48) 

r{P:P) = -C {x'{-P:-P)f ■ (2.49) 

They have to be antisymmetric with respect to the interchange of the two 
quarks and this entails 

x'J;!;\p,P)=-xt\-p^P) , , ■ (2-50) 

As any two quarks in a baryon have to belong to the color antitriplet represen- 
tation, the diquark-quark vertices are proportional to the antisymmetric tensor 
^ABD with color indices A,B for the quarks, and with D labelling the color of 
the diquark. Furthermore, the scalar diquark is an antisymmetric flavor singlet 
|0) represented by {T2)ab, and the axialvector diquark is a symmetric flavor triplet 
|1, k) which can be represented by {T2Tk)ab- Here, a and h label the quark flavors 
and k the flavor of the axialvector diquark. For their structure in Dirac space, 
we maintain only the dominant components, the antisymmetric matrix (7^C) for 
the scalar diquark and the symmetric matrices (7''C) for the axialvector diquark. 
Thus, with all these indices made explicit, the vertices read"^ 

xUp) = 9s{l'CUV{p\p.P) i^^, (2.51) 

X:,{P) = 9a{YCUV{p^p.P) (2.52) 

Normalizations of the flavor and color part have been chosen such that for the 
matrix elements in the respective spaces the following holds, 

(3, /|3, J)color = \^IMN^JMN — 5^'' , (2.53) 

(O|0)flavor = iTrr^Vs^l, (2.54) 

(l,i|l,i)flavor = |Tr(T2T,)^(T2T,)=5^^ (2.55) 

"^Symbolically denoting the totality of quark indices by the same Greek letters that are used 
as their Dirac indices should not cause confusion. The particular flavor structures are tied to the 
Dirac decomposition of the diquarks, color-3 is fixed. 
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The antisymmetry condition (2.50) translates into a symmetry condition for 
the scalar function V, 

V{p\p-P) = V{p\ -p-P)l^^,_^^ . (2.56) 

For (7 7^ 1/2 and thus for p :— p\^^^.^_^^ 7^ p, it is not possible to neglect the 
p-P dependence in the vertex without violating the quark-exchange antisymmetry. 
To maintain the correct quark-exchange antisymmetry, we assume instead that the 
vertex depends on both scalars, p^ and p ■ P in a specific way. In particular, we 
assume the diquark-quark vertex is given by a function that depends only on the 
scalar 

x:= {l-2a)p-P + p^ = -{I - 2a)p- P + p^ (2.57) 

with p = {1 — a)pa — crpi3 and P{a,i3} as given above. The two scalars that may be 
constructed from the available momenta p and P which have definite symmetries 
under quark exchange are given by the two independent combinations Pa ■ Pp 
(which is essentially the same as above x) and — P^- The latter may only 
appear in odd powers which are associated with higher moments of the Bethe- 
Salpeter vertex in a polynomial expansion in — p^. Hence, these are neglected 
by setting 

V{p\p-P)^V{x) (2.58) 

which obviously satisfies the symmetry constraint given by eq. (2.56) Vcr G [0, 1]. 
For the actual calculations we will restrict ourselves to equal momentum parti- 
tioning between the quarks in the diquark correlation, a — 1/2. The scalar x then 
reduces to the square of the relative momentum, x = p^. 

The overall strength of the diquark correlations given in eqs. (2.51,2.52) is 
hidden in the "diquark-quark coupling constants" Qs and Qa- To estimate their 
values, we can resort to the canonical normalization condition for Bethe-Salpeter 
vertex functions [34] which reads in our context, 

UW)"^'-^ 4. - 2mL, (2.59) 

1 1 W)~' ^^"^^^ ^'a^^^?^)"-'/^^ fe)T = 6mL. (2.60) 

Here, {x'')t = ~P^{P^X^) is the transverse part of the vertex ■ Note that the 
pole contribution of the axialvector diquark is transverse to its total momentum, 
and the sum over the three polarization states provides an extra factor of 3 on 
the r.h.s. of eq. (2.60) as compared to eq. (2.59). With normalizations as chosen 
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in eqs. (2.51,2.52) the traces over the color and flavor parts yield no additional 
factors, see eqs. (2.53-2.55). 

In the following, we will solve the baryon bound state problem with the cou- 
pling constants gg and Qa fixed in the indicated manner, keeping in mind that 
the Bethe-Salpeter normalization condition determines these values only on-shell. 
When calculating diquark contributions to the nucleon electric form factors, we 
will see that the Ward identity provides for an indirect off-shell constraint on gs 
and Qa- However, the ratio between the constants fixed either on- or off-shell will 
be of the order of unity. 

In the model sector with scalar diquarks only, we will explore the ramifications 
of several ansdtze for the scalar functions V{x), 

Kxp = cxp{-x/Xl^p) , (2.62) 
V^gau = exp{-{x-Xof/XlJ. (2.63) 

The last form, the Gaussian form for the diquark vertices, was suggested as a result 
of a variational calculation of an approximate diquark Bethe-Salpeter equation 
in ref. [49]. The nucleon calculations of ref. [46] used this form without fixing 
the strength Qg by a condition such as eq. (2.59). It turns out that the necessary 
value for Agau to obtain a reasonable nucleon mass is about an order of magnitude 
smaller than the value given in ref. [46]. Furthermore, we will show in sect. 2.5.2 
that the results for the electric form factors in the scalar diquark sector clearly rule 
out the Gaussian and the exponential form of the diquark vertices and favor the 
dipole ansatz. 



2.3 The diquark-quark Bethe-Salpeter equation 

Equipped with the separable form of the two-quark correlations, eq. (2.43), and 
the functional form of the scalar and axialvector diquark correlations in 
eqs. (2.51,2.52), we will have a closer look on the effective Bethe-Salpeter equa- 
tion derived in section 2.1 and set it up for the nucleon and the A in tum.^ To 
complete the model definition, we have to specify the functional form of the quark 
propagator. We take it to be a free fermion propagator with a constituent mass m^, 

S{p) — —5 ^ Cip^-i ITT' — iTT'o) with C(p^, m — mq) — 1. (2.64) 

p"* -|- rriq 

^We will sum up the four nearly mass-degenerate states A++, A+, A°, A~ under the label A 
in the following. 
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This incorporates from the beginning an effective quark mass in the order of mag- 
nitude of several hundred MeV which is believed to be dynamically generated 
in QCD. The property of confinement is neither manifest in eq. (2.64) nor in the 
form of the diquark propagators, eqs. (2.45,2.46). Effective parametrizations of 
confinement in the model propagators and its influence on results will be investi- 
gated in chapter 5. 



2.3.1 Nucleoli 

The Faddeev amplitude of eq. (2.35) takes the following form, 

V^p, = xl, D {^\)^ + x^, D'^'^ {^^u)^ . (2.65) 

The necessary antisymmetrization of r^/^-y, provided by the summation over even 
permutations of the multi-indices ai3-i as in eq. (2.35), will be understood implic- 
itly from now on. The spinorial quantities 

P) u{P) := u{P) (2.66) 

describe the effective nucleon-quark-diquark vertex functions. The quantity u{P) 
is a positive-energy spinor for a spin- 1/2 particle with momentum P. 

We define Bethe-Salpeter wave functions by attaching quark and diquark legs 
to the vertex functions, 

tin) := ^4,^)) . (2.67) 

"(P.-P) {t^(p,P) = S{p,)D(pi)('f)\[p,P)- (2.68) 



Vertex and wave function $ and ^ are Dirac matrices to be constructed from the 
quark-diquark relative momentum p, the nucleon momentum P and the set of 7 
matrices. The task of their decomposition will be taken up in the next section, 
showing that there are altogether eight components that encompass all possible 
spin-orbit couplings of quark and diquark to the nucleon spin. 

Inserting all these definitions into the Bethe-Salpeter equation (2.29), we find 
the equation 
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in which (G^ '^'i) ^ {p,k, P) is the inverse of the full quark-diquark 4-point func- 
tion. It is the sum of the disconnected part and the interaction kernel, 

{Gr''^)-\p,k,P)-K''^ , (2.70) 
{27tf 6\p - k) S-\p,) D-\p,) , (2.71) 
1 ( -X'ipi) S^ig) xhvl) V3 x'^'ivl) S^i,) 

(2.72) 

For completeness we list the definitions of quark, diquark and exchange quark 
momentum, see eqs. (2.32-2.34), 

Pg = r]P + p , (2.73) 
Pd = (l-r?)P-p, (2.74) 
q = -p-k + {l-2ri)P . (2.75) 

The diquark vertices and their conjugates depend on the relative momenta between 
a spectator quark and the exchange quark, 

pi = p + k/2-{l-37])P/2 , (2.76) 
P2 = -k-p/2 + {l-3r])P/2 . (2.77) 

Physical quantities should not depend on the momentum partitioning parameter 7] 
which expresses just the invariance under reparametrizations of the relative mo- 
mentum. As we explained earlier, this follows from relativistic translation invari- 
ance for the solutions of the Bethe-Salpeter equation. This implies that for every 
solution "^{p, P; r]i) of the Bethe-Salpeter equation there exists a family of solu- 
tions of the form ^{p + {r]2 — r]i)P, P; 772). The manifestation of this invariance 
in actual numerical calculations will be investigated in section 2.5.2. 

Color and flavor factors have already been worked out in eqs. (2.69,2.70). 
Their derivation has been relocated to appendix A.2.1. 



{G^-^^)-\p,k,P) = 



2.3.2 Delta 

The A(1232) resonance is a spin-3/2, isospin-3/2 state. Therefore only axial vector 
diquark correlations can be present in the 2-quark t matrix since the flavor state of 
the A has to be fully symmetric. Accordingly its Faddeev amplitude reads 

rS37 = (^"'^na . (2.78) 
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The effective A-quark-diquark vertices are now quantities which transform as 

vector-spinors: ^'^'^{p, P)u'^{P), where u^{P) is a Rarita-Sch winger spinor de- 
scribing a free spin-3/2 particle with momentum P. We may construct Bethe- 
Salpeter wave functions by attaching legs, 

-^^^{p, P) = S{p,) D^^^'ipa) ^^'""{p, P) . (2.79) 

Both vertex and wave function ^^'^ and ^!^^^ are Dirac matrices with two un- 
contracted Lorentz indices, respectively. Their complete decomposition in Dirac 
space will be shown in the next section. 

The Bethe-Salpeter equation for the A is in compact notation 

where the inverse quark-diquark propagator in the A-channel is given by 

(gX"'^) (p, K P) = (27r)^5^(p - k) S-\p,) (D^^')-\pa) + 

x''{pl)S^{q)r{pl). (2.81) 

The necessary color algebra that enters the explicit form of does not 

differ from the nucleon case as all diquarks are in a color antitriplet state. The 
computation of the flavor factor is given in appendix A.2.2. 

2.4 Decomposition of the Faddeev amplitudes 

At this point we will solve the problem of the Dirac space decomposition of nu- 
cleon and A vertex functions with quark and diquark. The exercise seems to be 
clearly straightforward in the sense that, using a finite basis of Dirac matrices, 
tensorial (A) or vector and scalar matrices (nucleon) have to be found that yield 
proper spin and describe positive parity and energy solutions. We will go a step 
further, however, and show that representations can be found which correspond 
to a partial wave decomposition in the baryon rest frame, i.e. a decomposition 
into eigenfunctions of spin and angular momentum. Thereby we recover all non- 
relativistic basis elements that would survive a corresponding reduction. In addi- 
tion, virtual components due to the off-shell axialvector diquark and relativistic 
"lower components" will be found. The results obtained in this section are thus 
an explicit illustration of the difference between relativistic and non-relativistic 
baryon wave functions. 
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Nucleon 

In eq. (2.65) we have defined the nucleon Faddeev ampUtudes in terms of the 
effective spinors P) u{P) and P) u{P) that describe the scalar and 
the axialvector correlations in the nucleon. In the following we will describe the 
complete decomposition of the vertex function $ in Dirac space such that the 
Faddeev amplitudes describe a spin- 1/2 particle with positive energy and have 
positive parity. 

First let us define the projectors onto positive energy, A+, and for negative 
energy. A", by 




(2.82) 



The nucleon bound state mass is denoted by Mn, and in the rest frame the nu- 
cleon momentum is given by P = (0, iMn). The projectors defined above form a 
complete set (A"*" -|- A~ = 1) and thus the effective spinors can be written as 

$ u{P) = $(A+ + A-) u{P) = $ A+ u{P) . (2.83) 

In the following, we will substitute the projected vertex function for the vertex 
function itself, ^A"*" Clearly it is then an eigenfunction of the positive- 

energy projector, 

$A+ = $ . (2.84) 

Next we consider the transformation properties of the Faddeev amplitudes un- 
der parity. To describe nucleons with positive parity the Faddeev amplitude must 
fulfill 

^T^p,{q,p,P) = {rfx>{fi)-i')0,D{S>,P){ri'^{S>,P)u{P))^ + 
(7V(g)7')/3T D^^{p,P) {^'^^{p,P)u{P))^ 
= r,^^(g,p,P), (2.85) 

with p — App, P = ApP and 

A^^ = diag4 {-1,-1, -1,1}. (2.86) 

Using the transformation properties of the diquark-quark vertices and a Dirac 
spinor, 

t'xY = X', (2.87) 

7V7' = -A^X^ (2.88) 

j'u{P) = u{P), (2.89) 
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one can deduce from eq. (2.85) a condition for the vertex function 



-$'^(p,P) 



(2.90) 



We see that apart from the Dirac indices $^ transforms Uke a scalar and Uke a 
pseudovector. 

The conditions (2.84) and (2.90) greatly restrict the number of independent 
components in the vertex function. The scalar correlations are described by 
two components and the axialvector correlations by six components. 



$5(p,P) 



^ j:S,{p^p-P)S,{p,P) ^ 

1=1 

.ZMp^p-p) ^,At{p,P), 
V=i / 



(2.91) 



The scalar functions Si and Ai depend on the two possible scalars (p^ and p ■ P) 
that can be formed out of the relative momentum p and the nucleon momentum 
P. The Dirac components describing the scalar correlations may be built out of 
A+ and |^ A+ and the Dirac part of the axialvector correlations can be constructed 
using the matrices P^A+, P'';^ A+, 7^A+, 7^/^A+, p^A+ and p^'^k^. 

In table 2. 1 the Dirac components Si and are given as certain linear combi- 
nations of these matrices. First, they obey a useful trace orthogonality condition. 
Defining the adjoint vertex function in terms of conjugate Dirac components. 



$5(p,p) 
$'^(p,P) 



^ j:Si{p,P)Si{p\p-P) ^ 

1=1 

E^f(p,p)75 A{p^p-P) 
1=1 



(2.92) 



we find for them 



S^{P,P) 

A'^iP, P) 



C{S,{-p,-P)f 
-C {Ati-p^-P))^ 



(2.93) 



This identification of the conjugated components follows by charge-conjugating 
the Faddeev amplitude and using the definition of the conjugate diquark-quark 
vertices in eqs. (2.48,2.49). The trace orthogonality property takes the following 
form. 



Tr SiSj 

Tr a'im; 



2f-l 



(5,: 



Tr 5,75^^ = Tr ^(^755, = 



(2.94) 
(2.95) 
(2.96) 
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51 A+ 

52 -iPTA+ 

At P^A^ 

A'i -iP''PTA+ 

•^3 ^ It 



Table 2.1: Basic Dirac components of the nucleon vertex function. With a hat we denote 
normaUzed 4- vectors, p = In the case of the complex on-shell nucleon momentum, 

we define P = The subscript t denotes the transversal component of a vector with 
respect to the nucleon momentum P, e.g. pr = P — {p-P)P. 

As a second attribute, the Faddeev amplitude written with these Dirac compo- 
nents is given as a complete decomposition into partial waves, i.e. into compo- 
nents that are eigenstates of the spin as well as of the angular momentum operator 
in the nucleon rest frame. This will be shown in section 2.4. 1. 

Delta 

The Faddeev amplitudes for the A are described by an effective Rarita-Schwinger 

spinor ^^''u'' , see eq. (2.78). As in the nucleon case we only need to consider 
the projected tensor vertex function ^^''P'^^ P^"^ is the Rarita-Schwinger 

projector [50] onto positive-energy, spin-3/2 spinors, 

/ 1 OP'^P" i P'^'y" - P'''y'^\ 

W := A+ S"" - -7^7^^ + . =: A+A"" . (2.97) 

V 3 ^ ' 3 Mi 3 Ma / 

It obeys the constraints 

The tensor vertex function must be an eigenfunction of the projector, therefore 



(2.99) 
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Due to the constraints (2.98), T'^'' must be either proportional to S'^'' or the trans- 
verse relative momentum px = p — P{p- P) may contract with one index of the 
Rarita-Schwinger projector and therefore 

j^np ^ xi'di'P + E'^fj. , (2. 100) 

with V and S'^ being Dirac matrices yet to be determined. 

To further constrain the vertex function, let us consider the effect of the parity 
transformation onto the Faddeev amplitude, 

Pr%^{q,p,P) = {^\''{qh')^,D^^{p,P){^%^''{p,PX{P)U 

= r%^{q,p,P). (2.101) 

We need to evaluate the effect of a parity transformation on the Rarita-Schwinger 
spinor. The latter may be built as follows [51], 



P, s)^Y. ^t^-s' ''miP) s') . (2.102) 



The are spin-1 polarization vectors (m = 0, ±1) that can be written in terms 
of a 3 -dimensional spin-1 polarization vector basis as 

V 'Ma / 

In eq. (2.102), these covariant polarization vectors are combined with a Dirac 

spinor u{P) of helicity s' and appropriate Clebsch-Gordan coefficients C^^ i^, to 

give a spin-3/2 object with helicity s. From these definitions it can be seen that 
under parity the u'^ transform as follows, 

-^^^t^u\P) = -u^{P) , (2.104) 

and therefore the vertex function must obey 

P$'^'^(p, P) = 7^ {A'^''<^>p\p, P)A^'') 7^ = ^^"{p, P) (2.105) 

which translates into conditions for the unknown Dirac matrices appearing in 
eq. (2.100), 

'D'{p,P) \ ,^fV'{p,P)\ 
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^8 





Table 2.2: Basic Dirac components of the A vertex function. 

These transformation properties suggest that V is a linear combination of the Si 
and can be composed using the that were found for the nucleon. Again we 
will reshuffle the matrices to find a decomposition that is orthogonal and yields 
partial wave components: 

8 

$^'^(p,P) = J2D^{P^P■P) T^riP^n ■ (2.107) 

i=l 

The eight components are given in table 2.2. To state the orthogonality con- 
dition we write the adjoint vertex function as 

8 

$^^(p,P) = Y,^i"i.P^P)Di{p\p-P) . (2.108) 

The definition of the adjoint Dirac components is inferred from the charge-con- 
jugated Faddeev amplitude and the corresponding relations for the axialvector 
diquark vertex and the Rarita-Schwinger spinor, 

pr(p, p) = c -p)f c-^ ■ (2.109) 

Orthogonality of the basic components can now be expressed as 



TrPfPf = 4(-l)^+^(^^^' 



(2.110) 
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2.4.1 Partial wave decomposition 

In a general moving frame, all the components of the Faddeev amplitude r|^'^ 
possess just the total angular momentum, 1/2 for the nucleon and 3/2 for the A, 
as a good quantum number. A further interpretation of the components is not ob- 
vious. In the rest frame of the bound state, however, the trispinor r|^'^ can be 
written as a sum of trispinor components each possessing definite orbital angu- 
lar momentum and spin, thus allowing a direct interpretation in terms of partial 
waves. 

First we set up the general formalism and apply it to nucleon and A in turn. In 
the rest frame the Pauli-Lubanski operator for an arbitrary tri-spinor il^ap^ is given 
by 

= ^e,,■fc£^■^ (2.111) 
whose square characterizes the total angular momentum. 

The tensor C is the sum of an orbital part, L, and a spin part, S, which read 

Saa' ® Spp' ia^%Y , (2.114) 

such that C = L+S. The tensor L is proportional to the unit matrix in Dirac space. 
The definition a^^ :— —^[y^,'y'^] differs by a minus sign from its Minkowski 
counterpart. The tensors L and S are written as a sum over the respective tensors 
for each of the three constituent quarks which are labelled a — 1,2,3 and with 
respective Dirac indices aa',/3/3',77'. 

With the definition of the spin matrix = \tijkCF^^ the Pauli-Lubanski oper- 
ator reads 

{W^)aa',l3f3','YY — (^aa' ® <H) (^77' + (5'*) aa',/3/3', 77' , (2.115) 



V = (-i)e. 



ijk 



Qpk Qqk 



(2.116) 



(^aa' ® V (S%7') ' (2.117) 
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where we have already introduced the relative momentum p between quark and 
diquark and the relative momentum q within the diquark via a canonical transfor- 
mation: 

P^p^+p^+p^, p = 77(p^ +/)- (1 q = \{p^-p^)- (2.118) 

In the diquark-quark model, the quark (a) and diquark (/?7) Dirac indices of 
z/'q,^^ do separate. In the parametrization of the diquark-quark vertices 
and xjj^^(g), cf. eqs. (2.51,2.52), we retained only the dominant Dirac compo- 
nent along with a scalar function of q. Therefore no orbital angular momentum is 
carried by the diquarks, 

L'X^^(?) = L'X^^(?) = , (2.119) 

and the second term in the definition of V- in eq. (2.116) depending on q may 
safely be dropped. The operator now takes the form 

= L2 + 2X-5 + S^ (2.120) 

\} = ('2p^^-p2A,+pvA^') , (2.121) 

\ ap^ ap^ op> J 

d 

5aa' ® ® (^77' + V ® (S%7'] ) ' (2-122) 

(S^)aa',/3/3',77' = \ {^^aa' ® ^0j3' ® (^77' + 

2 5,,,® (E0/3/3'«)(S%y) . (2.123) 

For evaluating the action of the spin operator on the diquark-quark vertices 
and the following identities proved to be useful, 

S^'(7^C') + (7^C)(E^y = 0, (2.124) 
E^(7^C)(E^y - -3(7^^) , (2.126) 
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Dirac Faddeev amplitude EV EV 

com- component /(/ + 1) s(s + 1) 

ponent of of 



^51 (n (75C)^, I 








4 



P^(°) (7^C)^, 3 



^^^^^^^ ' ^ 



15 
4 



Table 2.3: Classification of the components of the nucleon Faddeev amphtude in terms 
of eigenfunctions of and in the rest frame of the bound state. 'EV abbreviates 
'eigenvalue' and q denotes an arbitrary Pauli two-component spinor. Each Faddeev com- 
ponent is to be multiplied with the diquark width function V{x = q^). NormaUzed 4- 
vectors are p = pl\p\ and normalized 3-vectors are correspondingly p = pl\p\. 
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Nucleon 

We will make use of the Bethe-Salpeter wave function ^ defined in eq. (2.68). 
The Faddeev takes the following form, if the wave function instead of the vertex 
function is used, 

{^>^{p,P)u{P))^X^{(lh,- (2.128) 

We will seek eigenfunctions of the spin and orbital angular momentum opera- 
tors and S^, where the eigenfunctions shall be of the general form as indi- 
cated above. Thus they represent components of the Faddeev amplitude with the 
spectator quark leg attached, (S'r)^^^. We note that for the wave function ^ an 
analogous decomposition as for the vertex function holds, see eq. (2.91), 



m^{p,P) 



,EMp',p.p)j,Anp,p) , 

\i=i I 



(2.129) 



A new set of scalar functions, Si and Aj, has been introduced whereas the basic 
Dirac components are the ones from table 2.1. The reason why we choose the 
wave function for further treatment is that the Faddeev amplitude written with 
the vertex function depends explicitly on the axialvector diquark propagator D*^^ 
which for the Proca form induces mixing between space components of the vertex 
function and the (virtual) time component of the axialvector diquark (and vice 
versa) in which case no partial wave decomposition can be found. This can be 
circumvented by using the wave function as in virtue of its definition the diquark 
propagator has become absorbed in ^. 

In table 2.3 we have listed all eight components together with the respective 
eigenvalues of and S^. The Faddeev components that are given in the rest 
frame of the nucleon have been built using the basic Dirac components for the 
wave function shown in table 2. 1. 

There is one s wave associated with the scalar diquark, described by «Si, and 
two s waves associated with the axialvector diquark, one coimected with its virtual 
time component, Ai, and the other described by ^3. In a non-relativistic limit 
where diquark off-shell components are neglected, only the s waves Si, A3 and 
the d wave A5 would survive. It is remarkable that the relativistic description leads 
to accompanying four p waves, the "lower components", which are expected to 
give substantial contributions to the fraction of the nucleon spin carried by orbital 
angular momentum. These p waves would not be present in a non-relativistic 
model. 
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Delta 

The Faddeev amplitude written in terms of the tensor wave function ^'^^ defined 
in eq. (2.79) reads 

S^^,{r^P+p)V^,p^{q,p, P) = P)u''{P))^ x^iq)^! ■ (2-130) 

Again, the wave function decomposition proceeds in the same manner as done for 
the A vertex function 

8 

vI/^^(p,P) = J2D^{P^P■P) 'Dr{p,P) . (2.131) 

i=l 

Applying the operators and S^, eqs. (2.121,2.123), to the Faddeev amplitude 
{Sr^)a/3j, this decomposition is seen to yield eigenfunctions of the spin and or- 
bital angular momentum operators. Their explicit form in the rest frame is shown 
in table 2.4. The strength of the partial waves is determined by the scalar functions 

Diip^p-P). 

In contrast to the nucleon, only one s wave (described by Vi) is found. Two 
d waves that survive the non-relativistic limit are given by and Py, and one d 
wave can be attributed to the virtual time component of the axialvector diquark, 
P3. All even partial waves are accompanied by relativistic "lower" components 
that could be even more important as in the nucleon case for only one dominant s 
wave is present compared to three within the nucleon. 

To summarize, the relativistic decomposition of nucleon and A quark-diquark 
wave functions yields a rich structure in terms of partial waves. Well-known prob- 
lems from certain non-relativistic quark model descriptions are avoided from the 
beginning in a relativistic treatment. First, photoinduced N—A transitions that are 
impossible in spherically symmetric non-relativistic nucleon ground states (de- 
scribed by only one s wave) will occur in our model through overlaps in the axi- 
alvector part of the respective wave functions. Additionally, photoinduced transi- 
tions from scalar to axialvector diquarks can take place, thus creating an overlap 
of the nucleon scalar diquark correlations with the A axialvector diquark correla- 
tions. Secondly, although the total baryon spin will be mainly made of the quark 
spin in the s waves, the orbital momentum of the relativistic p waves contributes 
as well which could not happen in a non-relativistic description since the p waves 
are absent there. 
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Dirac Faddeev amplitude EV EV 

com- component /(/ + 1) s{s + 1) 

ponent of of 



^3 



-^6 





a\(Tp){p(;) 












%/5 







{f{cTp){p<;)-l[a\pc;) + {<Tp),^])J^ 
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2 

6 
2 

6 
2 

6 
12 



15 
4 

15 
4 



15 
4 



15 
4 



Table 2.4: Classification of the components of the A Faddeev amplitude in terms of 
eigenfunctions of and in the rest frame of the bound state. denotes an arbitrary 
2pauiiX 3 vector -component spinor that represents the non- vanishing components of in 
the A rest frame. Each Faddeev component is to be multipUed with the diquark width 
function V{x = q^). 
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2.5 Numerical solutions 



After having described the numerical algorithm that solves the Bethe-Salpeter 
equations for nucleon and A, we will first present solutions for the nucleon with 
scalar diquark correlations only. These results, published in ref. [23], will provide 
us with constraints on the scalar functions which enter the diquark-quark vertices. 
Especially results on the neutron electric charge radius and the proton electric 
form factor signal the superiority of using a dipole form for the scalar function 
P appearing in eqs. (2.51,2.52). Although presenting form factors at this point 
is an early anticipation of results which are not derived until chapter 3, it moti- 
vates the choice of the specific diquark-quark vertex in the full scalar-axialvector 
calculations of section 2.5.3. 

2.5.1 Numerical method 

We will solve the Bethe-Salpeter equations for nucleons, eq. (2.69), and for the A, 
eq. (2.80), in the baryon rest frame as a system of equations for vertex and wave 
function. In a concise notation the system of equations reads 



The product of quark and diquark propagators is given by Gg"'^'' and K^^ is the 
quark exchange kernel. 

The procedure to solve the equations can be divided into the following steps: 

• We will project out the unknown scalar functions describing the partial 
waves from the expansions for wave and vertex function by contraction with 
appropriate Dirac components. 

• As the scalar functions depend on two Lorentz scalars, and p - P = 
iM\p\ p-P, we shift the dependence on p-P into an expansion in Chebyshev 
polynomials and project the equations onto the expansion coefficients, the 
Chebyshev moments. We will see that this expansion is close to a hyper- 
spherical expansion. 



'iP{p,P) 



cy^^ip^p) ct>{p,p) 



(2.132) 




(2.133) 



• The interaction kernel can be approximated by a finite number of Cheby- 
shev polynomials. Employing appropriate orthogonality relations, the equa- 
tions are reduced to a coupled system of one-dimensional integral equations 
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for the Chebyshev moments. The numerical effort will be considerably re- 
duced compared to the brute force method which puts wave and vertex func- 
tion on a two-dimensional grid, since only a few Chebyshev polynomials 
will be needed. 

These steps are described in more detail in the following. 

The wave function ip and the vertex function have been expanded in Dirac 
space according to eqs. (2.91,2.129) in the nucleon case and for the A according 
to eqs. (2.107,2.131). We will use the generic labels Yi (i = 1, . . . , 8) for the 
scalar functions describing the partial waves in the vertex function and likewise 
Yi for the wave function. The scalar functions depend in turn on the scalars p^, 
p ■ P, and the latter scalar reduces in the baryon rest frame to z = cos where 
6 is the angle between p and the 4-axis, see appendix A.l for our conventions 
for hyperspherical coordinates. For numerical convenience, we redefine the scalar 
functions associated with the "lower" components, 

Yi[Yi] ^ -j=L= y,[y,], {i odd) . (2.134) 

The set of scalar functions is now expanded into Chebyshev polynomials of the 
first kind [52], 

oo 

Yi{p\z) = ^i"F,"(/)T,(;^), (2.135) 

n=0 

oo 

Y,{p\z) = ^i"F«(/)r„(^), (2.136) 

n=0 

thereby defining the Chebyshev moments F^" and K". Here, we employ a con- 
venient (albeit non-standard) normalization for the zeroth Chebyshev moment by 
setting To = 1/^2. 

In the first step, we project the left hand sides of eqs. (2.132,2.133) onto 
the scalar functions Y^ and Yi using the orthogonality relations (2.94-2.96) and 
(2.110), 

Y,{p',z) = {g',f{p\z)Y,{p\z) (2.137) 
Y.ip^z) = J-0-^(H'r(p^,k^y',z,z')Y,{k^z'). (2.138) 

The new propagator matrix (qq)''^ and the modified kernel (H'Y^ are the result of 
taking Dirac traces after multiplication with the respective adjoint Dirac compo- 
nents. As a consequence, the kernel {H'Y^ depends on the possible scalar products 
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between the vectors k, p, P which are can be expressed as k'^, p^, z' — k-P, z — pP 
and y' = kp. 

Using the orthogonality relation of the Chebyshev polynomials. 



/_ 



■™W..^|.„„, (2.139) 

the system of equations for the scalar quantities Yi and Yi is projected onto their 
Chebyshev moments, 

mp') = f dzha'oriP^z) "Y^^^-nYp^p^ f-^^^ (2.140) 

yrip') = j\ ^ / ^H'r^p\ k\ y') x 

i''-"'YJ'{k^) " y ^ "^^" "^ . (2.141) 

n=0 ~ 

Of course, in the practical calculation the Chebyshev expansion can be taken into 
account only up to a m^gx moment for the Y^ and up to a n^^^^ moment for the 
Y^. The propagator matrix {g^Y^ can be expanded into Chebyshev polynomials 
as well and thereby results an expression from eq. (2. 140) that involves an integral 
over three Chebyshev polynomials. The addition theorem 

Tn{z)T^{z) = \{T\„\{z)+Tr,+m{z)) (n, m ^ 0) , (2.142) 
To{z)T^{z) = ^Tr,{z) (2.143) 

simplifies this expression to a form where a subsequent application of the orthog- 
onality relation (2.139) reduces it to a matrix equation, 

^'r(p') = (^o)'^'""^(p') yr(p') ■ (2.144) 

The elements of the propagator matrix {goY^'™" are all real due to the explicit 
phase factor i" in the Chebyshev expansions. 

Let us now focus attention on eq. (2.141), the 4-dimensional integral over the 
Bethe-Salpeter kernel, 

—^d\k\ J ^ VT^dz' j ^dzj ^dy' {H'f{p\ k\ y\ z, z') 

X . (2.145) 

n=0 ^ 
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We will take care of the angular integrations over z and z' next. To this end 
we approximate the Bethe-Salpeter kernel with Chebyshev polynomials in the 
following manner, 

{H"f = {\-^^){H'f = E c''''' T,-i{z)Tt-i{z') . (2.146) 

s=l t=l 

The coefficients c*^ ** are not the Chebyshev moments of {H"y^ but rather 

n n mmax+1 Hmax +1 

X {H"r{p^k\y',zu,zl) . (2.147) 

The 2;„ [4] are the zero5 of the Chebyshev polynomial T„i^^^+iiz) [Tn,,,,^+i{z')]. 
Eqs. (2.146,2.147) are an example of the so-called Chebyshev approximation of a 
function which we have applied to both the arguments z and z'. As is known from 
the literature [53], the Chebyshev approximation of a function is very close to the 
approximation by the minimax polynomial, which (among all polynomials of the 
same degree) has the smallest maximum deviation from the true function. 

When inserting the kernel approximation back into eq. (2. 145), we see that the 
orthogonality relation (2.139) will take care of the integrations over z and z'. We 
find 



■n—m 



X T„(2„)T„(4) (J^'')«(p^P,!/',^.,4) ■ (2-148) 

Here the sum runs also over the label j and the Chebyshev moment label n. 

Finally we have succeeded to transform the original 4-dimensional integral 
equation into a system of coupled one-dimensional equations. In summary, the 
system reads, 

ynp") = (5o)^^'"'"(/) i"r(p') , (2.149) 

/■oo 

y.'"(p2) ^ / d|A;| yf(A;2) , (2.150) 
Jo 

with the definition 

1 7 3 /"l "^maxH"! ^max~l~l 

X T^{zu)T^{z',) {H"f{p'', k\ y\ 4) . (2.151) 
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The momentum k is discretized on a mesh with n\k\ points, with typically n\h\ — 
20, ... , 50, and the integration is performed as a Gaussian quadrature. Then the 
problem is equivalent to tuning the lowest eigenvalue of a matrix equation to one 
by readjusting the parameters of the model. 

We have introduced the numerical method from a merely technical point of 
view. Combining the Chebyshev expansion of the scalar functions describing the 
partial waves with the Chebyshev approximation of the quark exchange kernel 
replaces the angular integrations over z = cos 6 and z' = cos ©' by sums that 
are truncated by mmax and nmax, the highest Chebyshev moments of the vertex 
and wave function partial waves, respectively. The question is whether there is 
another, more intuitive argument that supports a quick convergence of the solution 
with increasing mmax [^max] • 

There exists a motivation using this kind of expansion that is related to sym- 
metries in a simpler dynamical system. The numerical solution of the ladder ap- 
proximated Bethe-Salpeter equation in the massless and massive Wick-Cutkosky 
model has been studied extensively in the literature, for a review see [54]. In this 
model, two massive mesons interact by exchanging a third (massless or massive) 
meson. In the case of a massless exchange meson, the Bethe-Salpeter equation 
exhibits an 0(4) symmetry [55]. By expanding the wave functions into hyper- 
spherical harmonics ynim (cf- appendix A.l), the solutions can be classified using 
the quantum numbers n (four-dimensional angular momentum) and I, m (three- 
dimensional orbital angular momentum and its third component). The zero orbital 
angular momentum states of the Wick-Cutkosky model (I — 0) are essentially 
given by 

^wick ^ ^wick(^2) ^ (2.152) 

with the Un (z) being Chebyshev polynomials of the second kind. They are related 

to the Tn{z) by 

Tn(z) = Un(z) - zUn-i(z) , (u > 0) . (2.153) 

Thus an expansion of wave functions into the T„'s is closely related to the ex- 
pansion into hyperspherical harmonics. Furthermore, the authors of ref. [56] have 
shown that in the case of a massive exchange meson, where the 0(4) symme- 
try is broken, the expansion into hyperspherical harmonics converges amazingly 
quickly although n ceases to be a good quantum number. Typically rimax < 6 was 
sufficient to obtain stable solutions, except for the very weak binding regime. 

After this short digression let us return to the Bethe-Salpeter equations in the 
diquark-quark model. The Bethe-Salpeter equations treated here are ladder equa- 
tions as well, though complicated by the mixing of components with differing 
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angular momentum. Nevertheless, the expansion into spinor hyperspherical har- 
monics in the model sector with pointlike scalar diquarks only has been shown to 
work well in ref. [57]. This finding has been confirmed and extended to point- 
like axialvector diquarks in refs. [58, 22]. Thus we can conclude that the 0(4) 
symmetry is also in the diquark-quark model approximately valid. 

Let us point out the kinematical situations in which there might arise problems. 
Due to the singularities in the propagator denominators of Gg"'^'', the allowed 
range of the momentum partitioning parameter 77 in the nucleon Bethe-Salpeter 
equation is restricted to 



if max > fT^sc > fnq is assumed. Singularities in the exchange quark propaga- 
tor and the diquark vertices employing n-pole scalar functions V, see eq. (2.61), 
impose the additional bounds 



Note that these bounds on the value of rj arise in the practical calculations when 
performed as outlined above. In principle, r] could be chosen arbitrarily between 
and 1, but beyond these bounds the connection of the respective Bethe-Salpeter 
equation in Minkowski space and the present one in Euclidean space is no longer 
given by a simple Wick rotation, but the Euclidean space equation picks up addi- 
tional residue terms. 

If rj is chosen to be close to one of the boundaries given in eqs. (2.154-2.156) 
then the propagators vary strongly due to the vicinity of the poles and their Cheby- 
shev expansion, employed to derive eq. (2. 144), will converge slowly. Thus, the 
number of Chebyshev moments for the partial waves in the wave function ip, n^g,^, 
needs to be increased in order to achieve sufficient accuracy for the vertex func- 
tion. We will also see in the numerical solutions that the wave function expansion 
itself will converge more slowly under these circumstances. 

2.5.2 Solutions I: The scalar diquark sector 

In this subsection, we will present results in the model sector with scalar diquarks 
only, summarizing the work described in ref. [23]. We recall the Faddeev ampli- 
tude for the nucleon, cf. eq. (2.128) and table 2.3, that has been decomposed into 



77 e [1 - rrisc/Mn, mq/Mn] , 



(2.154) 



77 e [1(1 -mg/M„), i(l + m,/M„)] , 
T] e [|(1-2A„/M„), |(1 + 2A„/M„] . 



(2.155) 
(2.156) 



^We will encounter the problem of relating Minkowski and Euclidean vertex and wave func- 
tions again in section 3.2.1 when calculating current matrix elements. 
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Fixed width I |p|=o.2M„ = 1/2 



V{x) 


[M„] 


rriq [M„] 


A[M„] 


9s 


n = 1 


0.7 


U.OoD 


U.ioz 


1 1 T 1 


n = 2 


0.7 


U.uzU 


U.zy4 


Q1 TO 

yl./y 


n = 4 


0.7 


U.UUJ 


U.'tJO 


o J / 


n = 6 


0.7 


0.600 


0.574 


84.37 


n = 8 


0.7 


0.598 


0.671 


83.76 


exp 


0.7 


0.593 


0.246 


82.16 


gau 


0.7 


0.572 


0.238 


71.47 


Fixed masses 








n = 1 


0.7 


0.62 


0.113 


155.8 


n = 2 


0.7 


0.62 


0.294 


91.79 


n = 4 


0.7 


0.62 


0.495 


81.08 


n = 6 


0.7 


0.62 


0.637 


78.61 


exp 


0.7 


0.62 


0.283 


74.71 



Table 2.5: Summary of parameters used for the various scalar functions V{x) in the 
diquark-quark vertices, cf. eqs. (2.61-2.63). 

partial waves. Using scalar diquarks only, it reads in the rest frame 

(^r)„^, = (^^^^^^^ J {g,V) (75C)^^ . (2.157) 

The two scalar functions Si and §2 (related to the wave function) describe up- 
per and lower component of the effective nucleon spinor. Likewise only the two 
functions 5*1 and 5*2 enter the decomposition of the vertex function. 

We will start out with investigations of the various forms of the scalar function 
V{x) that parametrizes the diquark-quark vertex, see eqs. (2.61-2.63). Further- 
more we pick (7 = 1/2 for the momentum distributing parameter in the diquark- 
quark vertex. 

We have chosen to consider two cases: First we fix the the ratios of quark 
and scalar diquark mass to the nucleon mass (mg/M„ and nis^/Mn) and tune the 
width parameter A in the scalar function V{x) until the normalization gg, cal- 
culated using eq. (2.59), provides enough binding to reach the physical nucleon 
mass. This corresponds to the constraint that iterating the system of equations 
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77 = 0.4 





71 

' ''max ^ 


L 


4 





cS 


^max 















0.7005 


0.9551 


0.9731 


0.9738 




2 


0.7806 




0.9967 


0.9979 




4 


0.7844 


0.9922 


0.9993 


0.9999 




6 


0.7846 


0.9928 


0.9999 


1.0000 










T] = 0.31 









0.6517 


0.8379 


0.9285 


0.9431 


0.9449 


2 


0.8755 


0.8988 


0.9743 


0.9940 


0.9965 


4 


0.9326 


0.9614 


0.9904 


0.9977 


0.9996 


6 


0.9406 


0.9714 


0.9978 


0.9992 


0.9998 


8 


0.9416 


0.9726 


0.9992 


0.9998 


0.9999 



Table 2.6: Convergence of the ground state eigenvalue of the Bethe-Salpeter equation in 
terms of mmax and rimax, shown for two values of the momentum distributing parameter 
77. For the kernel integration the Gauss quadrature with Uyi = 20 grid points was used, cf. 
eq. (2. 15 1). The momentum integration in the iteration was performed with =20 grid 
points, cf. eq. (2. 150). Higher accuracy in these integrations does not affect the eigenvalue 
in the precision given here. 



(2.149,2.150) yields an eigenvalue 1. The results for the ratios rUg/Mn = 0.62 
and nisc/Mn — 0.7 are given in the upper half of table 2.5. With this choice and 
the choice of a dipole form for V{x) we find a good description of the nucleon 

electric form factors, see below. In the second case, we have kept fixed the diquark 
mass {rrisc/Mn = 0.7 as before) and vary the quark mass and the width parameter 
A until a solution of the Bethe-Salpeter equation was found with the property 

^?(p')||p|=o.2M„ = l/2. (2.158) 

This condition fixes the width of the zeroth Chebyshev moment of the s wave 
amplitude 5"! which we normalize by setting 5'°(0) = 1. We can regard this as a 
condition on the extension of the whole nucleon-quark-diquark vertex in momen- 
tum space since the zeroth moment of 5*1 dominates in the expansion of the nu- 
cleon vertex function. The resulting parameter sets that solve the Bethe-Salpeter 
equation are displayed in the lower half of table 2.5. 

All of the parameter sets suffer from a general defect that results from the 
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Chebyshev Moments Chebyshev Moments 

Wave Function, S,, t|=0.4 Wave Function, S,, t|=0.31 




0.1 0.2 0.3 0.4 0.5 0.6 0.05 0.1 0.15 0.2 



|p|/M„ |p|/M„ 



Figure 2.6: The Chebyshev moments of the dominating scalar function in the ex- 
pansion of the wave function. In the left panel, the moments for the choice rj = 0.4 are 
shown. For the respective moments depicted in the right panel, the employed value of 
rj = 0.31 is close to the diquark propagator poles. 

restriction to scalar diquarks only. The constituent quark mass rriq has to be cho- 
sen rather large (corresponding to a strongly bound scalar diquark) such that the 
diquark normalization Qs provides enough binding energy for the nucleon. Fur- 
thermore, the diquark-quark vertices are rather narrow in momentum space as can 
be seen from the values for A. In configuration space^, the extension of the scalar 
diquark (measured by the FWHM value for the scalar function V) is roughly 1 
fm and thus the diquark is as broad as the nucleon appears to be when probed 
with electrons or pions [59] ! The inclusion of the axialvector diquark provides 
another attractive channel, therefore the scalar diquark normalization Qs could be 
adjusted to smaller values and consequently the constituent quark mass will drop. 
The results to be shown in section 2.5.3 indeed confirm this argument. 

We turn now to the convergence properties of the Chebyshev expansion, pick- 
ing as an example the set of parameters employing a dipole form (n — 2) for 
the scalar function V{x). Here the pole in the diquark propagator and the poles 
appearing in V{x) put the restriction r] G [0.3,0.52], cf. eqs. (2.154,2.156). We 
exemplify our results with the choice r] = 0.4 (far away from any poles) and 
7] = 0.31 (close to the diquark propagator pole). First we investigate the conver- 
gence of the eigenvalue (that is obtained in the numerical iteration of the Bethe- 



^We refer here to the Fourier transform of the relative momentum between the quarks within 
the diquark that corresponds to the distance between the two quarks. 
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Salpeter equation) in terms of the highest Chebyshev moment mmax employed in 
the expansion of the vertex function and rimax for the wave function. The results 
of table 2.6 show that for the "safe" value of 77 = 0.4 seven Chebyshev moments 
for both wave and vertex function are sufficient to determine the eigenvalue to a 
precision below 10~^. For the choice of r] near the diquark propagator pole two 
more Chebyshev moments are needed to arrive at the same precision. Further- 
more the results demonstrate that in order to obtain an ry-independent eigenvalue, 
the full dependence of the scalar functions Si, Si on the angle z = cos has to be 
taken into account, i.e. r^-invariance is not guaranteed for an arbitrary truncation 
in the Chebyshev expansion. 

In the determination of observables the explicit solutions for wave or vertex 
function enter. We will therefore have a look at their convergence properties as 
well. The Chebyshev moments of the scalar functions 5*1 and S2 building the 
vertex function show a very similar convergence pattern for the different choices 
of 77. Each subsequent moment drops in its magnitude by a factor of 7 . . . 10 
as compared to the magnitude of the previous one. The situation is completely 
different for the Chebyshev moments of 5*1 and S2, making up the wave function. 
This is depicted in figure 2.6. Whereas the convergence for the choice oir] = 0.4 
is still reasonably rapid (a factor of 10 between the subsequent even Chebyshev 
momenta), the moments of Si for the choice 77 = 0.31 become squeezed in mo- 
mentum space and only converge slowly with increasing Chebyshev order. 

In the calculation of form factors we will need both wave and vertex function, 
depending on the diagrams to compute. The slow convergence of the wave func- 
tion expansion will put certain restrictions on the numerical accuracy whenever 
the respective scalar functions are employed. 

Electromagnetic Properties 

In the following we will quote from ref. [23] some results for electric form factors, 
electric radii and magnetic moments which have been obtained with the parameter 
sets of table 2.5. The definitions for the current operator, the form factors and radii 
are collected in chapter 3. The necessary formalism for calculating the appropriate 
current matrix elements is also presented there, along with the proof of current 
conservation. Presenting the results at this point will exhibit the shortcomings of 
the scalar diquark sector and furthermore motivates the choice of the dipole form 
for V{x) in the calculations with scalar and axialvector diquarks. 

In table 2.7 we have compiled the results for the proton and neutron charge 
radii as well as for the nucleon magnetic moments. A first glance at the proton 
charge radius confirms that fixing the width of the dominant scalar function S^ of 
the nucleon-quark-diquark vertex function fixes roughly this radius as well. But 
whereas all parameter sets show for the proton charge radius a maximal deviation 
of about 15 % to the experimental value, the results for the neutron charge radius 
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Fomi of diquark Vp (fm) (fm^) Hp ji. 
vertex V{x) 



fixed 5'i_o-width : 


n — 


1 


0.78 


-0.17 


0.95 


-0.80 


n = 


2 


0.82 


-0.14 


1.09 


-0.93 




n = 


4 


0.84 


-0.12 


1.13 


-0.97 




exp 




0.83 


-0.04 


1.16 


-1.00 




gauss 


0.92 


0.01 


1.22 


-1.07 


fixed masses: 


n = 


1 


0.97 


-0.24 


0.98 


-0.86 




n — 


2 


0.82 


-0.14 


1.09 


-0.93 




n — 


4 


0.75 


-0.03 


1.12 


-0.95 




exp 




0.73 


-0.01 


1.13 


-0.96 


experiment 






0.84 


-0.11 


2.79 


-1.91 



Table 2.7: Electric radii of proton/neutron and the nucleon magnetic moments for the 
parameter sets having either the S'^-width or the quark mass fixed. 

differ drastically. The dipole and the quadrupole form for V{x) are here closest to 
the observed radius. 

We have plotted the nucleon electric form factors up to a momentum transfer 
of Q'^ — 2.25 GeV^ in figure 2.7. The proton form factor (left panel) has been 
normalized to the phenomenological dipole fit 

Gt"'" = ^ o , (2.159) 

(l + g7(0.71GeV^))' 

which represents for momentum transfer below 1 GeV^ a quite accurate descrip- 
tion of the experimentally measured Ge [60]. Here the results for the parameter 
set employing the dipole shape for V{x) follows the phenomenological dipole 
closest as compared to the other curves. Therefore we employ the dipole in the 
diquark vertex for further calculations. Regarding the overall shape of the neutron 
form factor, only the Gaussian parametrization for V{x) produces results grossly 
deviating from the scarce experimental information one has about this quantity. A 
discussion of the available experimental results is postponed until section 3.2.2. 

Returning to the discussion of the results in table 2.7, we see that the pro- 
ton and neutron magnetic moments are {i) much too small and {ii) their ratio is 
far from the experimentally measured value. Both defects are known from ex- 
ploratory studies within other forms of quark-scalar diquark models [62] and can 
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Figure 2.7: Proton (left panel) and neutron (right panel) electric form factors for the 
various parameter sets with fixed width of the zeroth Chebyshev moment 5i. The experi- 
mental data for the neutron is taken form ref. [61]. 

be attributed to the oversimplification of the two-quark correlations to just the 
propagation of a spinless scalar diquark that cannot exhibit a spin flip when struck 
by a photon. The problem is furthermore enhanced by the comparatively large val- 
ues of the constituent quark mass in all parameter sets. Remember that the main 
indication for a quark mass scale of around 0.3 GeV comes from a non-relativistic 
estimate of the nucleon magnetic moments. According to this argument, the nu- 
cleon magnetic moments are proportional to Mn/irig and thus quark masses being 
larger than 0.5 GeV are expected to yield magnetic moments that are far too small. 

Although the results on the electromagnetic properties have been reported at 
this stage without derivation, we regard the previously unspecified form of the 
diquark vertex as being fixed now by the dipole form. With the ingredients of the 
model being complete, we present in the next subsection solutions for nucleon 
and A using the scalar and the axialvector diquark channel. 



2.5.3 Solutions II: Scalar and axialvector diquarks 

Including axialvector diquarks, the Bethe-Salpeter equations can now be solved 
for nucleon and A. The full decomposition of wave and vertex functions, shown 
in table 2.3 for the nucleon and in table 2.4 for the A, has to be taken into account 
and the equations are solved for the eight scalar functions that describe vertex 
and wave function respectively. Altogether there are four parameters: the quark 
mass TUq, the diquark masses nisc and max, and the width A of the dipole-shaped 
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Set 


rriq 


rrisc 


miax 


A 


9s 


Qa 


Ma 






[GeV] 


[GeV] 


[GeV] 


[GeV] 






[GeV] 


[GeV] 


I 


0.36 


0.625 


0.684 


0.95 


9.29 


6.97 


1.007 


0.939 


n 


0.425 


0.598 


0.831 


0.53 


22.10 


6.37 


1.232 


0.939 



Table 2.8: The two parameter sets of the full model together with the values of diquark 
normalizations (couplings) and the bound masses Mn, Ma that arise for these sets. 

diquark-quark vertex. 

The solutions in the scalar diquark sector suffered from the large constituent 
quark masses that set a wrong scale within the nucleon, as seen by the results 
for the magnetic moments. The inclusion of the attractive axial vector diquark 
channel now allows for solutions with smaller quark masses and smaller binding 
energy for the scalar diquark. The diquark normalization (or coupling) constant gs 
drops for this reason and therefore the nucleon binding energy becomes smaller in 
the scalar channel, but the additional binding provided by the axialvector diquark 
compensates this effect to yield the physical nucleon mass. 

In the calculations presented here and in the following chapters we shall il- 
lustrate the consequences of the model assumptions with two different parameter 
sets which emphasize slightly different aspects. For Set I, we employ a constituent 
quark mass of niq — 0.36 GeV which is close to the values commonly used by 
non- or semi-relativistic constituent quark models. Due to the free-particle poles 
in the bare quark and diquark propagators, the axialvector diquark mass is below 
0.72 GeV and the delta mass below 1.08 GeV. On the other hand, nucleon and 
delta masses are fitted by Set II, i.e. the parameter space is constrained by these 
two masses. In particular, this implies ruq > 0.41 GeV. Both parameter sets to- 
gether with the corresponding values resulting for the effective diquark couplings 
and baryon masses are given in table 2.8. 

Two differences between the two sets are important in the following: The 
strength of the axialvector correlations within the nucleon is rather weak for Set 
n, as the ratio of axialvector to scalar diquark normalization (coupling), ga/gs — 
0.29, is much smaller than for Set I where ga/gs — 0.75. We will quantify the 
correlation strength more precisely when assessing the contributions of scalar and 
axialvector correlations to the normalization of the nucleon, cf. section 3.2.2. Sec- 
ondly, the 15 % difference between the constituent quark masses of the two sets 
will have an influence on the results for magnetic moments that is even stronger 
than could be expected from the non-relativistic argument. 
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Figure 2.8: The leading Chebyshev moments of the functions Si, Ai and A3 related to 
the nucleon s waves. All functions are normaUzed by the condition 5° = 1. 



Although the wave and the vertex function are no physical observables they 
do enter observable matrix elements via Mandelstam's formalism [63] (to be de- 
scribed in the next chapter) and therefore the strengths of the single components 
give a hint on their effect on observables. We have plotted the leading Chebyshev 
moments of the scalar functions describing the nucleon s waves in figure 2.8. 
These are and ^43, describing the s waves for scalar and axial vector diquark 
and Al that is connected with the virtual time component of the latter. We see for 
both sets that the functions are suppressed by a factor of 10^ compared to the 
dominating S^. The strength of the other s wave associated with the axialvector 
diquark is roughly proportional to the ratio ga/ 9$ for the respective parameter set. 

In appendix A. 3 the nucleon partial waves for Set I and the A partial waves for 
Set n (here the A bound state mass corresponds to the physical one) are shown 
and discussed. In summary, all the remarks from the previous subsection on the 
numerical accuracy and the convergence properties for wave and vertex functions 
apply also in the case of the full model. Especially for Set II, the A is an extremely 
weakly bound state and thus the Chebyshev expansion for the wave function con- 
verges quite slowly. 

Before we turn to the detailed discussion of the nucleon electromagnetic prop- 
erties we will sum up the results obtained so far. From a physical point of view, the 
parameters of the model could be chosen more intuitively after having included 
the axialvector diquark sector of the model. This especially applies to the value 
of the constituent quark mass which appeared to be rather large in the scalar di- 



2.5 Numerical solutions 



61 



quark sector. Still, the inclusion of the A resonance leads to rather large quark 
and axial vector diquark masses and consequently small axial vector correlations 
in the nucleon. This is an artefact of employing free quark and diquark propa- 
gators, therefore we will discuss appropriate modifications which mimic confine- 
ment later in chapter 5. 

From a technical point of view, we have succeeded to transform the full 4- 
dimensional problem into a manageable form. The chosen expansion in Dirac 
space into partial waves has been done covariantly and the results for the partial 
wave strengths identify the leading and subleading components. The Chebyshev 
expansion for the scalar functions which describe the partial wave strengths has 
proved to be quite efficient and accurate, furthermore the expansion has been done 
for a Lorentz invariant variable, z ~ p-P, with the consequence that the boost of 
the wave and vertex function solutions is rendered feasible. A similar method for 
treating mesons in a Bethe-Salpeter approach has been put forward in refs. [64, 65] 
where the efficacy of this technique has been demonstrated for the solution of the 
Bethe-Salpeter equations as well as in calculations of mesonic observables. 

As a conclusion, there is no technical reason why one should stick to some 
kind of non-relativistic approximation in solving a ladder Bethe-Salpeter equa- 
tion. Such an approximation could only be justified in a very weak binding situa- 
tion and results for physical observables should depend in a controlled manner on 
the non-relativistic truncation. We investigate this issue in chapter 6 for the case of 
the diquark-quark model where we compare results for vertex function solutions 
and observables obtained (i) with the exact treatment and {ii) in the Salpeter ap- 
proximation [66], a popular semi-relativistic approximation of the Bethe-Salpeter 
equation. We find that the Salpeter approximation strongly violates the approx- 
imate 0(4) symmetry present in the solutions of the full Bethe-Salpeter equa- 
tion and observables obtained with both methods deviate in a rather unpredictable 
manner. 
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Chapter 3 

Electromagnetic Form Factors 



Electromagnetic form factors describe one of the most basic dynamic observables 
of a composite system: the elastic scattering of an electron on the composite tar- 
get. Electron and target interact with each other by exchanging quanta of the 
electromagnetic field. For the target being a nucleon, the scattering is from a clas- 
sical point of view a function of the charge and magnetization distributions within 
the target which interact with charge and spin of the electron, respectively. The 
Sachs form factors Ge and Gm, to be introduced in eqs. (3.4,3.5), describe in the 
Breit system exactly the Fourier transforms of these distributions [67]. In the usual 
field theoretical treatment it is justified to assume that the scattering is mediated 
by one virtual photon (due to the smallness of the coupling constant in quantum 
electrodynamics), thereby the information on the electromagnetic structure of the 
probe is contained in the following matrix element, defined in momentum space, 

{J^){Q) = (target; Pf\ |target; P^{Q) , (3.1) 

where initial and final target state possess on-shell four- momenta Pj = P — Q/2 
and Pf = P + Q/2, respectively. The current operator is obtained from a 
suitable model Lagrangian where it is the coefficient of the linear term in A'*, the 
electromagnetic vector potential. For nucleon targets the matrix element is written 
in the most general form as 



(j'') = fi„,(p,) 



u.JPA . (3.2) 



Here, J^i and are the Dirac charge and the Pauli anomalous magnetic form 
factor, respectively. The constant k, describes the anomalous magnetic moment. 
No further independent structures exist if the current is assumed to transform like 
a vector under parity transformations and if current conservation, 

{J") = , (3.3) 
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is to hold [68]. The Sachs forai factors Ge and Cm are introduced as 

Ge = J'i-^,kT2, (3.4) 
Gm = J^i + kJ^2. (3.5) 

By definition, the values of and Gm at zero momentum transfer are the nu- 
cleon charge and magnetic moments. 

It is convenient in the following to introduce matrix elements which are them- 
selves matrix valued due to summations over initial and final spin, 

{Pf\ := u^^iPf) (J^) u,XP^) , (3.6) 

to remove the nucleon spinors in the parametrization of (J^). Using the Gordon 
identities [68], this matrix element is written in terms of and Gm as 



{Pf\ J'^ |P,) = -iA+(P/) 



rGM + iM—{GE-GM) 



A+iPi) . (3.7) 



The form factors can be extracted from this expression by taking traces as follows, 

Ge{Q') = ^Tr{Pf\J^\P,)P^, (3.8) 

GMm = -^^T{Pf\Jnp^)ir)T, {{r)T = i'-p'if) .0.9) 

This little introductory exercise provided us with the necessary formulae to 
obtain the form factors from the nucleon current matrix element. The more formi- 
dable task is to construct a suitable current operator for the diquark-quark model 
and to give a prescription how to obtain nucleon matrix elements from the current 
operator and the solutions of the Bethe-Salpeter equation. The current operator 
should respect gauge invariance and consequently the condition (3.3), i.e. current 
conservation, should hold. 

The following section aims at constructing the current operator. This amounts 
to clarify how the photon couples to the full quark-diquark propagator G'^~^'^ given 
in eq. (2.70), the 4-point function for quark and diquark in the channel with nu- 
cleon quantum numbers. The resulting 5 -point function has to fulfill a Ward- 
Takahashi identity which in turn guarantees that for the current matrix elements 
current conservation holds. Thereby we find non-trivial irreducible couplings of 
the photon to the diquark-quark vertices x' and the seagull vertices, besides 
the expected photon couplings to quark and diquark. In total we will find that the 
nucleon current operator consists of 
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• the photon coupUng to quark and diquark (impulse approximation), 

• the photon-mediated transitions from scalar to axialvector diquarks and the 
reverse transition (extended impulse approximation) and 

• contributions from the quark exchange kernel which contains the couplings 
to the exchange quark and the seagull couplings to the diquark-quark ver- 
tices. 

Unknown constants like the anomalous magnetic moment of the axialvector di- 
quark and the strength of the scalar-axialvector transitions will be calculated by 
taking into account the diquark substructure in a first approximation which does 
not violate gauge invariance for the nucleon current. 

The construction of the current operator by using Ward-Takahashi identities 
provides us furthermore with a proof of the equivalence of the canonical nor- 
malization condition for Bethe-Salpeter wave function with the normalization to 
correct bound state charge. The canonical normalization amounts to fixing the 
residue of the bound state pole in the full propagator C^^^"^. 

Equipped with the explicit form of the current operator we will then proceed 
by calculating the form factors for the two parameter sets introduced in section 
2.5.3. We find that current conservation is reflected in the numerical solutions 
very accurately. The electric form factors of proton and neutron are described 
very well by our solutions whereas the magnetic form factors fall short by some 
15 %, due to our simplified assumptions for the quark and diquark propagators. 
An interesting point is that the recent experimental data on the ratio Ge/Gm for 
the proton [69] strongly constrains the axialvector correlations in the nucleon to 
be rather moderate. We will conclude this chapter by discussing possible further 
improvements. 

3.1 The nucleon current operator in the diquark- 
quark model 

We shall make use of the "gauging of equations" formalism outlined in refs. [70, 
71, 72]. To simplify the treatment, we will derive the current operator by assum- 
ing that the diquark quasiparticle could in principle be represented by a diquark 
creation operator S — {d\, d'^) which creates scalar and axialvector diquarks. 
The resulting expressions are equivalent to the results when constructing the cur- 
rent operator by starting from the 3-quark problem, provided a Ward identity for 
the diquark propagator holds (cf. refs. [72, 73]). 

We employ a notation for Green's functions as in section 2.1. We remind that 
their multiplication is to be understood symbolically and involves summation over 
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Dirac indices and integration over the relative momenta as in eq. (2.6). We assume 
that, in contrast to eq. (2.70), the quark-diquark propagator in the nucleon channel 

(27r)V(fc, + A;,-p,-j9,) G">-'i^(A;„ = (3.10) 

F.T. (0| Tq{xk,)d{yk,) q{x,^)d\y,,) |0) 

is given as a matrix that describes the propagation of charge eigenstate combina- 
tions of quarks u, d and scalar diquark [ud) I axialvector diquarks [mi], [ud], [dd\, 
cf. appendix A.2.1. The part described by the free quark and diquark propagators 
is of course diagonal in the charge eigenstates whereas the possible (multiple) 
quark exchange describes transitions between them. Likewise the 5-point func- 
tion is defined to be 



G^{k,Mp,,Pd) = F.T. {Q\Tq{xk,)d{yka)Q{^pM^{yva)J^{^) |0) • O-U) 

We fix momentum conservation by kg + — Pq + Pd + Q, where Q denotes the 
photon momentum. It is evident from the application of Wick's theorem to this 
equation that contains (at least) all diagrams that can be constructed from G 
by attaching a photon to all propagators (like (0| Tqq\0)) and vertices (like the 
diquark-quark vertices (0| Tdqq\0)). An equivalent definition of G^ expresses 
the procedure of attaching photons on propagators and vertices by a functional 
derivative. 



G'^{kq,kd;Pq,Pd) = 



F.T. {0\Tq{x,^)d{y,Jq{xJd\y,Jexp{-Jd'xJ^{x)A^{x)) |0)|^,^^ .(3.12) 

"Gauging" is an operation which achieves the same result by the following simple 
rules [70]: 



1. Notation: = --^^^Fa 



where 

F = F.T. (0| q{xi) . . . q{yi) |0) and 

Fa - F.T. (0| q{xi) . . . q{yi) ... exp (-/ d'^x J^'{x)A^'{x)) |0). 



2. It is a linear operation, and for a constant c its action yields zero, = 0. 

3. It fulfills Leibniz' rule, (FiFa)'^ = F2 + F^F^. 

4. Gauging the inverse, is accomplished by applying rule 2 and 3, i.e. 
{F-^FY = ^ {F-^Y = -F-^F^'F-^. 
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Besides the above mentioned contributions to from attaching photons to 

propagators and diquark-quark vertices (which are obtained by gauging G'^~'^'^), 
there exist other contributions that cannot be obtained this way: Although a transi- 
tion from scalar to axialvector diquark cannot occur, (0| Td'^ dl\0) = 0, this tran- 
sition is possible if mediated by a photon, (0| Td^ d\J^ |0)irred 7^ 0. The special 
feature of this scalar-to-axialvector transition is that it satisfies gauge invariance 
on its own, i.e. its contribution vanishes upon contraction with Q'". 

We know from gauge invariance that the five-point function must obey a 
Ward-Takahashi identity, 

q,fG'^-'"'{k, - Q, ka;p„Pd) + qdfG'^-'"'{k„ k^ - Q;p„Pd) - 
q,,G'^-'"'{kg, kd;pg + Q,Pd) - qd,G'^-'"'{k„ kd;pg,pd + Q) ■ (3.13) 

The charges ^g^igj refer to the quark charges in the final [initial] state, likewise 
qdf{di] refer to the respective diquark charges. The notation is somewhat sloppy, 
since is a matrix that describes the propagation between the different charge 
eigenstates of quark and diquark. Therefore, eq. (3.13) is to be understood /or 
each matrix element of G'^^'^'^. 

We will show that gauging G'^~^^ leads to a that fulfills the Ward-Takahashi 
identity (3.13). First, we fix some notation for the basic quark and diquark ver- 
tices. Defining —{S~^Y — be the quark-photon vertex, we find for the 
gauged quark propagator, 

S^{k,-p,):=S{k,)V>;^S{p,). (3.14) 

The Ward-Takahashi identity for it reads 

Q^S^{k,-p,) = q, {S{k, -Q)- S{p, + Q)) . (3.15) 

In the same manner we gauge the diquark propagator, given by eq. (2.67), and 
write its Ward-Takahashi identity as 

D'^ikd.Pd) := D{kd)T^^^D{pd) (3.16) 
(D{kd) \ /r^, \ (D{pa) \ 
- 1^ D'^'^'ikd)) \ r^f'/^V 1^ D^'^ipa)) 

Q^b^{kd;pd) = ^\{b{kd-Q)-b{pd + Q)), (3.17) 

\ ^ qax J 

where we have introduced the photon vertex with the scalar diquark, F^^ and the 
vertex with the axialvector diquark, V^^^. For the gauged free quark-diquark 
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propagator we obtain after applying Leibniz' rule, 

{Gl'^Y := Glj-^^ (3.18) 
rj; = D-^ + S-^ V'ag . (3.19) 

The rules 2 and 3 can now be applied to the inhomogeneous quark-diquark 
Bethe-Salpeter equation in order to find G^/ 

(^q-dq ^ (^q-dq ^ ^q-dq ^BS ^q-dq (g ^Q) 

After simple algebraic manipulations using the above equation we find 

= {Gl-'^'^Y + (gI"^'' K^^ G'^-^'^Y , (3.21) 
= ^q-dq (p/^ + [K^^Y) G'^-^'^ . (3.22) 

Since the Ward-Takahashi identities for the free diquark and quark vertices hold, 
eqs. (3.15,3.17), a Ward-Takahashi identity for Fq follows trivially. Now we con- 
sider: j/the following identity for [K^^Y holds, 

Q^^{K''''Y{Kkd■,p,,Pd) = 
qq^K^^{kq - Q, kd;Pq,Pd) + qdfK^^{kq, kd - Q]Pg,Pd) - 

qq,K''^{kq, kd;pq + Q , Pd) ' QdiK^'^ikq, kd] Pq, Pd + Q) , (3.23) 

then the identity 

(r^ + (K^'^Y) = 

q,AG''-''T\kq,kd;Pq + Q,Pd)+qdAG''-''T\kq,kd;Pq,Pd + Q) - 
qq^iG^-'^'^Y'ikq - Q, kd;pq,pd) ' ^d, (G^"'*^) kd - Q;Pq,Pd) (3.24) 

is valid, as can be inferred from eq. (3.22) and (G^-'ii)-i = (G^^'^'^y^ - K^^, 
eq. (2.70). The Ward-Takahashi identity for G^ itself, eq. (3.13), follows imme- 
diately after left and right multiplication with G'^~^'^. Indeed, eq. (3.23) which we 
will call the Ward-Takahashi identity for the quark exchange kernel is derived 
from a Ward-Takahashi identity^ for {G^"^'^K^^G2'"^Y and eqs. (3.15,3.17). 
Therefore we have established eq. (3.13). 

Matrix elements of the current between bound states are obtained from the ex- 
pression for G^, eq. (3.22). The full quark-diquark propagator, which sandwiches 

'As for G'i~^i, the kernel K^^ has not exactly the form as given in eq. (2.72) but is rather writ- 
ten as a larger matrix describing the quark exchange between all charge eigenstate combinations 
of quark and diquark. 

^Any perturbative graph with external legs and an external photon can be shown to fulfill a 
corresponding Ward-Takahashi identity. 
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the full vertex + (K^^^, can be approximated at the nucleon bound state pole 
by 

^q-dq ^ ^ - ^ + regular terms (3.25) 

as will be shown in the beginning of section 3.1.3. Inserted into eq. (3.22), we 
find for the current matrix element (the residue term of the leading pole) 

{Pf\ \Pi) = (r(; + {K^y) . (3.26) 

This expression describes how to obtain bound state matrix elements from its 
Bethe-Salpeter wave function and is usually referred to as Mandelstam's formal- 
ism [63] although in the original paper only contributions in impulse approxima- 
tion have been discussed. 

When contracting eq. (3.26) with Q^^, we can use the Ward-Takahashi identity 
(3.24). Current conservation, Q^{Pf\ |Pj) — 0, follows immediately upon 
using the Bethe-Salpeter equations 

We have arrived at an important point in this chapter. The nucleon current defined 
in eq. (3.26) above is conserved. To guarantee gauge invariance, we must include 
in the definition of the current operator the gauged Bethe-Salpeter kernel (K^^)^ 
besides the photon coupling to the free propagators Fq . The latter contribution is 
commonly referred to as the impulse approximation. Incidentally, the necessity 
of including current contributions from the Bethe-Salpeter kernel for a gauge in- 
variant matrix element was first discussed in ref. [74] in a deuteron bound state 
equation where the binding is provided by meson exchange. 

We turn now to the construction of explicit expressions for the current operator 
in impulse approximation, F^, and the gauged kernel [K^^Y. 

3.1.1 Impulse approximation 

Let us specify the electromagnetic vertices of the basic constituents, quark and 
diquark. For the quark-photon vertex, we employ the perturbative expression as 
the model propagator in eq. (2.64) is the tree-level one, 

= -iq,Y ■ (3.28) 
Of course it satisfies the corresponding Ward identity (3.15). 
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Figure 3.1: Resolved vertices: photon-scalar diquark, photon-axialvector diquark and 
anomalous scalar-axialvector diquark transition. 

The diquark propagators are the ones of a free spin-0 and spin- 1 particle, there- 
fore we also employ the corresponding free vertices 

r^c = -qscikd + PdT, (3.30) 

^ (g/3 _ ga . (3.31) 

As the spin-1 vertex is not necessarily textbook material, we refer to ref. [75] 
for its derivation. The constant k which appears in eq. (3.31) denotes a possible 
anomalous magnetic moment of the axialvector diquark. We obtain its value from 
a calculation for vanishing momentum transfer (Q"^ = 0) in which the quark sub- 
structure of the diquarks is resolved, i.e. in which a (soft) photon couples to the 
quarks within the diquarks, cf. the second graph in figure 3.1. The calculation of 
K, is provided in appendix B.1.1. Although the parameters presented in table 2.8 
differ for the two sets, especially the diquark width parameter A, we obtain in both 
cases K — 1.0 (cf. table B.l in appendix B.1.1). This might seem understandable 
from nonrelativistic intuition: the magnetic moments of two quarks with charges 
qi and q2 add up to (gi + q2)/mq, the magnetic moment of the axialvector di- 
quark is (1 + K){qi + q2)/max and if the axialvector diquark is weakly bound, 
niax — 2mq, then k, 1. 

As already mentioned, there is one contribution to the current which cannot be 
obtained by the gauging method: photon-induced anomalous^ transitions from a 
scalar to axialvector diquark and vice versa. The corresponding vertex describing 



•'if the product of the parity eigenvalues tt; of the three particles at some vertex that describes 
their interaction is (—1) we call the vertex anomalous. For the scalar-axialvector transitions, we 
have TT-y = TTsc = —TTax = 1, thus the corresponding vertex is anomalous. 
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the transition from axialvector (with index (3) to scalar diquark must have the form 

f f = e^^^\p, + k.YQ' , (3.32) 

and that for the reverse transition from an scalar to axialvector (index a) is given 

by, 

f r = e''''\pd + kdYQ^ . (3.33) 

The tensor structure of these anomalous diagrams is derived in appendix B.1.1 by 
resolving the diquarks in a way as represented by the right diagram in figure 3.1. 
The explicit factor 1/M„ was introduced to isolate a dimensionless constant Kga- 
Its value is obtained roughly as Uga — 2.1 (with the next digit depending on the 
parameter set, cf. table B.l). Since the only possible transition for nucleons in- 
volves the scalar (ud) diquark and the axialvector [ud] diquark, there is only a 
constant charge (or flavor) factor in eqs. (3.32,3.33). It is (g^ — gu)/2 = —1/2. 
We see immediately that nucleon current contributions using these vertices are 
transversal and thus do not affect current conservation. 

The current operator associated with the anomalous transitions reads 

^/T^c-ax^. = «^M] {War-'-' + {JD'""'') (3.34) 
W] ( ) / V *M •(3-35) 



As for nucleons only the transition between scalar (ud) diquark and axialvec- 
tor [ud] diquark is possible, the flavor weight factor W[ud] for the strength of 
the [ud] diquark correlations in the total axialvector correlations is included, cf. 
eqs. (A.18,A.20). 

To conclude this subsection, we write down the explicit expressions for the 
proton and neutron current matrix elements in the extended impulse approxima- 
tion, 

^ {{j^ar-^ + (j^sr-^'') , (3.36) 
m^-^) - ^ {w-^ + {j^sr-'") -(3.37) 

The superscript 'sc-sc' indicates that the current operator is to be sandwiched 
between scalar nucleon amplitudes for both the final and the initial state. Likewise 
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Figure 3.2: Impulse approximate contributions to tlie electromagnetic current. For the 
scalar-axialvector transition, a diagram analogous to the third one (with initial and final 
nucleon states interchanged) has to be computed. 

'sc-ax' denotes current operators that are sandwiched between scalar amplitudes 
in the final and axialvector amplitudes in the initial state, etc. As an example we 
give 

- / (|^^"(^/'^+ (1 - (D-'r^{ka) *'^(P,,p) .(3.39) 

The weight factors multiplying the single matrix elements arise from flavor al- 
gebra: we have chosen to define the matrix elements to be independent of the 
charge of the struck particle, cf. eq. (3.38). The nucleon is decomposed into its 
charge eigenstates of quark and diquark as in eqs. (A.18,A.20). Then the weight 
factors for the matrix elements are composed of the weight factors of the charge 
eigenstates and the charge of the struck quark or diquark. 

We note that the axialvector amplitudes contribute to the proton current only 
in combination with diquark current couplings. 

For visualization, the diagrams of the extended impulse approximation are 
shown in figure 3.2. 

3.1.2 Seagulls 

A matrix element of the quark exchange kernel (between charge eigenstates) is 
given by 

Kf! = -ha X\K, q) S^'iq) r{q,P,) (3.40) 

cf. figure 3.3. The corresponding flavor factor which is unimportant in the follow- 
ing is given by ft,a. As in the actual solutions of the nucleon Bethe-Salpeter equa- 
tion we have assumed that the diquark-quark vertices depend on just the relative 
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Figure 3.3: Left panel: The quark exchange kernel. Right panel: The seagull vertices. 
The charges belonging to the incoming and outgoing particles are enclosed in brackets, 
the other labels denote their momenta. 



momentum between the two "participating" quarks, x"(/i;,p) x"((^ — p)/2). 
Gauging the kernel leads to three terms. 



The gauged diquark-quark vertices and are called (M")'* and (M")^, re- 
spectively. These are the seagull vertices mentioned in the beginning of the chap- 
ter. 

The Ward-Takahashi identity (3.23) which leads to a conserved nucleon cur- 
rent constrains the longitudinal part of the seagull vertices. Contracting (3.41) 
with Q^, inserting on the left hand side eq. (3.23) and using the Ward-Takahashi 
identity (3.15) for the gauged exchanged quark on the right hand side, we find 

Q^{M'^r{Pa,Pp) = ?a [x"(Pa-Q,P/3)-X"(Pa,P/3)] + 

?/3 [x"(Pa,P/3 - Q) - X"(Pa,P/3)] (3.42) 
(with Pd + Q^Pa+Pp) , 

Q^iMTiPa,Pp) = -qa[riPa + Q,Pp)-riPa,Pp)]- 

\x"iPa,Pp + Q)- riPa,Pp)] (3.43) 

(with Pd-Q^Pa+Pp) ■ 

The notation is explained in the right panel of figure 3.3. Of course, such con- 
straints do not fix the seagull vertices completely. We choose a form which solves 
for this constraint and is regular in the limit Q — > 0. With p — {pa — Pp)/'^ the 
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vertices read 

{M-np,Q;q^,qp) 



{M^r{p,Q;q^,q^) 



In the limit Q ^ they simplify to 

{M-np,0;q^,qf,) = — {Qa - Qp) p"" 
(M")'^(p,0;g„,g^) = ~{qa-qp)p^ 

We add as a side remark that these seagull terms were derived a little bit dif- 
ferently in ref. [23]. There the Ward-Takahashi identity for the matrix element 
(0| TqqqqJf^ |0) has been considered. Using the separable ansatz for the two- 
quark correlations, eq. (2.43), for expressions like (0| Tqqqq\Q), the constraints 
(3.42,3.43) have been derived under the assumption that the photon-diquark ver- 
tices contain no quark loop, i.e. they are given by the perturbative expressions 
(3.30,3.31). However, in the limit Q ^ the seagulls must have the form as in 
eqs. (3.46,3.47), as here the Ward-Takahashi identities reduce to the (differential) 
Ward identities which fix the vertices completely. In this kinematical situation a 
dressing quark loop in the photon-diquark vertices must yield the diquark charge 
(for an on-shell diquark). 

Furthermore we note that seagull vertices of a similar kind were first consid- 
ered in ref. [76]. The author discusses extended meson-baryon vertices and derives 
meson-baryon-photon vertices from a Lagrangian point of view by applying the 
minimal substitution rule. 

For completeness we give the explicit expressions for the seagull and exchange 
quark contributions that are obtained when recombining the quark-diquark charge 
eigenstates to the correct nucleon flavor state. The symbolic matrix multiplica- 
tions are to be understood as accompanied by an integration over the relative 
quark-diquark momenta in the initial and final states, respectively, and a sum- 
mation over Dirac and Lorentz indices. First we write down the current term for 



i^^^ [X"(P - Q/2) - X"(P)] + 



^P-Q-Q^ 
(4p + QY 

ip-Q + Q^ 

i^p-QT 



Ap-Q-Q'' 
(4p + Q)^ 



[x"(p + Q/2)-x"(p)] , (3.44) 
[x"(p + g/2)-x"(p)] + 
[x'^(p-g/2)-x"(p)] • (3.45) 



dp^ 



(3.46) 
(3.47) 
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the proton that corresponds to the diagram in the lower right of figure 3.4, 

\ sg,as sg,aa ' ^ ^ 

K^a^ss = {M'np[,Q;q^,q,)S^{q')f{p2) 



K%,sa = -^V3(M^)'^(p;,g;3g„,2g„ + g,)5^(g')x'(P2), 



to be followed by the conjugated seagull contributions as depicted in the diagram 
to the lower left in figure 3.4, 

^jM^proton ^ _1 ( ^5 ) ^ j^^^^^ j^^^^^ ^ / ^5 ^^ (349) 

sg,as sg,aa / \ 



K^g,sa = V3 x'^iPi) S^{q) {M'r{p'„ Q- + 2g„, g„ + 2g,) , 
A^a. = xHpi) S^iq) (M")^(P2, g; 2g„ + g,, 3g,) , 

K^,aa = X^'lPi) -^^C?) iM''nP2, Q; 4?. - qu + 2g.) , 

and to be completed finally by the exchange quark contribution, cf. the upper 
diagram in figure 3.4, 

^jM)proton ^ _1 ^ ^5 W W (3-50) 

\ ex, as e2:,aa / \ / ' 

K^.,sa - (2?« + ?<^) X^Pi) S^iq) K)^ 5^(g') X'(P2) , 

K.,as = -^V3(2g„ + g,)x^(pO^^(?)(r3^^^(gOr(P2), 

i^r.,aa = -^(4?„-g.)x^(Pi)5^(?)(r9^5^(g')r(P2). 



The respective contributions to the neutron current are obtained by interchanging 
the charges of the u and d quarks, qu qd- The relative momenta in the diquark- 
quark vertices and seagulls are according to figure 3.4: Pi = kg — q, p\ = kg — q', 
P2^q' - Pq andp'2 ^q-Pq. 
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Figure 3.4: Exchange quark and seagull diagrams. 



The nucleon current operator is now complete. We recapitulate that the current 
matrix element is written as, 

(P/l \p,) - r(^,, = (r^ + r^,_,, + {K^^y) , (3.51) 

including all the discussed elements. The impulse approximation operator Fg in- 
volves the couplings to quark and diquark, cf. eqs. (3.19,3-28-3.31). The extended 
impulse approximation comprises additionally the anomalous transitions between 
scalar and axialvector diquarks. Their operator F^^-ax is given in eq. (3.35). Fi- 
nally we have discussed the coupling to the exchange quark kernel, [K^^Y, in the 
last subsection. 



3.1.3 Normalization of the Bethe-Salpeter wave function and 
the nucleon charges 

A bound state Bethe-Salpeter equation is a homogeneous integral equation, there- 
fore its solutions need subsequent normalization. The canonical normalization 
condition is derived by demanding that the bound state contributes a pole to the 
full two-body propagator with residue 1. As the nucleons are spin- 1/2 particles, 
we demand that near the nucleon pole the full quark-diquark propagator behaves 
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like a free fermion propagator, 

{G--'-U. = (-2^n)^^^l|^ • (3-53) 

Now we write the inhomogeneous Bethe-Salpeter equation (2.70) for the full 
quark-diquark propagator in the form 

^q-dq^^q-dq (^(G^-dq)-! _ ^^^5^ G<q-dq (3 54) 

and insert the pole form (3.53) on both sides to find the condition 

('/^q-dq\-i _ j^BS 1 
lim = u^^ ^ ° / , = -— — . (3.55) 

Summing over spins, applying 1' Hospital's rule and writing the result covariantly 
yields the final normalization condition 



2M„ 



Tr * 



* = 1 , (3.56) 



ip=p„ 

where Pn denotes the on-shell nucleon momentum. 

On the other hand, the electric form factor is normalized to yield the phys- 
ical nucleon charges in the soft limit, Q — > 0. Using eq. (3.8) that extracts Ge 
from the current matrix element and eq. (3.26) for the definition of the current 
operator"^ we find another condition, 

Tr ^ (r(; + [K^'^ f) ^ = g^v • (3.57) 



2M, 



In fact, the two conditions (3.56,3.57) are equivalent for the proton which has 

non-zero charge ^at. In the literature, a proof of this statement can be found [77] 
where the authors treat the case of one interaction channel and the momentum 
partitioning fixed by = 1/2. Crucial for their proof is the differential form 
of the Ward-Takahashi identity for the gauged full 2-particle propagator which 
corresponds in our case to the gauged quark-diquark propagator, eq. (3.13). The 
proof is easily accomodated for our purposes and here it comes. 

In our notation for the full quark-diquark propagator indicating all four particle 
momenta, G'^^'^°'{kq, kd;Pq,Pd), the momenta are not independent. We employ 

''As the scalar-to-axialvector transitions are transversal, they do not contribute to the nucleon 
charges. 
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therefore the equivalent notation '^^{k, p; P) with P — kq + kd— Pq+Pd and 
p, k being the relative momenta. Our usual conventions hold: 

kM] =VP + k[p] , kd[pd] = (1 - V)P - k[p] . (3.58) 

For the gauged propagator G'^{kq, kd;Pq,Pd) the four momenta are independent, 
and here we define 

Pf^Pg+Pd, Pi ^ kq + kd. (3.59) 

The Ward-Takahashi identity (3.13) can now be formulated as 

{Pf-PiYG^' = (3.60) 
^.^^^'^^(fc - (1 - r^){Pf - P,),p- P,) + qd.G^-^'^k + v{Pf - P^),P■, P^) - 
qq,G^-'''^{Kp+ (1 - r^){Pf - P,); Pf) - q<i,G^-^''{Kp - v{Pf - Pi); Pf) ■ 

As we remarked earlier, this identity holds for each matrix element of G'^"'^'^ that 
connects initial and final charge eigenstates of the quark-diquark system. We 
differentiate both sides with respect to P^ and subsequently assume the soft limit, 
p. = =: p, to find 

= -{qg. + qd.)^pj^ik,p;P)+ (3.61) 

{rjqdf - (1 - v)Qqf)—Q^{k:P; P) + (m - (1 - r])qq.) {k,p; P) . 

Since qq- + g^- = 1 (proton charge) for all the matrix elements of G*^""^^, we 
employ for the first term on the right hand side 

_^ ^^q-dq^^q-dq^-l] ^ q ^ 

_ r'q-dq ^(<^'' '^'^) V q-dq CX f,0\ 

On the other hand, we can use eq. (3.22) for on the left hand side of eq. (3.61). 
Therefore we obtain 

^q-dq ^pM _|_ (^^BS^Ai^ (^q-dq _ _Qq-dq_\ (^q-dq _|_ 

imdf - (1 - ^)%)-^^ + ivQd, - (1 - v)qgi)—g^ ■ (3.63) 

Near the bound state the full quark-diquark propagator can be written in the pole 
form of eq. (3.53). We see that the first terms on the left and right hand side 
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of eq. (3.63) have a double pole at = — whereas the remaining terms can 
only contribute a single pole. Fortunately these drop out when comparing residues 
and we are left with 

^^^p^J ' = (^0 + (K^^r) ■ (3.64) 

Putting (7 = a', performing the spin summation and contracting with P^, we 
indeed find the equality of the canonical normalization condition (3.56) and the 

charge normalization condition (3.57). Also the neutron has the correct charge: 
Although the overall normalization cannot be fixed in this case, eq. (3.61) guaran- 
tees that its electric form factor is zero in the soft limit since qq^ + qdi = ^ for all 
quark-diquark channels. 

It is possible to prove the equivalence of eqs. (3.56,3-57) without explicitly re- 
curring to the Ward-Takahashi identity (3.13). The proof is slightly more technical 
and for the scalar diquark sector of the model it is described in ref. [23]. 



3.1.4 Summary 

In conclusion we want to emphasize the power of the general Ward-Takahashi 
identity for the full quark-diquark propagator, eq. (3.13). Its implications allowed 
us to derive (i) conservation of the nucleon current, (ii) constraints on the longi- 
tudinal part of all contributing graphs, notably of the irreducible quark-diquark- 
photon vertices (seagulls) and finally (in) the equivalence of the canonical and the 
charge normalization condition. It is interesting to note that actually current con- 
servation holds for each single diagram depicted in figures 3.2 and 3.4. This has 
been proved in ref. [62] by using just the transformation properties of the respec- 
tive matrix elements under time reversal and parity. But as we know, gauge invari- 
ance provides more information than conserved current matrix elements on- shell. 
This additional information is manifest in the (off-shell) Ward-Takahashi identi- 
ties and it is certainly wrong to restrict oneself to assorted current contributions 
since at least the correct nucleon charges cannot be obtained in general. This is a 
shortcoming of a quark-diquark model investigated in refs. [78, 79]. The authors 
employ parametrized nucleon Faddeev amplitudes and calculate form factors in 
a generalized impulse approximation which, besides containing overcounted con- 
tributions [73], does not respect a Ward-Takahashi identity like eq. (3.13). The 
immediate consequence turned out to be that the correct nucleon charges have 
been obtained only for a special choice of parameters. 
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3.2 Numerical calculations 

Before we present the results for the nucleon electric and magnetic form factors, 
we will discuss shortly the numerical method. As we have calculated the Bethe- 
Salpeter wave functions ^ in the rest frame of the nucleon, we have to discuss their 
boosting to a moving frame. Thereby we will find that an accurate description of 
the boosted wave function is linked to the accuracy of its Chebyshev expansion in 
the rest frame. Since the vertex function $ show better convergence in the Cheby- 
shev expansion, it would be desirable to use $ instead of the wave function ^. We 
find that this is feasible only for the diagrams of the extended impulse approxi- 
mation but that due to their singularity structure the diagrams evaluated with the 
vertex function acquire additional residue contributions in a certain kinematical 
region. In this region the (Euclidean) vertex function $ is not connected to the 
(Euclidean) wave function ^' by simply cutting quark and diquark legs off ^. 



3.2.1 Numerical method 

The numerical computation of the form factors is done in the Breit frame where 

= (Q,o), 

Pt = {-Q/2,zu;q), (3.65) 
P; = {Q/2,tUQ), 
P^^{P, + Pf)/2 = (0,zcuq), 

withcuQ = + Q74. 

The total current operator, eq. (3.51), employs the wave functions ^ in the 
initial and final state. These are determined by the basic covariants and the corre- 
sponding scalar functions Si{p'^, z = p-Pn), Ai{p'^, z = p-Pn), cf. eq. (2.129). We 
recall that we absorbed the angular (z) dependence in the Chebyshev expansion, 
eq. (2.136). In the rest frame, P„ = (0, iM„), z is real and \z\ < 1, therefore 
the expansion is always convergent. In the Breit frame we have the real relative 
momenta in the initial state, p, and the final state, k, which are integrated over. 
Therefore the angular variables z^ = k-P and Zp = p- P are real and not larger 
than 1. The angular variables that enter the initial and final state wave functions 
as an argument are, however. 



Zi^p-Pi = 
Zf^k-Pf = 



(3.66) 
(3.67) 
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Therefore, in order to use the rest frame solutions, analytical continuation to com- 
plex values for Zi and Zf is necessary. This can be justified for the bound-state 
Bethe-Salpeter wave functions ^. These can be expressed as vacuum expecta- 
tion values of local (quark field) and almost local operators (diquark field) and we 
can resort to the domain of holomorphy of such expectation values to continue 
the relative momenta of the bound-state Bethe-Salpeter wave function ^{p, P) 
into the 4-dimensional complex Euclidean space necessary for the computation of 
Breit-frame matrix elements from rest frame nucleon wave functions. The nec- 
essary analyticity properties are manifest in the expansion in terms of Chebyshev 
polynomials with complex arguments. 

There is a practical obstacle in the outlined procedure. As the maximum of the 
squared absolute values |2;i[/]|max = 1 + Q^/(4M^) rises with increasing Q"^, we 
encounter convergence problems for the current matrix elements as the Chebyshev 
moments of higher order n become more and more important (T„ is a polynomial 
of degree n). Consider the strengths of the wave function Chebyshev moments 
presented in appendix A.3 (Set I). The fourth Chebyshev moment of the dominat- 
ing function Si is about 7 % of the zeroth moment. With Ti{z) = 8z^ — 8z'^ + 1 
one can estimate the critical for which contributions from Sf are comparable 
to the ones from 5"°, 

\Stn{zf)Us.r.\S°,To{zf)U^ Q^^lAGeV^. (3.68) 

As a consequence, the grid size in the angular integrations has to be increased be- 
yond this value of because the higher Chebyshev polynomials introduce rapid 
oscillations in the integrand. Since the Chebyshev moments of the vertex func- 
tion converge much more rapidly ( 1 5*^1 max = 0.01 |S'°|max for Set n), it would be 
desirable to use $ instead of ^ in the expressions for the current, eq. (3.51). Hav- 
ing ^ = Gq"'^'* the free quark and diquark propagators with the corresponding 
singularities are introduced in the integrand. For sufficiently small Q"^ these are 
outside the complex integration domain. For larger Q'^, these singularities en- 
ter the integration domain, therefore, a proper treatment of these singularities is 
required when they come into the domain of interest. 

For the impulse approximation contributions to the form factors we have been 
able to take the corresponding residues into account explicitly in the integration. 
Here the integration is reduced to one four-dimensional integral over either p or k 
due to the inherent 5-function of the free quark-diquark propagator Gq~'^'^. Being 
technically somewhat involved, we have placed the matter into appendix B.2. For 
these contributions, one can compare the calculations with boosted ^ and $ and 
verify numerically that they yield the same, unique results (up to several Q^). This 
is demonstrated in appendix B.2. 

For the current contributions from the gauged Bethe-Salpeter kernel, cf. figure 
3.4, two four-dimensional integrations over p and k are required. The complex 
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singularity structure of the integrand did not allow a similar treatment as described 
in appendix B.2 for the diagrams of the impulse approximation. Therefore we 
have resorted to employ the boosted wave functions As a result, the numerical 
uncertainty (for these diagrams only) exceeded the level of a few percent beyond 
momentum transfers of 2.5 GeV^. Due to this limitation (which could be avoided 
by increasing simply the used computer time) the form factor results presented in 
section 3.2.2 are restricted to the region of momentum transfers below 2.5 GeV^. 

In addition to the aforementioned complications, there is another bound on the 
value of above which the exchange and seagull diagrams cannot be evaluated. 
It is due to the singularities in the diquark-quark vertices x^^^\Pi) in the ex- 
change quark propagator. The rational n-pole forms of the diquark-quark vertices, 
^n-poie(p^) — {^n/ ('^n + P^))" for example yield the following upper bound, 

A free constituent propagator for the exchange quark gives the additional con- 
straint, 

Q'<*{-(T^,-^n) ■ (3.70) 

It turns out, however, that these bounds on are insignificant for the model 
parameters employed in the calculations described herein. 

3.2.2 Results 

For the scalar diquark sector, we have presented results for the nucleon electric 
form factors, radii and magnetic moments already in section 2.5.2. There we 
found a good description of both the proton and neutron electric form factor, 
whereas the magnetic moments turned out to be less than half of the respective 
experimental value. We attributed this failure to the absence of axialvector di- 
quarks (which carry spin and therefore contribute to the spin-flip currents) and to 
the high constituent quark mass employed for the solutions. Furthermore there is 
a substantial contribution to the electric form factors by the exchange quark and 
the seagull diagrams [23] with the latter being of the order of 30 %. We attribute 
this finding to the relatively small diquark width in momentum space and the quite 
sizeable nucleon binding energy (which is around 300 MeV for the parameter sets 
from table 2.5). 

We turn now to the form factors for the two parameter sets of table 2.8 which 
include axialvector diquark correlations. First we want to judge the strength of 
the axialvector correlations within the nucleon by evaluating the contributions 
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0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 0.2 0.4 0.6 0.8 1 1.2 



Q' [ GeV' ] q' [ GeV' ] 

Figure 3.5: Electric form factors of both nucleons for the parameter sets of table 2.8. 
Experimental data for the proton is taken from [60]. The older neutron data analysis [61] 
contains more systematic uncertainties (due to specific nucleon-nucleon potentials) than 
the more recent data from [80, 81]. 

of and to the norm integral, eq. (3.56). The strength of the axialvector 
correlations is rather weak for Set II, since the scalar diquark contributes 92 % 
to the norm integral of eq. (3.56) while the axialvector correlations and scalar- 
axialvector transition terms together give rise only to the remaining 8 % for this 
set. Scalar diquark contributions are terms ~ 'J'^il/^, axialvector contributions are 
~ ij/^'vl/'^ and the transition terms account for the two cross terms. For Set I, the 
fraction of the scalar correlations is reduced to 66 %, the axialvector correlations 
are therefore expected to influence nucleon properties more strongly for Set I than 
for Set II. 

The results for the nucleon electric form factors are presented in figure 3.5. 
For the two sets we show here the sum of all diagrams, formed by the extended 
impulse approximation and the exchange kernel graphs. In table 3.1 we have col- 
lected the corresponding radii. ^ The phenomenological dipole behavior of the 
proton G^; is well reproduced by both parameter sets. In the neutron case, the 
slope at the soft point and the general behavior are well accounted for although 
the data points of ref. [61] have to be interpreted with care since here the elec- 
tric form factor has been extracted from electron-deuteron scattering data using 
specific nucleon-nucleon potentials. Data from more recent experiments which 



radius rp corresponding to a form factor F is defined by rj^ — — 6^^|q2^o/F(0) if F{0) 
is non-zero. For the neutron electric charge radius, the quotient F{0) is simply omitted. 
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Set I Set n experiment 



(rp)ei [fm] 

(r^)ei [fm^] 

('^p)mag [fm] 

(rn)mag [fm] 



0.88 0.81 

-0.12 -0.10 

0.84 0.83 

0.84 0.83 



0.836 ± 0.013 
-0.113 ±0.007 
0.843 ± 0.013 
0.840 ± 0.042 



Table 3.1: Nucleon electric and magnetic radii for the two parameter sets compared to 
the experimental values from [82] (for (r^)ei) and [60] (for the remaining radii). 



largely avoid this model dependence are still rare. Our results lie here at the lower 
experimental bound. 

In figure 3.6 we evaluate the contributions from the single diagrams for Set 1. 
For the notation of the single current matrix elements we refer to eqs. (3.36,3-37) 
and (3.48-3.50). As can be seen, the overall behavior of the impulse approxima- 
tion diagrams is close to the total form factor. Clearly the ones sandwiched be- 
tween the scalar nucleon correlations (straight thin line) dominate over the matrix 
elements of the axial vector correlations (dashed thin line). Seagulls and exchange 
quark represent just a correction below 10 % but are nevertheless important to 
guarantee the correct behavior for Q —>^ 0. On the other side, the neutron elec- 
tric form factor is much more sensitive to the behavior of the single diagrams as 
it arises as a result of delicate cancellations. Especially we want to point to the 
contribution from the impulse approximation graphs which are proportional to the 
nucleon axialvector correlations in the initial and final states, i.e. to the matrix ele- 
ments {J^y-^'''^'' and (J^g) (dashed thin line): At Q"^ = 0, they add positively 
to the neutron charge but for increasing their contribution quickly drops. As 
a consequence, stronger axialvector correlations tend to suppress the peak in the 
form factor which the data suggest to be located at Q"^ = 0.2 . . . 0.5 GeV^. Again, 
the seagulls and the exchange quark graph provide for the correct neutron charge. 

The results for the proton magnetic form factor are shown in figure 3.7. Cor- 
respondingly, we have collected the single diagram contributions to the magnetic 
moment /i — G'm(O) in table 3.2. There, the influence of two parameters is vis- 
ible. First, in contrast to non-relativistic constituent models, the dependence of 
the proton magnetic moment on the ratio M„/mg is stronger than linear. As a 
result, the quark impulse contribution to /Xp with the scalar diquark being specta- 
tor, which is the dominant one, yields about the same for both sets, even though 
the corresponding nucleon amplitudes of Set I contribute about 25 % less to the 
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Figure 3.6: Contributions of the single diagrams to the nucleon electric form factors for 
Set I. They are obtained by applying eq. (3.8) to the current matrix elements. 




Figure 3.7: The proton magnetic form factor for Set I and II along with the dominat- 
ing contributions from the impulse approximation. Experimental data are provided by 
ref. [60]. 
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Set 1 


^4- TT 

Set 11 


experiment 






1.35 


1.33 






/ \ax— ax 


0.44 


0.08 






{J^jsc-s. + ( JM )ax-sc 


0.43 


0.24 








0.25 


0.22 






SUM 


2.48 


1.92 


2.79 




SUM 


-1.53 


-1.35 


-1.91 


jJ'p ~\~ fJ^n 


isoscalar 


0.95 


0.57 


0.88 




isovector 


4.01 


3.27 


4.70 



Table 3.2: Magnetic moments of proton and neutron. The contributions from the indi- 
cated current matrix elements are to be understood as coming with the appropriate flavor 
factors, cf. eq. (3.36). 

norm than those of Set II. Secondly, the scalar-axialvector transitions contribute 

equally strong (Set 1) or stronger (Set 11) than the spin flip of the axialvector di- 
quark itself. While for Set 11 (with weaker axialvector diquark correlations) the 
magnetic moments are about 30 % too small, the stronger diquark correlations of 
Set I yield an isovector contribution which is only 15% below and an isoscalar 
magnetic moment slightly above the phenomenological value. 

Stronger axialvector diquark correlations are favorable for larger values of the 
magnetic moments as expected. If the isoscalar magnetic moment is taken as an 
indication that those of Set 1 are somewhat too strong, however, a certain mismatch 
with the isovector contribution remains, also with axialvector diquarks included. 
Looking at the left panel of figure 3.7, we note that the axialvector diquark dia- 
gram (dashed thin line) falls very slowly with increasing Q"^ and eventually causes 
a violation of the observed dipole shape. For Set II, this effect is absent due to the 
weakness of the axialvector correlations. 

We can add another twist to the story by comparing our results to the recently 
measured ratio HpGe/Gm for the proton [69], cf. figure 3.8. The ratio obtained 
from Set 11 with weak axialvector correlations lies above the experimental data, 
and that for Set 1 below. The experimental observation that this ratio decreases sig- 
nificantly with increasing Q"^ (about 40 % from = to 3.5 GeV^), can be well 
reproduced with axialvector diquark correlations of a certain strength included. 
The reason for this is the following: The impulse-approximated photon-diquark 
couplings yield contributions that tend to fall off slower with increasing Q'^ than 
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t JLab Hall A Coll, PRL 84 (2000), 1398 

Set I 




0.5 1 1.5 2 2.5 3 3.5 
Q' [GeV'] 

Figure 3.8: The ratio {jip Ge)/Gm compared to the data from ref. [69]. 



those of the quark. This is the case for both, the electric and the magnetic form 
factor. If no axialvector diquark correlations inside the nucleon are maintained, 
however, the only diquark contribution to the electromagnetic current arises from 
( Jfp)^^~^^, see eqs. (3.36,3.37). Although this term does provide a substantial con- 
tribution to Ge, its respective contribution to Gm is of the order of 10~^. This 
reflects the fact that an on- shell scalar diquark would have no magnetic moment 
at all, and the small contribution to Gm may be interpreted as an off-shell ef- 
fect. Consequently, too large a ratio /j^pGe/Gm results, if only scalar diquarks 
are maintained [23]. For Set II (with weak axialvector correlations), this effect is 
still visible, although the scalar-to-axialvector transitions already bend the ratio to- 
wards lower values at larger Q^. These transitions almost exclusively contribute to 
Gm, and it thus follows that the stronger axialvector correlations of Set I enhance 
this effect. Just as for the isoscalar magnetic moment, the axialvector diquark 
correlations of Set I tend to be somewhat too strong here again. 

3.2.3 Conclusions 

The numerical results for the electric form factors support the already emphasized 
importance of using the Ward-Takahashi identities for constructing a suitable cur- 
rent operator. The correct value for Ge{0) = gjv, the nucleon charge, has been 
obtained not only for the two parameter sets that have been presented here but also 
for a general choice of parameters and of the momentum partitioning r/. Further- 
more, the proton and neutron electric form factors are described very well in the 
diquark-quark-model. 
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The strength of axialvector correlations within the nucleon turns out to be 
constrained by two observations: 

• The impulse approximation graphs involving tend to suppress the neu- 
tron electric form factor. This becomes effective only for stronger correla- 
tions than presented here, though. 

• Although larger axialvector correlations are favorable for the magnetic mo- 
ments as expected, the ratio UpGE/Gm provides tight bounds. As a result, 

should contribute to the nucleon norm around 20 %. 

As has been mentioned, a considerable improvement in the description of the nu- 
cleon magnetic moments has been achieved compared to the scalar diquark sector. 
The anomalous transitions between scalar and axialvector diquark are in this con- 
nection equally important as the contributions of the axialvector diquark itself. 
We hereby confirm the conjecture from ref. [41] where the nucleon magnetic mo- 
ments have been ascertained in a simple additive diquark-quark picture including 
scalar and axialvector diquark. 

Nevertheless, the isovector magnetic moments are still too small by about 
15 %. We discuss two possibilities that might account for the mismatch. 

• As we will see in the next chapter on the pionic and the axial form factor 
of the nucleon, vector diquarks need to be included in the parametrization 
of the two-quark correlations on theoretical grounds. Although their con- 
tributions to the nucleon binding energy are expected to be negligible, cf. 
the discussion in the beginning of section 2.2, the photon coupling of the 
vector diquarks could be strong enough to compensate the underestimated 
magnetic moments. 

• The tree-level form of the employed quark and diquark propagators is cer- 
tainly too naive an assumption. As an example, a non-trivial quark prop- 
agator would lead to some effective anomalous magnetic moment of the 
(constituent) quark via the associated non-trivial quark-photon vertex. The 
effect of the latter would be visible, e.g., in an enhanced magnetic moment 
contribution of the quark impulse approximation diagram. In chapter 5 we 
will therefore discuss some aspects of a dressed quark propagator and em- 
ploy an especially simple parametrization (which avoids the free-particle 
threshold) in calculations of the octet-decuplet spectrum and of form fac- 
tors. 



Chapter 4 

Strong and Axial Form Factors 



In this chapter we will investigate the pion-nucleon form factor as an example of 
a strong form factor and the isovector axial form factor which is of interest in the 
description of nucleon weak interactions. In the first case, we are interested in the 
matrix elements of a pseudoscalar current operator Jg and in the second case it 
is the pseudovector' current operator J^'^ whose matrix elements are the object 
under scrutiny. The superscript a labels the component of the isovector triplet 
which both current operators belong to. 

If the squared momentum transfer is equal to the pion mass, = ~"^^> 
then the pseudoscalar matrix element is entirely determined by the pion-nucleon 
coupling constant q-^nn which is accessible experimentally. There appears to be 
no experimentally feasible scattering process which is dominated by one-pion ex- 
change such that the form factor could be determined away from the pion mass 
shell. Therefore corresponding theoretical results from microscopic nucleon mod- 
els such as the diquark-quark model only serve as a guideline for effective nucleon 
potential models which incorporate one-pion exchange. 

The matrix element of the pseudovector current at = is determined by 
the (5 decay and the corresponding weak decay or axial coupling constant qa is 
known quite accurately^. For spacelike Q^, the form factor may be determined by 
neutrino-deuteron scattering [85] but these experiments are difficult to perform. 
Thus the axial structure of the nucleon is known far less precisely than its electro- 
magnetic structure. 

As it is well-known, these seemingly unrelated constants qa and g-nNN are 
linked with each other by the famous Goldberger-Treiman relation, see below. 
It can be derived from the assumption of chiral symmetry of an underlying La- 

' We reserve the label pseudovector for the external currents to distinguish them from axialvec- 
tor diquarks. Both objects have of course the same Lorentz and parity transformation properties. 

^Although in 1959 the constant was assumed to be known precise enough, (?a|i959 = 1-17 ± 
0.02, it is quite different from the presently known value, S'a|20oo = 1.267 ± 0.0035 [83, 84]. 
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grangian which conserves the pseudovector current. Additionally it is assumed 
that this chiral symmetry is spontaneously broken and pions appear as mass- 
less Goldstone bosons. Now the physical pion is massive, albeit its mass is 
small, so chiral symmetry is explicitly broken and the pseudovector current can 
be only partially conserved (PCAC hypothesis). However, it can be shown that the 
Goldberger-Treiman relation is still valid under the assumption that the axial form 
factor is slowly varying in the vicinity of = 0. Furthermore, the Goldberger- 
Treiman relation holds experimentally on the five per cent accuracy level which, 
among other observations, lends support to the concept of spontaneously broken 
chiral symmetry and the partially conserved pseudovector current. 

We tried to show by these introductory arguments that the pion-nucleon cou- 
pling constant and the weak decay constant are related with each other in a model 
independent way. Thereby we imply that if one were to start with a chirally 
symmetric quark theory one should always recover the Goldberger-Treiman re- 
lation for the nucleon, similar to the case of electromagnetic interactions where 
the nucleon charges should always come out correctly. Again similar to the elec- 
tromagnetic case, there are powerful chiral Ward identities which constrain the 
pseudovector current vertices. 

Unfortunately we will not be able to comply with the demand of a conserved 
pseudovector current and to actually derive the Goldberger-Treiman relation in 
the diquark-quark model in its present form. As will become clear in section 
4.1, vector diquark correlations are inevitably needed to conserve the pseudovec- 
tor current. These are neglected in our model, for their influence on the binding 
energy is assumed to be small and they introduce six additional scalar functions 
in the nucleon vertex function therefore being a substantial numerical com- 
plication. Our calculations presented in this chapter aim therefore at quantitative 
results for the expected discrepancy and at possible explanations for the measured 
weak coupling constant. In contrast to the anomalous magnetic moments, the non- 
relativistic quark model estimate = 5/3 is 31% larger than the measured value, 
qa = 1-27, and therefore a relativistic treatment such as taken up here might offer 
an explanation for this difference. 

We will begin with fixing our conventions for the matrix elements and re- 
derive quickly the Goldberger-Treiman relation. Section 4.1 explains the neces- 
sity of vector diquarks for pseudovector current conservation. In section 4.2 we 
are concerned with the construction of the pseudoscalar and pseudovector current 
operator for our model. Here we start with a suitable quark-pion vertex which in 
conjunction with a chiral Ward identity yields the quark-pseudovector vertex. The 
structures and strengths of the couplings of the diquarks to the pion and the pseu- 
dovector current are obtained from resolving the diquarks in a way similar to their 
electromagnetic couplings. The numerical results, again obtained by employing 
the two parameter sets from table 2.8, are presented in section 4.3. We will find 
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a slight overestimation of qa and a violation of the Goldberger-Treiman relation 
of around 15 %. The axial form factor results for finite Q'^ are consistent with the 
presently available experimental results. 

Let us start with a few definitions. The (spin-summed) matrix element of the 
pseudoscalar current J5 is parametrized as 

{Pf\ Jt |P,) = A+(P^) rV^/.TVivlQ') A+(P.) , (4.1) 

which defines the pion-nucleon form factor g-^NNiQ"^)- As before, Q = Pf — 
holds. For numerical evaluation, we will choose the neutral component, J|, and 
sandwich between proton flavor states, indicated by the subscript p. Straightfor- 
ward algebra allows to extract the form factor as the following Dirac trace, 

g.NNiQ') = Tr [p(P^| J| |P,)p 7^] . (4.2) 

We take for granted in the following that g-n-NNi^) is close to the pion-nucleon cou- 
pling constant g^^NN = gnNNi^'iTt-l) as consistent with the assumption of PCAC. 

The matrix elements of the pseudovector current are parametrized by the form 
factor gA{Q^) and the induced pseudoscalar form factor gp{Q^), 

{Pf\ \Pi) = A+(P/) y [iri'gAm + Q'l'gpm] ^"-(Pi) ■ (4.3) 

For — > 0, the induced pseudoscalar form factor will be dominated by the pion 
pole, 

.i. ,p(g^, = 2^*^ . (4.4, 

We denote the pion decay constant by = 92.7 MeV. The Goldberger-Treiman 
relation follows from current conservation, Q^^{J^'^) — 0, and after replacing gp 
by eq. (4.4) we find, 

gA(0) = Ug.NNiO)/Mn . (4.5) 

The regular part of the matrix element, gA{Q^), and the pion pole part, gp{Q^), 
can be extracted from eq. (4.3) as follows: 



^ I 4M2 



(4.6) 



gpm - ^ {gAm - ^ Tr [p(P,| 4- |P.)p Q^^'] ) . (4.7) 

Having specified our conventions, we turn now to a short discussion about the 
consequences of excluding vector diquarks from our considerations. 
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4.1 Scalar and vector/axialvector diquarks as chiral 
multiplets 

The necessity of including vector diquarks along with the axialvectors was ob- 
served only very recently in ref. [86] within the context of an NJL model. Re- 
member that diquarks as separable 2-quark correlators are in the NJL model solu- 
tions of the t matrix integral equation, with the interaction kernel given in lowest 
order in the coupling constant. The separable interaction kernel describes here the 
4-quark point interaction: two quarks in a channel with certain diquark quantum 
numbers interact locally with two antiquarks with respective antidiquark quantum 
numbers, see e.g. eq. (2.41) for the scalar channel. The statement of ref. [86] 
is summarized as follows: Chiral invariance of the NJL Lagrangian requires the 
inclusion of quark-quark interaction kernels with both vector and axialvector di- 
quark quantum numbers and the same coupling constant. Solutions for the separa- 
ble vector diquark t matrix do not have the same form as the axialvector t matrix, 
e.g., a possible vector diquark mass pole will appear at larger invariant momentum 
values than an axialvector pole. Nevertheless vector diquarks have to be carried 
all the way through the nucleon Faddeev problem in order to comply with chiral 
invariance. 

In our model, we have parametrized the diquark correlations without recur- 
ring to a specific model Lagrangian. To adapt the above argument to our case, 
we introduce interpolating diquark fields as done in section 3.1. Let dl denote 
scalar diquark creation operators and dt^\ d^J^^ the ones for vector and axialvector 
diquarks. Scalar and vector diquarks are flavor singlets, the axialvector diquark 
belongs to a flavor triplet, with its components denoted by the index h. The appro- 
priate color indices are suppressed as they are unimportant for the argument. The 
interaction of (pointlike) diquarks with quarks can now be written as an interaction 
Lagrangian, 

Qa {q^ Ci^'^ t'^t^ q ) d^^^ + herm. conj. , (4.8) 

where q denotes as usual the quark field. A combined infinitesimal chiral and 
isospin variation of the quark field reads 

Sq ^ -iUhT^-fU , (4.9) 

with LOb denoting a set of three infinitesimal parameters. Straightforward evalua- 
tion to first order in the variation yields 

5{q'^C^^T^q) = 0, (4.10) 

5{q^CiYl^r^q) = -iub {q^ Ci^ r^r\) , (4.11) 

5 {q^CiYr'^r^q ) = -iUb {q^CiYl^r'^q) . (4.12) 
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We demand that the interaction Lagrangian (4.8) be invariant under an infinitesi- 
mal chiral variation. Substituting eqs. (4.10^.12) into eq. (4.8), we see that this 
leads to the requirements 

5d'f ^iujbd':\ Sd'i^ = iui,dl^^ and ga = 9v. (4.13) 

Vector and axialvector diquarks mix under a chiral (isovector) transformation. The 
scalar diquark channel is unaffected by the transformation as its variation is zero. 
This is equivalent to the statement that a pseudovector (isovector) current does not 
couple to the scalar diquark. Summarizing these findings we have found that the 
scalar diquark on the one side and vector/axial vector diquarks on the other side 
belong to different chiral multiplets [86]. 

We learn something else. The above argument is only valid for a local quark- 
diquark interaction. If we introduced extended diquarks by a non-local interaction 
with some momentum dependent vertex we would find additional chiral seagull 
vertices as a result of the variation. This is very similar to the electromagnetic 
case. Although we derived the seagulls in section 3.1.2 by requiring gauge invari- 
ance via a Ward-Takahashi identity for the Bethe-Salpeter interaction kernel, they 
can also be obtained from an interaction Lagrangian with momentum-dependent 
vertices [76] in much the same way as indicated here. 

In the following we will neglect both the chiral seagulls and the vector di- 
quark correlations. As we remarked earlier, the expected high mass of the vec- 
tor diquarks would suppress the strength of the corresponding vector correlations 
within the nucleon. Another argument is concerned with the numerical feasibil- 
ity. In the Dirac algebra, the vector correlations could be decomposed similarly to 
the axialvector correlations, leading to six additional scalar functions. This would 
increase the numerical effort considerably and since the numerical calculations 
are already very involved at this stage, the inclusion of vector diquarks is beyond 
the scope of this thesis. On the other hand, we have no serious argument for the 
neglect of the seagulls. However, we can estimate from the numerical results pre- 
sented in section 4.3 that their contribution to the weak coupling constant should 
be below ten per cent. 



4.2 Pseudoscalar and pseudovector current opera- 
tor 



First we will specify the quark vertices with pseudoscalar and pseudovector cur- 
rent operators and then continue with a discussion of the effective diquark vertices 
that can be obtained by coupling the currents to an intermediate quark loop. 
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As we use a free fermion propagator for the quark, eq. (2.64), the scalar func- 
tions A, B in the parametrization = —i'pA — B reduce Xo A—1 and B — niq. 
This suggests that we use for the pion-quark vertex^ 

r,,.-,=£M^)=_,a^, (4.14, 

and discard the three additionally possible Dirac structures. The momenta of out- 
going and incoming quark are k and p, respectively. The reason is that in the 
chiral limit eq. (4.14) represents the exact pion Bethe-Salpeter amplitude (on- 
shell, i.e. = 0) for equal quark and antiquark momenta, since in this limit, 
with constant A, the Dyson-Schwinger equation for the scalar function B agrees 
with the Bethe-Salpeter equation for a pion of zero momentum. Of course, the 
subdominant amplitudes should in principle be included for physical pions (with 
momentum = —ml), when solving the Dyson-Schwinger equation for A and 
B and the Bethe-Salpeter equation for the pion in mutually consistent truncations 
[48, 87, 64]. The generalization to pion off-shell momenta as used in eq. (4.14) is 
by no means unique but certainly the most simple one. 

We use chiral symmetry constraints to construct the pseudovector-quark vertex 
r^^. In the chiral limit, the Ward-Takahashi identity for this vertex reads, 

Q^Tl^ = S-\k)^' + ^'S-\p) , {Q = k-p). (4.15) 

To satisfy this constraint we use the form of the vertex proposed in ref. [88], 

r^,, = -nV + ^2Ar5,, . (4.16) 

The physical picture behind this form of the vertex is clear. In the spontaneously 
broken phase of chiral symmetry the quark acquires mass and at the same time 
the Goldstone boson pole must appear in the full quark-pseudovector vertex, and 
the Goldstone boson pole in turn is proportional to the generated quark mass. The 
second term which contains the massless pion pole does not contribute to qa as 
can be seen from eq. (4.6). 

The pion and the pseudovector current can couple to the diquarks by an in- 
termediate quark loop. As for the anomalous contributions to the electromagnetic 
current, we derive the Lorentz structure of the diquark vertices and calculate their 
effective strengths from this quark substructure of the diquarks in appendix B.1.2. 
For momentum definitions at these vertices we refer to figure 3.1. 

•'We drop the flavor indices on the vertices in the following. Each pseudoscalar vertex comes 
with a factor of r" and each pseudovector vertex with a factor of as these are widely used 
conventions. 
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As already mentioned in the last section, no such couplings arise for the scalar 
diquark. This can be inferred from parity or covariance or, simpler, one considers 
the flavor trace of an isovector current between two isoscalar states (scalar diquark 
vertex) which yields zero. The axialvector diquark and the pion couple by an 
anomalous vertex. Its Lorentz structure is similar to that for the photon-induced 
scalar-to-axialvector transition in eq. (3.32), 

rfa. = ^^^"'''^(^^'^ + ^'^)'^^- (4.17) 

Here, a and /3 are the Lorentz indices of outgoing and incoming diquark, respec- 
tively. The factor rriq/ comes from the quark-pion vertex (4.14) in the quark 
loop (see appendix B.1.2), and the nucleon mass was introduced to isolate a di- 
mensionless constant k^^. 

The pseudovector current and the axialvector diquark are also coupled by 
anomalous terms. As before, we denote with a and (3 the Lorentz indices of 
outgoing and incoming diquark, respectively, and with ji the pseudovector index. 
Out of three possible Lorentz structures for the regular part of the vertex, 

Pd e'^^^^PdQ^ , e^'^'^PQP and e'^°^^(pd + 

only the last term contributes to qa in the limit Q ^ 0. We furthermore verified 
numerically that the first two terms yield negligible contributions to the form fac- 
tor also for nonzero Q. Again, the pion pole contributes proportionally to Q'^, and 
our ansatz for the vertex thus reads 

For both coupling strengths in the vertices Fg^^. and Fg*^^ we roughly obtain ~ 
'^liax — '^■^ slightly dependent on the parameter set, see table B.2 in appendix 
B.i.2. 

Scalar-to-axialvector transitions are also possible by the pion and the pseu- 
dovector current. An effective vertex for the pion-mediated transition has one free 
Lorentz index to be contracted with the axialvector diquark. Therefore, two types 
of structures exist, one with the pion momentum Q, and the other with any combi- 
nation of the diquark momenta and A;^. If we considered this transition as being 
described by an interaction Lagrangian of scalar, axialvector and pseudoscalar 
fields, terms of the latter structure would be proportional to the divergence of the 
axialvector field which is a constraint that can be set to zero. We therefore adopt 
the following form for the transition vertex. 



Ksa=-^<a^Q'- (4.19) 
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Here the Lorentz index of the participating axialvector diquark is given by (3. This 
vertex corresponds to a derivative coupUng of the pion to scalar and axialvector 
diquark. 

The pseudovector-induced transition vertex has two Lorentz indices, denoted 
by n for the pseudovector current and P for the axialvector diquark. From the mo- 
mentum transfer Q'* and one of the diquark momenta altogether five independent 
tensors can be constructed, 

5^^, Q^Q^ , Q^p^,, p^.QP and • 

Note that a term constructed with the totally antisymmetric tensor would have 
the wrong parity. We assume, as before, that all terms proportional to Q^^ are 
contained in the pion part. They do not contribute to qa anyway. From the diquark 
loop calculation in appendix B. 1 .2 we find that the terms proportional to p^^Q'^ and 
p^p'l can again be neglected with an error on the level of one per cent. Therefore, 
we use a vertex of the form, 

^ta - iM^<sa + ^'^U^isa- (4-20) 

For the strengths of these two transition vertices we obtain k^^ ~ 3.9 and (on 
average) ^ 2.1, cf. table B.2 in appendix B.1.2. 

The vertices Fg^^, Fg^^^ for the reverse transitions are obtained by simply re- 
versing the sign of Q. 

Having constructed all vertices, we can now perform the flavor algebra and 
write down the pseudoscalar current matrix element in the extended impulse ap- 
proximation. As done in section 3. 1 . 1 we introduce quark current matrix elements 
by 

( Js,,)''"'' = D-^ Fs,, , (4.21) 

where this time no flavor matrix factor comes with F5 ^ as we compute the flavor 
factors separately. Likewise scalar/axialvector diquark and the transition matrix 
elements are defined. With these definitions we find for the matrix element of the 
extended impulse approximation 

^ {{\saY" + . (4.22) 

As we neglect seagull contributions, the only contribution from the Bethe-Salpeter 
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kernel is due to the coupling to the exchange quark and it reads explicitly 

p(P/| j||P,)f = (J5,e.) (4.23) 



ex,ss 




ex,sa 




ex, as 


= 1^3 x'iPi) S^iq) (Fs,,)^ S^q') (P2) , 


ex,aa 


= ^x^(pi)^^(g)(r5,,r^V)r(P2). 



The diagrams of the extended impulse approximation are equivalent to those dis- 
played in figure 3.2 and likewise the exchange quark diagram is similar to the 
first diagram in figure 3.4. The black dot which represents the quark-photon or 
diquark-photon vertex should be replaced with the appropriate quark-pseudoscalar 
or diquark-pseudoscalar vertex. The matrix elements for the pseudovector current 
differ from eqs. (4.22,4.23) only by an overall factor of 1/2 after replacing the 
pseudoscalar vertices by their pseudovector counterparts. 

4.3 Numerical results 

Having defined the pseudoscalar and pseudovector current operators in the last 
section and having fixed all unknown coupling constants by appropriate quark 
loop calculations, cf. appendix B. 1 .2, we can now compute the form factors for the 
two parameter sets of table 2.8. Regarding the numerical procedure, the method 
and its intricacies have been lined out in section 3.2.1 and all the remarks made in 
this section also apply here. We summarize: impulse approximation diagrams are 
computed in the Breit frame employing boosted vertex functions $ and taking into 
account the residue contributions, cf. appendix B.2. The exchange quark diagram 
can only be computed using boosted wave functions and therefore the related 
numerical accuracy becomes very limited beyond momentum transfers Q'^ > 2.5 
GeV^. 

Examining g^^jv at (0) which is assumed to be close to the physical pion-nucleon 
coupling and the weak coupling constant gAi^), we find large contributions to 
both arising from the scalar-axialvector transitions, cf. table 4.1. As mentioned 
in the previous section, the various diquark contributions violate the Goldberger- 
Treiman relation. Some compensations occur between the small contributions 
from the axialvector diquark coupling in impulse approximation and the compar- 
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Set I Set II 

5'7rivjv(0) gA{0) g7TNN{0) 5'a(0) 



ax— ax 
ax— ax 



/ tIWVsc— sc 
I tM \sc-ax _|_ / tM \£ 

SUM 



7.96 


0.76 


9.25 


0.86 


0.50 


0.04 


0.10 


0.01 


1.44 


0.18 


0.34 


0.04 


5.66 


0.39 


3.79 


0.22 


1.69 


0.12 


2.70 


0.22 


17.25 


1.49 


16.18 


1.35 





Q-kNN 




9a 


experiment 


13.14± 0.07 


[89] 


1.267 ± 0.0035 [84] 




13.38±0.12 


[90] 





Table 4.1: Various contributions to g-j^NNi^) and (7^(0), labelled as in table 3.2. 

atively large ones from scalar-axialvector transitions which provide the dominant 
effect to yield 5fA(0) > 1. 

Summing all these contributions, the Goldberger-Treiman discrepancy, 

^ ^ M[Ar(0) _ 

amounts to 0. 14 for Set I and 0. 18 for Set II. The larger discrepancy for Set II (with 

weaker axialvector correlations) is due to the larger violation of the Goldberger- 
Treiman relation from the exchange quark contribution in this case. This contribu- 
tion is dominated by the scalar amplitudes, and its Goldberger-Treiman violation 
should therefore be compensated by appropriate chiral seagulls as already men- 
tioned in section 4.1. For Set II which is dominated by scalar correlations within 
the nucleon, their contribution to gAiS^) can be estimated to be around 0. 1 because 
the seagulls together with the contributions from {J^}q)^'^~^'^ and {J^).x) should 
approximately obey the Goldberger-Treiman relation. 

The strong and weak radii are presented in table 4.2 and the corresponding 
form factors in figure 4.1. Experimentally the axial form factor is known much 
less precisely than the electromagnetic form factors. In the right panel of figure 
4.1 the experimental situation is summarized by a band of dipole parametrizations 
of gA that are consistent with a wide-(energy)band neutrino experiment [85]. Be- 
sides the slightly too large values obtained for which are likely to be 
due to the PCAC violations of axialvector diquarks and of the missing seagulls as 
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Nucleon Axial Form Factor 
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Figure 4.1: The pion-nucleon form factor g-KNN{Q'^) and the axial form factor (^^(Q^). 
'Diquark diagrams' labels the sum of the impulse-approximate axial vector contributions 
and scalar-axialvector transitions. The shaded region in the right panel represents the 
uncertainty in as determined from quasi-elastic neutrino scattering when a dipole form 
is fitted to both, the vector and the axial form factor [85]. 



discussed in section 4. 1, our results yield quite compelling agreement with the ex- 
perimental bounds. As in the case of the electromagnetic form factors we observe 
a dominance of the diquark contributions for larger Q^. 

Turning to the pion-nucleon form factor we have fitted our results to an n-pole, 

The form factor of Set I is close to a monopole, n = 1, with cut-off Ai = 0.6 GeV, 
whereas the one of Set II is better approximated by a dipole, n = 2, with cut-off 





Set I 


Set II 


experiment 


r^NN [fm] 


0.83 


0.81 




rA [fm] 


0.82 


0.81 


0.70±0.09 



Table 4.2: Strong radius r.n-NN and weak radius va, the experimental value of the latter 
is taken from [85]. 



100 



Strong and Axial Form Factors 



A2 = 0.9 GeV. The exponent n is determined by the strength of the diquark contri- 
butions whereas the cut-off has to comply with the slope at = (the radius). It 
is not a big surprise that the proton electromagnetic radii and the pionic/axial radii 
are very close to each other since the dominating impulse approximation diagrams 
have a very similar behavior for all observables. 

As we have mentioned in the beginning of the chapter, pion-nucleon form fac- 
tors are used in one-boson exchange potential ansdtze to describe the A^A^ force. 
Most of the potentials that fit the data employ a monopole with a large cut-off, 
Ai > 1.3 GeV, see e.g. ref. [91]. However, other nucleon models, lattice and 
QCD sum rule calculations indicate a monopole behavior with a much smaller 
cut-off Ai ^ 0.5 — 0.95 GeV, cf. ref. [92] and references therein. Our calculations 
support this point of view. 

Summarizing the results of this section, we have found moderate violations 
of PCAC and the Goldberger-Treiman relation. These violations are estimated to 
be partitioned equally between the effects of neglecting chiral seagulls and vector 
diquarks. The inclusion of axialvector diquarks is essential to obtain values for the 
weak coupling constant qa > ^- Oi particular importance for this result are the 
transitions between scalar and axialvector diquarks, additionally they restrict the 
axialvector correlations within the nucleon to be small. For the pion-nucleon form 
factor, our model calculations support a monopole- or dipole-like fall-off with a 
cut-off well below 1 GeV. 



Chapter 5 

Effective Confinement 



In this chapter we will investigate the possibility of incorporating confinement into 
the diquark-quark model by a suitable modification of quark and diquark propa- 
gators. As will become clear in the following, this will enable us to calculate 
the mass spectrum of octet and decuplet baryons. Electromagnetic form factors 
recomputed with the modified propagators reveal a drastic increase on the mag- 
netic moments, being almost too much of an improvement. We have done the 
calculations for pointlike and for extended diquarks to demonstrate the necessity 
of the latter for a consistent description of masses and electromagnetic structure. 
Although the results for the spectrum and the magnetic moments are encouraging, 
there are serious difficulties associated with this approach of implementing con- 
finement, both motivated from "theoretical" and "phenomenological" arguments. 

Remember that the propagators along with the diquark vertices are the only 
ingredients of the model after the full 3 -quark problem has been reduced, cf. the 
summary at the end of section 2.1. By choosing the most simple ansdtze for them, 
i.e. free spin- 1/2 and spin-O/spin-1 propagators for quark and diquark we have 
succeeded to describe quite successfully various spacelike nucleon form factors. 
However, shortcomings remain which might be cured by allowing for non-trivial 
propagators. These shortcomings are: 

• To incorporate the A into the description, rather large quark and axialvector 
diquark masses have to be chosen, as their sum must be larger than the 
mass of the A. This problem extends to the strange quark mass and strange 
diquarks when aiming at including the whole octet-decuplet spectrum. 

• The isovector part of the nucleon magnetic moments is too small, even after 
inclusion of the axialvector diquarks. The difference is substantial for pa- 
rameter sets which fit the A mass. A reason for this, besides too large a con- 
stituent quark mass, might be the insufficiency of the bare quark-photon ver- 
tex to account for the quark contributions to the magnetic moments. Gauge 
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invariance demands, however, that a nonperturbative quark-photon vertex is 
intimately connected with a nonperturbative quark propagator. 

• Hadronic reactions where timelike momenta of the order of 1 GeV are de- 
posited onto nucleons, e.g. meson production processes, cannot be de- 
scribed within the model in its present form. The free-particle poles of 
quark and diquark cause unphysical thresholds in these processes. 

As we see, the lack of confinement which is manifest in the poles of the quark 
and diquark propagators used so far is the main motivation to consider alterations 
in the model premises. At this stage we hope that the necessary modifications in 
the quark-photon vertex which are required by gauge invariance will improve the 
results on the magnetic moments as well. 

5.1 Confinement in model propagators 

Confinement is understood as the absence of colored states in the spectrum of 
observed particles. From a phenomenological point of view (as adopted here), 
calculations of S matrix elements using diagrams with internal quark loops must 
not have any imaginary parts which are associated with the singularities of the 
quark propagators. So, either they are absent or their contributions cancel in some 
manner [21]. In this way, one is led to consider the following possibilities, which 
are not necessarily the only ones: 

a. Propagators are entire functions in the whole complex plane [93]. If they 
are to be analytic, they must possess an essential singularity in the infinite. 
An example is the function f{x) — e^. 

h. Propagators have complex conjugate poles [94]. In some sense the poles 
correspond to virtual excitations that cancel each other in physical ampli- 
tudes. Here an example is provided by the function f{x) — x/ (a^ + x^). 

We will investigate the consequences of propagators being entire analytic func- 
tions. We insist here on analyticity in order not to loose the relativistic reparame- 
trization invariance related to the choice of the momentum partitioning parame- 
ter Tj.^ A word of caution is in order here. The Wick rotation which connects 

'As explained in section 2.3. 1, i] invariance requires the following: If 5'(p, P; 771) is a solution 
of the Bethe-Salpeter equation then ^'(p + (772 — r]i)P, P; 772) with 772 ^ rji is also one. In the rest 
frame of the bound state, P = (0,zM„), this can be shown to be valid if the integration of the quark 
exchange kernel K^^{p, k, P) over the component k'^ can be shifted as k"^ ^ k'^ + (772 — 77i)iM„. 
This is of course possible, if the propagators are analytic in the complex domain of fc** given by 
(Refc4,ImA;^) e ((-00, +00), [0, (772 - ??i)Afn]). 
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Minkowski and Euclidean metric is no longer applicable as the propagators pos- 
sess the essential singularity at infinity. Furthermore we are aware that this prop- 
agator contradicts general arguments about the behavior of propagators in gauge 
theories with confinement. Following ref. [95], it should vanish for ^ oo in all 
directions of the complex -plane faster than the free propagator (p is the particle 
momentum). Furthermore the propagators should have cuts along the real pr- 
axis. Unphysical colored states associated with these cuts are then to be removed 
from the observable spectrum by other conditions such as a suitable definition 
of a physical subspace (cf. the long-known Gupta-Bleuler mechanism to remove 
scalar and longitudinal photons from the physical particle spectrum in QED [96]). 
It seems to be impossible to implement propagators which behave as indicated 
above and such a projection onto physical states in actual matrix element calcula- 
tions at present, though. Therefore we believe it is instructive to investigate some 
phenomenological advantages and deficiencies of an "entire" propagator. In prac- 
tical calculations, the quark propagator will be needed only on a limited domain in 
the complex plane and the essential singularity thus poses no practical difficulty, 
at least for the form factor calculations presented here. 

The second possible form for the propagator mentioned above suffers from 
only one principal deficiency, as there are complex conjugate pole associated with 
it. The asymptotic conditions in the complex plane are in this case easy to incor- 
porate. Actual calculations with this form are harder to perform and are currently 
under investigation. 

Before we turn to the propagator parametrization chosen here, we discuss 
some aspects of the quark propagator that are known so far. We regard the quark 
propagator to be given in its most general form, — —i'^A—B, with A, B being 
scalar functions. Information about the quark propagator may be obtained from 
calculations within the framework of Dyson-Schwinger equations, for a recent 
review see ref. [21]. In the near future, lattice results will provide us with more 
"empirical" information, for a recent study in quenched QCD see ref. [97]. Unfor- 
tunately, the quark propagator is usually computed in both schemes only for pos- 
itive real p^, where p is the Euclidean quark momentum, therefore the immediate 
benefit for phenomenological calculations is limited. However, these calculations 
have established an important generic feature. In rainbow truncation of the Dyson- 
Schwinger equations as well as in improved truncation schemes it could be shown 
that the quark acquires an effective constituent mass, i.e. M(0) = B{0)/A{0) 
is non-zero for a vanishing current quark mass. This feature of dynamical chiral 
symmetry breaking has also been observed in the already mentioned lattice cal- 
culation [97]. For phenomenological calculations, the numerical results for A, B 
have been fitted to entire functions which has been successfully used in the de- 
scription of (spacelike) meson [98] and nucleon properties [78, 79]. The most im- 
pressive success in this direction has been achieved in a combined rainbow/ladder 
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truncation scheme for quarks and mesons employing a specific gluon propagator 
which displays significant enhancement at momenta of the order of a few hun- 
dred MeV. Here the quark propagator has been computed from Dyson-Schwinger 
equations on all complex momentum points^ where it has been needed to solve 
the Bethe-Salpeter equation for mesons and to compute their electric form fac- 
tors [64, 65, 99]. Additionally the quark-photon vertex has been obtained from 
the solution of a corresponding inhomogeneous Bethe-Salpeter equation [100]. 
In these studies the solution for the quark propagator exhibited no poles for the 
(complex) momenta sampled by the meson calculations. These calculations sup- 
port the working hypothesis of a quark propagator having no poles, at least in the 
interesting momentum regime. 

In eq. (2.64) we have already indicated a possible parametrization of the quark 
propagator by multiplying the free one with a single dressing function C{p^ , m). 
We choose a function for C which is entire and analytic and removes the free- 
particle pole in the denominator of the free propagator, 



The essential singularity occurs here for — > — oo. The parameter d regulates 
the modification as compared to the free propagator. If d is large, the spacelike 
properties remain nearly unchanged, whereas for timelike momenta the propaga- 
tor blows up quickly. For small d <1, even the spacelike behavior is substantially 
modified. 

For simplicity, we apply the same modification to the propagators of scalar 
and axialvector diquarks, cf. eqs. (2.45,2.46). Additionally we choose { = 1 in 
eq. (2.46). This parameter is a gauge parameter when describing the interaction of 
a vector field with photons [75] . The photon vertex with the axialvector diquark 
has to be modified correspondingly. 

We note that due to the absence of poles the "masses" m^, m^c, ruax have no 
physical interpretation. They are merely width parameters of the propagator. 



Equipped with propagators which effectively mimic confinement we can calcu- 
late the mass spectrum of octet and decuplet baryons. This investigation has been 
reported in ref. [24]. In doing so, we have to extend all previous considerations 
to flavor SU{?)) and introduce explicit symmetry breaking as visible in the exper- 
imental spectrum. However, we will confine ourselves to the limit of unbroken 

^This was rendered possible through the analytical ansatz for the gluon propagator. 




(5.1) 
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isospin and, accordingly, the strange quark mass is the only source of symmetry 
breaking, ^ m„ = m^, where m^, nid denote the mass parameters of up and 
down quark. 

In order to limit the number of parameters we assume the scalar and axialvec- 
tor diquark mass parameters to be equal and introduce a diquark mass coefficient 
Cby 

^"c = T^ax = C{ma + nik) (a, h G {u, d, s}) . (5.2) 

Thus ab denotes the flavor content of the diquark. It will be interesting to see 
whether the mass differences between spin- 1/2 and spin-3/2 baryons can be ex- 
plained without an explicitly larger axialvector diquark mass parameter. 

Diquark 3x3 flavor matrices {t[ab))jf. are chosen antisymmetric for scalar di- 
quarks. 

Here we use for the A's the standard Gell-Mann matrices [40]. Axialvector di- 
quark matrices (t[ab]}jj^ are symmetric in their indices j, k, 

^M = ^A^ ^M = ^A^ ^[d.] = ^A^ ^^^^^ 

t[uu] = Sjl^kl ) t[dd] = ^j2^fc2 ) t[ss] = SjS^kS ■ 

To derive the Bethe-Salpeter equations for the octet (N, A, E, S) and decuplet 
baryons (A, S*, S*, 0,), one starts with their flavor wave functions prescribed by 
the Eightfold Way. The octet baryons have scalar diquark correlations, therefore 
this part of the flavor wave function is mixed antisymmetric. The axialvector 
diquark correlations for octet members are described in flavor space by a sum of 
two mixed antisymmetric states. For the decuplet members with only axialvector 
diquarks contributing these are either totally symmetric (A++, A^, Q) or mixed 
symmetric states. These wave functions may contain different mass eigenstates 
of quark-diquark pairs. For these mass eigenstates the system of coupled Bethe- 
Salpeter equations has to be solved, where the decomposition in the Dirac algebra 
according to section 2.4 can be used. As an example we consider the hyperon 
with its flavor wave function 

u{su) 

S /flavor = /T r 1 /i" r 1 ■ i^-^) 



The scalar diquark part, u{su) is clearly mixed antisymmetric and the axialvector 
diquark part is the normalized sum of {us — su)u and uus — suu, both also mixed 
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experi- 
ment 


pointlike 


extended 
diquark 








V{p) = 1 


A = 0.5 GeV 


rriu 
rUs 

Qa/gs 


[GeV] 
[GeV] 




0.50 
0.63 
0.73 
1.37 


0.56 
0.68 
0.60 
0.58 


Ma ^ 

Ms 

Ms 

Ms* ^ 
Ms. 
Mn ^ 


[GeV] 


1.116 
1.193 
1.315 
1.384 
1.530 
1.672 


1.133 
1.140 
1.319 
1.380 
1.516 
1.665 


1.098 
1.129 
1.279 
1.396 
1.572 
1.766 



Table 5.1: Octet and decuplet masses. 

antisymmetric. The Dirac-flavor wave function contains three independent tenns 
'^^ , X , 1 and ^!^, 1 with a total of 2+6+6= 1 4 unknown scalar functions which 

u{us)^ s[uu\ u[us\ 

shows the necessary increase of computer resources. 

We have listed the Bethe-Salpeter equations for all octet and decuplet baryons 
in appendix A.4 and quote from the results given there just one remarkable fact. 
Due to the symmetry breaking, there occurs mixing of the flavor singlet wave 
function with the scalar diquark correlations within the A. This mixing is absent 
in quark potential models. 

As the diquarks here can never be on-shell, we are not allowed to use the 
normalization conditions (2.59,2.60) to fix the quark-diquark coupling strengths 
Qs and Qa- Therefore we adjust them to the masses of nucleon and A, M„ = 0.939 
GeV and Ma = 1.232 GeV. 

In table 5.1 we show results for the masses of the remaining baryons {i) for 
pointlike diquarks and {ii) for a quadrupole-like diquark vertex function. For both 
cases the octet-decuplet mass difference is a result of solely the relativistic dynam- 
ics; it is even overestimated for extended diquarks. The parameter set using the 
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extended diquarks allows a considerable reduction of the ratio Qa/ Qs as compared 
to the ratio for pointlike diquarks. Still, axialvector correlations are somewhat too 
strong regarding the constraints which have been obtained in chapters 3 and 4. 

As can be seen from table 5.1, the calculated baryon masses are higher for the 
spin- 3/2 states than the sum of the quark and diquark mass parameters. Due to the 
absence of poles, it is possible to smoothly cross the tree-level thresholds given by 
the boundaries for the momentum partitioning parameter rj in eqs. (2.154,2.155), 
cf. the discussion in refs. [58, 22]. 

These results have been obtained with the choice = 1 in the propagator mod- 
ification function C^'^p, cf. eq. (5.1), corresponding to a quite substantial modifi- 
cation of the spacelike propagator properties. Since mainly these are important for 
solving the Bethe-Salpeter equation as well as for the calculation of form factors, 
we expect visible changes for the latter as well. We will now turn to this issue. 



5.3 Electromagnetic form factors 

In order to use the machinery developed in chapter 3 for calculating the form 
factors, we have to re-specify the basic quark and diquark vertices. To satisfy 
gauge invariance, the quark vertex must fulfill the Ward-Takahashi identity 

{k-pr = q, iS-\k) - S-\p)) , (5.6) 

with the quark propagator being dressed by the function C^'^p of eq. (5.1). The 
scalar functions A,B are given in terms of C^'^p by 

It turns out that the longitudinal part of the quark vertex is unambiguously given 
by the Ball-Chiu ansatz [101], labelled by the superscript BC, 

rj = * (rf''^ + t-^t) , (5.8) 

™,BC ^ + (5.9) 



. {k + pY 



{A{k')-A{p'))t±l-i{B{k')-B{p')) 



k^ — p 

Its exclusivity results if the following four conditions are valid 
1. It satisfies the Ward-Takahashi identity (5.6). 



2. It does not possess any kinematical singularities. 
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3. It transforms under parity, charge conjugation and time reversal as the free 
vertex. 

4. It reduces to the free vertex for A = 1 and B — niq. 

The transverse part of the vertex, V^j,, remains undetermined. We will put it to 
zero in the following which is certainly an oversimplification regarding the infor- 
mation about it which has been obtained during the last years. It has been shown 
in ref. [102] that the requirement of multiplicative renormalizability greatly re- 
stricts the tensor structure of the transversal vertex. For large spacelike momenta 
Q, where no resonances are present, this leads to a phenomenologically unim- 
portant vertex modification. For the small and moderate momentum transfers at 
which we investigate the form factors, a transversal term which contains the ef- 
fect of the p meson pole might turn out to be of importance for the quark-photon 
vertex. The already mentioned study in ref. [100] has obtained the full quark- 
photon vertex as a solution of an inhomogeneous Bethe-Salpeter equation. The 
quark-antiquark scattering kernel in the vector channel must enter the equation for 
the quark-photon vertex as well, thereby the solution for the full vertex contains 
the vector meson pole as the pole also appears in the full quark-antiquark scatter- 
ing amplitude. Besides having found this vector meson pole contribution (being 
transversal) to the quark-photon vertex, the study has numerically confirmed that 
the longitudinal part of the vertex is given by the Ball-Chiu recipe. Within this 
model context, the remnants of the vector meson pole have been shown to con- 
tribute about 30 % to the full pion charge radius r-,^, with the remainder being 
attributable to the Ball-Chiu vertex. 

In order to create a link to the old and, at this time, fairly successful picture 
of vector meson dominance, it would be interesting to investigate the influence of 
the vector meson pole also for the nucleon form factors. This requires, however, a 
tremendous effort to use an effective gluon interaction like the one from ref. [100] 
for the quark-antiquark and quark-quark scattering kernels to describe mesons and 
diquarks, to obtain the quark-photon vertex self-consistently, and to calculate or 
reliably parametrize the quark propagator in the complex plane where it is needed 
for calculations of nucleon properties. This is beyond the scope of this work. 

Having introduced the Ball-Chiu vertex as an approximation to the full quark- 
photon vertex and having discussed possible shortcomings, let us turn to the 
diquark-photon vertices. We make use of the Ward-Takahashi identities for scalar 
and axialvector diquarks, 

(k-pT = Qscia.] {{D-'Y^'^^k) - {D-'y-^\p)) , (5.10) 
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to fix the longitudinal part of the respective photon vertices. Neglecting again 
transversal parts, we find for the vertices. 



The last term in eq. (5.12) is transversal and therefore no dressing can be inferred 
from the Ward-Takahashi identity. For this reason, we do not modify this term. 
Compared to the bare vertex, eq. (3.31), the transversal term appears with a mod- 
ified prefactor, (1 + k) instead of k where k denotes the anomalous magnetic 
moment. This is a consequence of choosing the gauge parameter { = 1 in the 
dressed propagator [75]. 

Anomalous transitions between scalar and axialvector diquarks are of course 
possible as well. Their vertices remain unchanged, cf. eqs. (3.32,3.33). 

We can determine the unknown anomalous magnetic moment n appearing 
in eq. (5.12) and the strength of the scalar-to-axialvector transitions Kga from 
eqs. (3.32,3.33) by much the same procedure as described in section 3.1 and ap- 
pendix B.1.1. That is, we calculate the resolved diquark vertices from figure 3.1, 
of course by using the dressed quark propagator and quark-photon vertex given 
above, and re-adjust the diquark normalization (coupling) constants Qa g^^^^ 
and gs gl^^^ such that the differential Ward identity is fulfilled, cf. appendix 
B.1.1. Using these rescaled quantities g^^^ and gl^^'^, we determine k and Kga- 

For pointlike diquarks as used by the first parameter set of table 5.1 the calcu- 
lation of the resolved vertices leads to divergent integrals and therefore prohibits 
their use. We have performed the calculation for the second parameter set from 




(5.11) 



+ 



(5.12) 



K 



'sa 



extended diquarks 0.73 



1.33 



0.69 1.76 



Table 5.2: Rescaled diquark normalizations and constants of photon-diquark couplings 
for the parameter set with extended quarks, cf. table 5.1. 
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Figure 5.1: Electric form factors of both nucleons for the parameter sets of table 5.1. 

table 5.1, z.e. the set which employs extended diquarks. From the results found in 
table 5.2 we see that the anomalous magnetic moment and the transition strength 
are reduced compared to the previous calculations with free propagators, cf. table 
B.l. 

The form factors are now calculated as the sum of the extended impulse ap- 
proximation diagrams, cf. figure 3.2, and of the quark exchange kernel diagrams 
depicted in figure 3.4. The relevant formulae for the impulse approximation are 
given by eqs. (3.36,3.37) and by eqs. (3.48-3.50) for the exchange kernel contri- 
butions. 

Turning to the results for the electric form factors, shown in figure 5.1, we see 
the favorable influence of the diquark-quark vertices with finite width. They lead 
to a basically correct description of the proton electric form factor. The neutron 
electric form factor is also much better described by the parameter sets which 
employs extended diquarks. Here we see the already mentioned quenching effect 
of the axialvector correlations, though. They are quite strong as compared to the 
parameter sets from table 2.8, with the scalar diquark contributions to the norm 
integral reduced to 51 % (the numbers are 66 % for Set I and 92 % for Set II). 

Regarding the overall shape of the calculated magnetic form factors, depicted 
in figure 5.2, we obtain the same conclusion as before. The parameter set with 
extended diquarks fairly well describes the empirical dipole shape whereas point- 
like diquarks assume the nucleons to be far too rigid. It is interesting to observe 
that for the extended diquarks the calculated magnetic moments, fip = 3.32 and 
Hn = —2.14, now exceed the experimental values by 19 and 12 per cent, respec- 
tively. We find that the magnetic moment contribution of every single diagram 
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Figure 5.2: Magnetic form factors of both nucleons for the parameter sets of table 5.1. 
Experimental data for the neutron is taken from refs. [103, 104, 105, 106]. 

weighed with its contribution to the total nucleon charge is enhanced when com- 
paring with the corresponding ratios for the Sets I and II which employ bare prop- 
agator and vertices. Clearly the dressing of propagators and vertices has a quite 
substantial effect on the magnetic moments and a more refined choice for them 
may provide for an exact description of both electric and magnetic form factors. 



5.4 Discussion 

We have employed a simple parametrization of quark and diquark propagators 
which renders them to be devoid of poles in the complex momentum plane ex- 
cept for an essential singularity at infinity. Thereby we could demonstrate that 
an overall satisfying description of both the electromagnetic form factors and the 
mass spectrum of octet and decuplet baryons has been achieved. However, these 
results have to be taken with a grain of salt as has been already alluded to. Gen- 
eral arguments about the functional form of the quark propagator in QCD simply 
prohibit such a behavior for timelike momenta as has been assumed here. Since 
mass spectrum and form factor calculations test the propagator only to a moder- 
ate extent in the timelike region (Re(p^) > —0.3 GeV^) one might argue that the 
chosen form for the propagator parametrizes the momentum regime needed well 
enough for practical calculations of physical processes and only the continuation 
to large timelike momenta needs to be adapted. 

On the other hand, there are processes which permit to test the quark propaga- 
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Figure 5.3: The dominating processes in the diquark-quark model for kaon photopro- 
duction, assuming small axialvector correlations. 

tor properties for larger timelike momenta. Consider the photoproduction of kaons 
off the proton (7p K+A), measured by the SAPHIR collaboration [107]. The 
dominant contributions to this process in the diquark-quark model are depicted in 
figure 5.3. The threshold photon energy for this process to occur (with the proton 
being at rest) is = 0.913 GeV and a broad, dipped maximum for the total cross 
section is found for = 1.1 . . . 1.5 GeV. Now consider the left diagram in fig- 
ure 5.3. For the momentum p of the quark propagating between the points of the 
photon absorption and kaon emission we find the bound Re(p^) > —{7]Mn + E-y^ 
[108], with Mn being the nucleon mass and r] the momentum partitioning pa- 
rameter between quark and diquark. For photon energies larger than 1 GeV a 
propagator dressed as given by eq. (5.1) will show considerable enhancement and 
this has indeed been found [109, 108]. The calculation of the total cross section 
using the parameter set with pointlike diquarks from table 5.1 reveals an upshot 
of the cross section beyond E^. = 1.2 GeV, clearly in contrast to the data. Only 
with a (non-analytical) modification function, 

(1 |r)^ -I- 772^ I \ 

which damps the propagator for timelike momenta, the experimental cross sec- 
tion could be reproduced. Unfortunately solutions of the Bethe-Salpeter equation 
and results for observables employing such a propagator do not exhibit invari- 
ance under shifts of the momentum partitioning parameter rj anymore due to the 
non-analyticity of C. Thereby relativistic translation invariance seems to be lost. 

Another phenomenological argument against the pole-free form for the quark 
propagator comes from the analysis of Deep Inelastic Scattering in the diquark- 
quark picture. In the Bjorken limit the leading contribution to the hadronic tensor 
is obtained by calculating the imaginary part of the "handbag diagram" shown in 
figure 5.4. This diagram constitutes the leading part of the Compton tensor in 
forward scattering. By Cutkosky's rule the imaginary part is obtained by cutting 
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Cutkosky's rule 



Figure 5.4: Handbag diagram for Virtual Compton Scattering in the diquark-quark 
model. The imaginary part of the forward amplitude is obtained by Cutkosky's rule, 
i.e. by setting the intermediate quark and diquark propagators on-shell, indicated by the 
dashed Une. 

the diagram as indicated in the figure and setting the propagators on-shell at the 
cuts. This procedure is of course not possible when employing pole-free prop- 
agators and in this case the imaginary part of the Compton scattering amplitude 
could only be obtained by calculating more complicated diagrams with intermedi- 
ate mesons and baryons. In the limit of the squared photon momentum going 
to infinity, each of these diagrams will be suppressed and the hadronic tensor van- 
ishes, in obvious contradiction to the measured finite structure functions. 

Therefore we are in the peculiar situation that results obtained by using free 
quark and diquark propagators compares favorably with experiments when cal- 
culating quark distributions. For the quark-diquark model with pointlike scalar 
diquarks this has been done in ref. [57]. A simple quark-diquark spectator model 
which parametrizes the nucleon vertex function in terms of the "s wave" covari- 
ants Si and A3, cf. table 2.1, has been used in ref. [110] to calculate various 
distribution and fragmentation functions. These calculations indicate that the ap- 
proximation of the plethora of intermediate (colorless) states possible in Compton 
scattering by a simple (colored) quark and a (colored) diquark works amazingly 
well. However, in the beginning of the chapter we have emphasized that exactly 
these thresholds are sought to avoid in other calculations such as those related 
to kaon photoproduction. It therefore remains an interesting task to resolve this 
puzzle in a consistent manner. 
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Chapter 6 

The Salpeter Approximation 



In this chapter we will critically compare results for vertex functions and nucleon 
observables obtained (i) by using the full solution as described in chapters 2-4 and 
(ii) in the Salpeter approximation [66] to the (ladder) Bethe-Salpeter equation. 
The Salpeter approximation, to be explained below, is a popular 3-dimensional 
reduction of the original 4-dimensional problem. We mention that the here ob- 
tained results are summarized in ref. [25]. 

We can identify two motivations to study semi-relativistic reductions to the 
Bethe-Salpeter equation. First, the gain in technical simplicity is substantial when 
the equation is transformed into a 3-dimensional problem which can be solved 
with the tried and tested methods from quantum mechanics. Secondly, excited 
states which have negative norm are absent in the spectrum of a Salpeter ap- 
proximated equation. We will shortly dwell upon this point. Let us consider 
the spectrum of the Bethe-Salpeter equation in ladder approximation within two 
well-studied models, the Wick-Cutkosky model [111] and the positronium sys- 
tem [112, 113]. They have a class of "normal" solutions, i.e. solutions with a 
non-relativistic quantum mechanical analogue, which exist already for infinitesi- 
mally small coupling constants. On the other side, there are solutions which begin 
to exist only for a finite coupling constant, and these are called "abnormal" states. 
In the Wick-Cutkosky model, they are excitations in the relative time between 
the two constituents and clearly in a non-relativistic description with no separa- 
tion in the time coordinate between the two constituents corresponding states do 
not exist. Furthermore, these abnormal states have partly negative norm [35]. 
Therefore, the interpretation of excited states within a given model employing the 
Bethe-Salpeter equation for modelling bound states is seriously hampered. Here, 
the Salpeter approximation allows a transformation of the problem into a Hamilto- 
nian form with the Hamiltonian being hermitian and thus no negative norm states 
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can arise . This transformation is possible as the interaction kernel in the Bethe- 
Salpeter equation is assumed to be independent of the relative time coordinate or, 
equivalently, independent of the fourth component of the relative momenta {equal 
time approximation). Studies employing this technique in solving for ground state 
and excited meson masses can be found in refs. [1 15, 1 16]. 

Appreciating the advantages of the Salpeter approximation in the computation 
of excited states one could proceed and calculate dynamical observables for the 
bound states. However, we will show in the following that nucleon form factors, 
the pion-nucleon and weak coupling constants calculated in the diquark-quark 
model differ substantially for the full solution and the Salpeter approximation. 
Therefore, results from the latter might lead to a misinterpretation of the model's 
parameters. As there exist studies for the diquark-quark model using the Salpeter 
approximation [62, 117], we will employ model parameters as chosen there to 
furnish the comparison. 



6.1 Solutions for vertex functions 

We will use the diquark-quark model as defined in chapter 2. Following refs. [62, 
1 17], we employ free quark and diquark propagators, with the gauge parameter ^ 
for the axialvector diquark chosen to be 1, cf. eq. (2.46), 

Sip) - (6.1) 
-|- rriq 

D{p) = D^^{p) = - . (6.2) 

The Dirac part of the diquark-quark vertices is given as in eqs. (2.51,2.52) by 

= 9s il'C) V{q^) , x' = 9a irC) V{q') , (6.3) 

with the scalar function V chosen to be 

y(g^) = exp(-4A2g2) . (6.4) 

The relative momentum between the quarks is denoted by q. We recapitulate the 
Bethe-Salpeter equation in short notation. 



*(p,P) = S{p,)D{pa)^p,P) (6.5) 
Hp,P) = J-0^,K^^(p,k,P)^{k,P). (6.6) 



'Abnormal states may still exist, cf. ref. [114] for an investigation on the persistence of abnor- 
mal states in 3-dimensional reductions of a two-fermion Bethe-Salpeter equation. 
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All terms appearing here are defined in section 2.3.1. Let qt — q — P{q-P) denote 
the transversal part of any momentum q with respect to the bound state momentum 
P. The Salpeter approximation consists in reducing the momentum dependence 
of the quark exchange kernel in the following manner, 

K^^ip,k,P)^K^^ipT,kT) . (6.7) 

In the rest frame of the bound state, K^^(pt, kx) = K^^{p, k), and hence the 
(other) names equal time or instantaneous approximation. It follows immediately 
from eq. (6.6) that P) = $(p). Although covariance is lost by assuming 
instantaneous interactions, the solutions for $ are invariant against shifts in the 
parameter r] which defines the distribution of momentum between quark and di- 
quark. 

Again following ref. [62], the reduced Bethe-Salpeter equation can be pro- 
cessed by introducing 

np)-i' j ^^(p^p)^ (6-8) 

integrating eq. (6.5) over likewise, and thus arriving at an equation with the 
formal structure, 

= . (6.9) 

The "Hamiltonian" H is an integral operator which contains "kinetic" parts stem- 
ming from the propagators and a "potential", being the instantaneous quark ex- 
change. 

Instead of solving this Schrodinger-type equation, we will employ the numeri- 
cal method developed in section 2.5.1. Having expanded the wave function \1/ and 
the vertex function $ according to eqs. (2.91,2.129), we solve the Bethe-Salpeter 
equation for the sets of scalar functions {Si,Ai} and {Si,Ai}. As described in 
section 2.5.1, each of the scalar functions depends on and z — p-P. The depen- 
dence on z will be absorbed by the Chebyshev expansion of eqs. (2.135,2.136). 
Note that in the Salpeter approximation the scalar functions Si^ Ai pertaining to 
the vertex function depend only on the single variable = p^{l — z^), thus their 
Chebyshev expansions contains just even moments. 

As stated above we adopt the parameters of refs. [62, 117]. The first study in- 
vestigates one parameter set in the scalar diquark sector only and the second one 
includes the axialvector diquark channel using another parameter set. In these 
studies, the author did not fix the diquark-quark coupling constants by a normal- 
ization condition like in eqs. (2.59,2.60), but rather adjusted Qs so as to obtain 
the physical nucleon mass for given values of quark and diquark mass, and of 
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Set 




rrig 
[GeV] 


rrisc 
[GeV] 


[GeV] 


A 
[fm] 


9a/ Qs 


9s 




la 


Salpeter 
full 


0.35 


0.65 




0.18 




20.0 
16.5 


(20.0) 


Ha 


Salpeter 
full 


0.35 


0.65 


0.65 


0.24 


1 


12.3 
9.6 


(11.5) 



Table 6.1: The two parameter sets of the model as taken from refs. [62, 1 17]. 



9a/ 9s when including the axialvector diquarks. This procedure is convenient as 
the Bethe-Salpeter equation is an eigenvalue problem for the eigenvalue gg. In ta- 
ble 6.1 we have listed the two sets of parameters with their corresponding eigen- 
values gs obtained by us in the full calculation and the Salpeter approximation. 
The values in parentheses are the ones from refs. [62, 117]. Please note that due 
to a different flavor normalization these values had to be multiplied by to be 
directly comparable to ours. 

Although we could reproduce the eigenvalue for the model case with scalar 
diquarks only, this is not the case for Set 11a. We observe that the calculations 
of [117] involved only 4 instead of 6 axialvector components of namely the 
projected ones onto zero orbital angular momentum and the corresponding lower 
components. Still one would expect a higher eigenvalue in the reduced system. 
The more striking observation is the amplification of the eigenvalue by about 20 
... 25 % in the Salpeter approximation although the binding energy is small, be- 
ing only 6% of the sum of the constituent masses. This is in contrast to results 
obtained in the massive Wick-Cutkosky model [118] where the Salpeter approx- 
imation leads to a reduction of the eigenvalue. This may be attributed to the 
exchange of a boson instead of a fermion as here in the diquark-quark model. 

The substantial difference between the two approaches is also reflected in the 
vertex function solutions themselves. Figure 6.1 shows the Chebyshev moments 
of the dominant scalar function 5*1 for both methods using the parameters of Set 
la. Only the even momenta are given, since the odd ones are zero in the Salpeter 
approximation (but are present, of course, in the full calculation). Two things are 
manifest: the Salpeter amplitudes have a much broader spatial extent than the full 
amplitudes. Secondly, the expansion in Chebyshev polynomials that relies on an 
approximate 0(4) symmetry converges much more rapidly for the full solution 
but is hardly convincing in the Salpeter approximation. Again, since the scalar 
functions depend in the Salpeter approximation on just one variable, p^(l — z'^). 
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Chebyshev Moments of S, Chebyshev Moments of 



Salpeter approximation Full calculation 




1 2 3 4 5 ■ 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



iPl [GeV] Ipl [GeV] 

Figure 6.1: The Chebyshev expansion of the dominant scalar function <S'i(p, z) in the 
Salpeter approximation and the full calculation. 

our expansion is a cumbersome way of visualizing the solution but makes clear 
that the Salpeter approximation does not exhibit an approximate 0(4) -symmetry. 

We remind ourselves that Salpeter in his original work [66] devised the instan- 
taneous approximation of the Bethe-Salpeter kernel as the first term in a pertur- 
bation series for the full solution. We clearly see that one should not stop after 
the first term as the character of the exact and the approximated solution differ so 
profoundly. 

6.2 Results for observables 

Let us turn first to the calculation of the electromagnetic form factors Ge and 
Gm- They can be extracted from the current matrix element as prescribed by 
eqs. (3.8,3-9). In impulse approximation the current operator for proton and neu- 
tron reads 

,{Pf\Jnp^)^ = l{JT"'' + l(Jscr"'+{j^.r-^\ (6.10) 

(6.11) 

Here we adopted the notation from section 3.1. Diagrammatically the single ma- 
trix elements correspond to the left two graphs of figure 3.2. The quark and di- 
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0.2 0.4 0.6 0.8 1 1.2 0.2 0.4 0.6 0.8 1 1.2 



Q'[GeV'] ' Q'[GeV'] 



Figure 6.2: Nucleon electric form factors in the Salpeter approximation and in the full 
calculation. The curves for the Salpeter approximation have been obtained by including 
the Chebyshev moments of the vertex function up to order mmax- 

quark vertices which are needed in these diagrams are given by eqs. (3.28-3.31). 
These diagrams are the only ones which have been calculated in refs. [62, 117], 
with the following justification given: as already mentioned in the summary of 
section 3.1, the diagrams of the impulse approximation separately conserve the 
current. Furthermore, peculiar identities at zero momentum transfer hold in the 
Salpeter approximation, 

{J^{Q' = 0)y'-'' = {J^^{Q^ = 0))" (6.12) 
{J^{Q' = 0)r-^ = {JtiQ' = 0)r-^ , (6.13) 

which are proved by straightforward calculation. These identities, in conjunction 
with eqs. (6.10,6.1 1) and the normalization condition (3.56) for the nucleon wave 
function, guarantee that proton and neutron have their correct charge, i.e. the 
electric form factors satisfy G'g"°*°''(Q2 = 0) = 1 and "*"'°"^(g2 = 0) = 0. 

We have seen that in the full calculation we need the current contributions from 
the quark exchange kernel to establish the correct nucleon charges. This followed 
as a consequence of the general Ward-Takahashi identity, eq. (3.13), for the quark- 
diquark propagator. We will therefore compare the impulse approximation for the 
Salpeter calculation to the calculations with full vertex functions with and without 
the exchange diagrams of figure 3.4. 

In figure 6.2 the electric form factors of proton and neutron are displayed, 
using the parameters of Set Ila. The first observation is that in the Salpeter ap- 
proximation we could not obtain convergence with the expansion in Chebyshev 
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— l^p — 



Set 






/ jfi \ax— ax 
\'^ax/ 




SUM 


la 


Salpeter 


1.57 






1.58 (1.58) 




full 


2.04 




0.31 


2.38 


na 


Salpeter 


0.92(1.08) 


1.37 (1.7) 




2.29 (2.78) 




full 


1.09 


0.98 


0.37 


2.45 



Table 6.2: The most important contributions to the proton magnetic moment from the 
single diagrams. The notation is as in table 3.2. In parentheses the values of refs. [62, 117] 
are given. 

polynomials beyond Q"^ ^ 0.4 GeV^. As described in section 3.2.1, we com- 
puted the form factors in the Breit frame where Q is real but Zi = p ■ Pi and 
Zf — k ■ Pf, the angular variables which enter as their arguments the Chebyshev 
polynomials for the boosted vertex functions, have imaginary parts. Their abso- 
lute values may exceed one, except for the case of zero momentum transfer, cf. 
eqs. (3.66,3.67). So this expansion that works in the rest frame, and it does barely 
so for the Salpeter approximation, will not generally work in a moving frame. On 
the other hand, the decrease of the higher Chebyshev moments is good enough for 
the full four-dimensional solution that the numerical values for the form factors 
converge up to several GeV^. However, as can be seen from fig. 6.2 the Salpeter 
approximation badly fails above 0.5 GeV^ thereby revealing its semi-relativistic 
nature. 

The second finding concerns the electromagnetic radii. The Salpeter approx- 
imation tends to underestimate the proton charge radius and to overestimate the 
absolute value of the neutron charge radius, see the first two result lines of table 
6.3. The axialvector correlations tend to suppress the neutron electric form factor 
much more in the full calculation than in the Salpeter approximation. Remember 
that this suppression is due to the axialvector correlations which have been chosen 
to be rather strong for Set Ila, as Qs — Qa- 

Turning to the magnetic moments, the contributions of the various diagrams 
are tabulated in table 6.2 for the proton. Following ref. [117] we ascribed to the 
axialvector diquark a rather large anomalous magnetic moment of k, = 1.6 which 
was needed in that study to fit the proton magnetic moment. As already observed 
earlier, we could reproduce the magnetic moment for Set la, however, for Set Ila, 
the values differ and especially the coupling to the axialvector diquark is much 
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Set 




la 






Ha 








Salpeter 


full 




Salpeter 


full 








lA 


lA 


all 


lA 


lA 


all 




[fm] 


0.89 


1.00 


0.99 


0.88 


1.01 


1.01 


(^n)el 


[fm^] 


-0.28 


-0.21 - 


-0.23 


-0.06 


-0.04 - 


-0.04 


I^P 




1.58 


2.05 


2.58 


2.29 


2.07 


2.45 


(^p)mag 


[fm] 


1.04 


1.09 


1.06 


0.85 


0.99 


0.99 


IJ'n 




-0.77 


-1.02 - 


-1.65 


-0.98 


-1.02 - 


-1.26 


('"n)mag 


[fm] 


1.05 


1.09 


1.00 


0.94 


1.06 


1.05 


g%NN 




8.96 


11.71 


15.34 


6.03 


6.95 


9.36 


TtiNN 


[fm] 


1.04 


1.07 


1.04 


1.15 


1.21 


1.15 


9A 




0.93 


1.16 


1.46 


0.53 


0.65 


0.82 


ta 


[fm] 


0.93 


0.99 


1.02 


1.08 


1.14 


1.10 



Table 6.3: Some static nucleon observables. The acronym 'lA' indicates that the results 
in these columns have been obtained in the impulse approximation whereas 'all' means 
that also the diagrams of the quark exchange kernel have been taken into account. 

weaker in our Salpeter calculation. Rather more interesting is the comparison 
between the full calculations using either Set la or Set Ila: The axialvector di- 
quark improves the magnetic moment only marginally, due to the neglect of the 
anomalous transitions between scalar and axialvector diquarks. The Salpeter ap- 
proximation tends to overestimate the contribution of {Jax)^~^^ quite drastically. 

We have also calculated the pion-nucleon constant, the weak coupling con- 
stant and the corresponding radii for the two sets in a first approximation, taking 
into account only the diagrams where the pseudoscalar and pseudovector current 
couple to the spectator or the exchanged quark. Please refer to table 6.3 for the 
numerical results. Of course, the numbers for gA and g^NN obtained for the pa- 
rameters of Set Ila signal the missing contributions from the resolved axialvector 
diquark and from the scalar-to-axial vector transitions. But we want to emphasize 
another point: Already in impulse approximation the results vary strongly (on the 
level of 20 %) between the full and the Salpeter calculation with the Salpeter ap- 
proximation underestimating the couplings. But as we have seen for the magnetic 
moments for Set Ila, the Salpeter approximation overestimates them, so we arrive 
at the conclusion that observables vary quite strongly and unpredictably when 
employing the equal time approximation. In this case, it is hardly of use to relate 
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parameters to the physical content which the model hopefully describes. 

Let us summarize the results obtained in this chapter. The Bethe-Salpeter 
equation for the nucleon has been solved in a fully covariant way and in the in- 
stantaneous Salpeter approximation. As for the model with scalar diquarks only 
we have verified the results of ref. [62] whereas discrepancies remain if the ax- 
ialvector diquark is included. These differences are partly due to the fact that 
in ref. [117] not all (ground state) axialvector components have been taken into 
account. Additionally, we take our result as an indication that the calculations 
presented in ref. [117] might suffer from some minor error. 

However, the main purpose has been the comparison of observables calculated 
in the Salpeter approximation to the ones obtained in the fully four-dimensional 
scheme. The first very surprising observation is the overestimation of the Bethe- 
Salpeter eigenvalue in the Salpeter approximation. Phrased otherwise, for a given 
coupling constant the binding energy would be much too small in the Salpeter 
approximation. We have also demonstrated that the character of the vertex func- 
tion solutions in the full treatment and for the Salpeter approximation is grossly 
different, since the full solutions exhibit an approximate 0(4) symmetry which 
is absent in the Salpeter approximation. This has drastic consequences for the 
resulting nucleon electromagnetic form factors if the photon virtuality exceeds 
0.4 GeV^. Although the non-convergence of the form factors in the Salpeter ap- 
proximation in this momentum regime can be solely ascribed to the Chebyshev 
expansion which is inappropriate for the instantaneous approximation, we find 
nevertheless differences for static observables. On the one hand, different nucleon 
radii differ only mildly in these two approaches, with the exception of the neutron 
charge radius. On the other hand, one sees very clearly that the results (obtained in 
the Salpeter approximation) for the magnetic moments, the pion-nucleon coupling 
and the weak coupling constant deviate sizeably from the full calculation. 
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Chapter 7 

Summary and Conclusions 



In this thesis we have presented a framework for a fully covariant diquark-quark 
model of baryons and obtained solutions for bound states and certain nucleon ob- 
servables having employed simple ansdtze for the quark and diquark propagators, 
and the diquark-quark vertices. We derived the basic Bethe-Salpeter equations 
governing the dynamics of the quarks and the diquark quasi-particles by a suit- 
able reduction of the relativistic three-quark problem. This reduction consists in 
neglecting three-quark irreducible interactions and assuming separable two-quark 
correlations, the diquarks. Scalar and axialvector diquarks have been included 
into the description as the presumably most important diquark correlations. 

We solved the Bethe-Salpeter equation which sums up the continuous quark 
exchange between quark and diquark in a fully covariant manner, thereby ob- 
taining relativistic baryon wave functions which contain far more structure than 
non-relativistic wave functions. These wave functions have been subsequently 
employed to calculate electromagnetic form factors. Here, we constructed the 
current operator for the electromagnetic field using Ward-Takahashi identities and 
thereby maintaining gauge invariance. Covariance and gauge invariance are man- 
ifest in the numerical results. For the nucleon electric form factors we found 
good agreement with the experimental data. The magnetic moments fall short by 
approximately 15 . . . 30 %, depending on whether the A resonance is included 
in the description or not. A possible source of this failure is found in the free 
constituent quark and diquark propagators which have been employed in these 
calculations. The lack of confinement forced us to choose rather large quark and 
diquark masses to obtain a bound A, and the perturbative quark-photon vertex 
is probably insufficient to account for the observed magnetic moments. The re- 
cently measured ratio of electric to magnetic form factor of the proton restricts the 
strength of the axialvector correlations within the nucleon to be small, on the level 
of 20 %, if our results are to match the data. 

In an attempt to circumvent the problem posed by the inclusion of the A reso- 
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nance and the underestimated magnetic moments, we have investigated a possible 
parametrization of confinement by modifying quark and diquark propagators to 
render them pole-free in the complex momentum plane. These modifications al- 
lowed us to calculate the mass spectrum of octet and decuplet baryons in fair 
agreement with experiment and the nucleon magnetic moments received addi- 
tional contributions which even rendered them to be larger by 15 % than the ex- 
perimental values. Nevertheless, there are serious doubts about the justification 
of such a modification. The modified propagators have an essential singularity at 
timelike infinity which is not in accordance with general theorems about the be- 
havior of propagators in quantum field theory. Additionally, this singularity causes 
a spurious enhancement of the calculated cross section in kaon photoproduction, 
and would render the calculation of the leading twist contributions to Deep Inelas- 
tic Scattering impossible. Therefore, refined methods to incorporate confinement 
are needed and work on this problem is currently done. 

The isovector axial form factor for the nucleon and the pion-nucleon form 
factor have also been calculated. Chiral symmetry could be used to construct the 
pseudovector and the pseudoscalar quark current operators. As it turns out, chiral 
symmetry requires also the inclusion of vector diquarks besides the axialvector 
diquarks, since they form a chiral multiplet. The neglect of the former leads to a 
moderate violation of the Goldberger-Treiman relation, and the pion-nucleon and 
the weak coupling constant are found to be somewhat too large compared to their 
experimental values. Since the transitions between scalar and axialvector diquarks 
proved to be the main source of the calculated Goldberger-Treiman discrepancy 
and the large pion-nucleon coupling constant, we again conclude that axialvector 
correlations within the nucleon should be small. 

We also compared results for nucleon vertex function solutions and observ- 
ables between the fully covariant treatment of the diquark-quark Bethe-Salpeter 
equation and the semi-relativistic Salpeter approximation. The vertex functions 
for the latter differed strongly from the full solution, on a qualitative as well as 
on a quantitative level. Furthermore, we found rather erratic deviations between 
the results for static observables obtained with both methods and therefore con- 
clude that the semi-relativistic treatment is inadmissible within the diquark-quark 
model. 

Regarding the simplicity of the quark-diquark picture, the model reproduces 
fairly well the basic nucleon observables investigated here. It should be con- 
sidered a success that maintaining covariance is feasible in calculating the form 
factors and this is a good starting point to investigate more involved hadronic 
processes. For kaon photoproduction, 7p K^A, and associated strangeness 
production, pp — > pK^A, preliminary results are already available, cf. refs. [109, 
119]. 

To check the viability of the chosen approach within QCD, some more in- 
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formation about the quark propagator and the two-quark correlations are needed. 
Whereas there is justified hope that lattice or Dyson-Schwinger calculations will 
provide us with more refined results for the quark propagator in the near future, 
it is not clear at all whether a non-perturbative characterization of the two-quark 
correlations can be found. To proceed further in this direction, one probably has 
to resort in a next step to an approximate model for the quark-gluon dynamics. 
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A.l Conventions 

Throughout the text we work in EucUdean space with the metric g^'^ — 5^". Ac- 
cordingly, a hermitian basis of Dirac matrices {7, 7^} is used which is related to 
the matrices {7?), 7d} of the Dirac representation by 

l^^lD, 7 = -ilD ■ (A.1) 

Furthermore we define 7^ = — 7^7^7^7''. The charge conjugation matrix is de- 
fined by C = —7^7^. The free quark propagator S{k) is determined by the rela- 
tion, valid in the theory for a free, massive fermion, 

{0\Tq{x)g(,j)\0) = j J0-exp(ik-{x - t))S{k) , (A.2) 

We arrive at this equation, if we formulate the continuation between Minkowski 
and Euclidean space using the well-known rules 

-'ix'^ , Xm ^ X , (A.4) 
kf^ ^ ik , ^ ^ (A. 5) 

for a formal transcription of position variables xm and momentum variables /cm in 
Minkowski space to their Euclidean pendants x and k. Additionally we prescribe 
for the Euclidean momentum variables k —k. Furnished with an additional 
transcription rule for the metric tensor in Minkowski space. 



(A.6) 
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we find for the scalar and the Proca propagator. 



D{p) 



1 



+ mj 



p^ + 



(A.7) 



Four-dimensional momentum integrals, appearing in the Bethe-Salpeter equa- 
tion or in the expression for current matrix elements, are evaluated by using hyper- 
spherical coordinates. We parametrize the Cartesian components of an Euclidean 
vector q^^ in the following way. 



f sin sin 6 sin Q\ 
cos 4> sin 6 sin B 
cos 9 sin © 
\ cos © / 



(A.8) 



With the abbreviations z = cos and y — cos 9 the four-dimensional volume 
integral measure is 



Hyperspherical harmonics are defined by 



/22m (n + l)(n-/)!(/!)^ 



TT 



{n + i + iy. 



sin'eClt\{cose)Yimi9,(f>), 



(A. 10) 



where the polynomials C^l;(cos ©) are the Gegenbauer polynomials [52] and the 
Yim.{9^ 0) are the familiar spherical harmonics. For vanishing three-dimensional 
angular momentum, / = 0, the Gegenbauer polynomials reduce to the Chebyshev 
polynomials of the second kind, = and the corresponding hyperspherical 
harmonics are simply 



nOO 



[/„(cos©) 



(A. 11) 



Thus, the expansion of relativistic ground state wave functions (/ = 0) in terms of 
Chebyshev polynomials, as done in section 2.5 for the baryons, corresponds to an 
expansion into hyperspherical harmonics. 
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A.2 Color and flavor factors in the Bethe-Salpeter 
equations 

A.2.1 Nucleon 

For the derivation of the color factor, let us introduce the color singlet state of the 
nucleon as being proportional to a 3 -dimensional unit matrix (3 3 = 1 ® 8) 



|0)coior = yU>^'')AB := ^JlSAB, (0|0) = 1 . (A.12) 

Since both diquarks are in an antitriplet representation, the color part of the Bethe- 
Salpeter kernel reads (see figure A.l) 

-^TOlor = 1 ^LKA £JNA = —^{SlN^KJ — ^LJ^Kn) ■ (A. 13) 

Sandwiching between the nucleon color singlet state yields 

(0|i^e''oL|0) = -ISjk{SlnSkj - SljSkn)Sln = -1 . (A.14) 

The nucleon flavor state distinguishes of course between scalar and axialvec- 
tor diquarks. We rewrite the diquarks as eigenstates of their total isospin and the 
projection onto the respective z axis. They are easily expressed by the combina- 
tions of u and d quark, (ud) = —^J\{ud — du) has isospin zero (scalar diquark). 

[uu], [ud] = ^J^{ud + du) and [dd] have isospin one with the respective third 
{z) component being +1, 0, —1. Expressed with spherical isospin matrices Tm, 
T±i = =t ^"T"^) aiid To = r^, they read 

{ud) = — ^ir^ (A. 15) 
v2 

[uu\\ / ITiT^ \ 

[ud\ \ = —\ iror^ (A. 16) 
[dd]J V2 \,r.,Ty 

Remember the expression for the nucleon Faddeev amplitude, 

r„^, = x'p, D + x^, D^'^ . (A.17) 

The two summands on the right hand side represent scalar and axialvector correla- 
tions, respectively, and the Greek indices contain also the flavor indices. We will 
write the nucleon flavor state as a 2-vector, with scalar and axialvector correla- 
tions as components, and keep in mind that they add up to the Faddeev amplitude 
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\p)tV \^\u[ud]-^ld[uu] 



as indicated above. Then the proton flavor wave function can be constructed with 
the appropriate Clebsch-Gordan coefficients as 

1/^/ flavor — \ \^\abc 



(A.18) 



The scalar and axialvector flavor components are separately normalized to yield 
unity, 

:c[L](p|p):c[ax] = l, (A. 19) 

where the conjugate flavor states are obtained by Hermitian conjugation. The 
absolute strength of the scalar and axialvector correlations is determined by the 
solution to the Bethe-Salpeter equation. A redefinition of the diquark flavor ma- 
trices appearing in eq. (A.18) also affects the diquark normalization conditions, 
eqs. (2.59,2.60), and the Bethe-Salpeter kernel such that the contributions of the 
scalar and axialvector correlations to observables remain unchanged. The neutron 
flavor state accordingly writes 

(d{ud) \ 
-^^d[ud]^^^u[dd]) ^^-^^^ 

The matrix representation of the Bethe-Salpeter kernel with respect to the diquark 
isospin eigenstates defined above is given by (see figure A.l) 

/=0 1=1 
-^flavor = o I / o /. o^ . o I ■ (A.21) 



1=1 



kn/ 



Sandwiching between either the proton or the neutron flavor state yields the iden- 
tical flavor matrix 



7=0 7=1 

miLr\p) = MLM) = I ( "_f) . (A.22) 

Eqs. (A.14,A.22) provide for the explicit factors in the Bethe-Salpeter kernel of 
eq. (2.70). 
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Figure A.l: Color and flavor indices appearing in the Bethe-Salpeter kernel. Capital 
letters run from 1 ... 3 and label the color state of quarks (3) and of diquarks (3). Small 
letters n,k = 1,2 label quark isospin and m, m' = ±1, label the third component of the 
axialvector diquark isospin. 



A.2.2 A 



Taking over the notations from the previous subsection, we write the flavor states 
of the four nearly mass-degenerate A resonances as 



u\uu\ 



\ 



^ulud] + -i/ ^dluu] 
d[dd] ) 

The flavor matrix of the Bethe-Salpeter kernel involves only the axialvector di 
quark (7 = 1) channel. 



) flavor 



{irm'T')ir'rl)),^ ■ (A.24) 
This yields for all A states i G {++,+,0,—} identical flavor factors, 

(A«|(i^f)flavor|A«) = l. (A.25) 



A.3 Wave functions 

On the next pages the scalar functions describing the partial waves of nucleon 
and A are displayed. In figure A.2 we show the nucleon functions Si, Ai for the 
parameter Set I of table 2.8. The momentum partitioning parameter 77 = 0.36 
has been chosen. This value is roughly half way between the bounds set by the 
constituent poles, cf. eq. (2.154), 



T] e [0.334, 0.383] 



(A.26) 
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The even Chebyshev momenta dominate the expansion of the scalar functions and 
therefore we show the three leading even momenta for each function. The two 
dominating s waves §1 and A3 converge quickly in their expansion. The p waves 
^2, A2 and which follow next with respect to the overall magnitude display 
a more moderate convergence. One can estimate that these p waves will have 
an influence on observables on the level of below 10 %. All other functions are 
subdominant, especially the d wave A5 turns out to be extremely small. 

The A partial waves depicted in figure A. 3 have been calculated using the 
parameter Set II which fits the experimental A mass. The binding energy is 
only 24 MeV and consequently the range of allowed rj values is restricted by 
eqs. (2.154,2.155) to 

77 e [0.328, 0.345] . (A.27) 

We have chosen r] = 1/3. The small binding energy is reflected in the scalar 
functions as follows. The only s wave Di clearly dominates and the p waves D2, 
D4 and De represent corrections of the order of one per cent. All other partial 
waves are suppressed. Although the convergence of the Chebyshev expansion is 
still moderate for Di, the remaining functions possess a number of higher mo- 
ments which are larger than one tenth of the zeroth moment. In figure A.3 we 
have picked the moments 0, 4 and 8 to illustrate this behaviour. 

A.4 Bethe-Salpeter equations for octet and decuplet 

Here we list the Bethe-Salpeter equations for octet and decuplet baryons needed 
in section 5.2. First we introduce the notation: 



Spin- 1/2 wave function for flavor components with spectator 
quark c and scalar (ab) or axialvector diquark [ab], decomposed 
asineq. (2.129) 

Spin-3/2 wave function for flavor components with spectator 
quark c and axialvector diquark [ab] , decomposed as in eq. (2. 1 3 1 ) 

Quark propagator as in eq. (2.64) with rUg = rua 

Diquark propagators as in eqs. (2.45,2.46) with rrisc = rn^ab) and 
max = rniab] 

Quark exchange kernel with momentum definitions as below 
eq. (2.72): = -| (a, b = 1, . . . , 5) 



c(a6) ' c[ab] 

c[ab\ 
Sa 

^{ab),Df^ab] 
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As we consider the isospin symmetric case, the flavor labels a,b,c will be chosen 
among up and strange (u, s). We recapitulate the nucleon equation (2.69), 



(A.28) 



The S hyperon is described by 



/s- 


'D7\ 

(us) 








\ 
















u[us\ 


\ 








(^fn) 


"7 







-Kf 


















u\us\ 


. {A29) 






) 









The equation for the S hyperon is obtained by interchanging it <-> s. SU{?,) 
symmetry breaking couples the eightfold way state of the A to the flavor singlet, 
s{ud) + u{ds) + d{su). We introduce the flavor states /i = [d{us) — u{ds)]/ \/2, 
f2 = s{ud) and / = {d[us] — u[ds])/\/2 and the equation for the physical A hy- 
peron reads 








o-i r)-i 








s-^ (D- )-7 



^Kf ^i^f 



The equation for the A reads, cf. eq. (2.80), 



Wry 



(A.30) 



(A.31) 



By interchanging u s one finds the equation for the Finally the E* obeys 
the equation 
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(A.32) 

In eq. (A.32) the interchange u s leads to the H* equation. 



A.4 Bethe-Salpeter equations for octet and decuplet 
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Figure A.2: Chebyshev moments of the scalar functions Si and Ai describing the nu- 
cleon wave function for parameter Set I. 
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Figure A.3: Chebyshev moments of the scalar functions Di describing the A wave 
function for parameter Set n. 
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Form Factor Calculations 



B.l 



Resolving diquarks 



B.1.1 



Electromagnetic Vertices 



Here we adopt an impulse approximation to couple the photon directly to the 
quarks inside the diquarks obtaining the scalar, axialvector and the photon-induced 
scalar-axialvector diquark transition couplings as represented by the 3 diagrams 
in figure 3.1. For on-shell diquarks these yield diquark form factors and at the 
soft point {Q = 0) the electric form factors of scalar and axialvector diquark are 
equivalent to the normalization conditions (2.59,2.60). Hereby it is assumed that 
the diquark-quark vertices could be obtained from a quark-quark scattering 
kernel which is independent of the total diquark momentum, cf. also section 3.1.3 
where we discussed the equivalence of canonical and charge normalization in the 
nucleon case. 

Due to the quark-exchange antisymmetry of the diquark amplitudes it suffices 
to calculate one diagram for each of the three contributions, i.e. , those of figure 3.1 
in which the photon couples to the "upper" quark line. The color trace yields unity 
as in the normalization integrals, eqs. (2.59) and (2.60). To perform the traces 
over the diquark flavor matrices with the charge matrix acting on the quark line 
it is advantageous to choose the diquark charge eigenstates, cf. appendix A.2.1. 
With the definition 



(B.l) 
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we then obtain the following flavor factors 

(f^c)flavor = lTrT^Qy = l, (B.2) 



( A 
' 3 » 



(B.3) 



(f^f)flavor = -iTV r2g,(r^r2) = -5^0 i ■ (B.4) 

For the flavor factors of the scalar and axialvector diquarks we obtain half their 
total charges. The factor 1/2 is due to the antisymmetry of the vertices and also 
appears in the normalization conditions (2.59,2.60). 

We are left with the quark loop that involves just a trace over the Dirac indices. 
Including the minus sign for fermion loops, we find for the Lorentz structure of 
the resolved vertices 

' = ^^^^^ ^^^^^ ^'^^^^ ' ^^-'^ 

= -Tr/(|^X^(^)^fe)r^5(pOx^(^)5^(p^ (B.7) 

The quark momenta herein are, 

+ Q , Pd + kd Q , 

Pd + kd „, 
P3 = — ^ q ■ (B.8) 

Even though current conservation can be maintained with these vertices on-shell, 
off-shell rf4 ^iid r^i"'^ do not satisfy the Ward-Takahashi identities for the free 
propagators in eqs. (2.45,2.46). Thus they cannot be directly employed to cou- 
ple the photon to the diquarks inside the nucleon without violating gauge invari- 
ance. For Q — 0, however, they can be used to estimate the anomalous magnetic 
moment k of the axialvector diquark and the strength of the scalar-axialvector 
transition, denoted by in (3.32), as follows. 

First we calculate the contributions of the scalar and axialvector diquark to the 
proton charge, i.e. the contribution of the first diagram in figure 3.2 to Ge{0), 
upon replacing the vertices F^^ and F^^"^ given in eqs. (3.30,3.31) by the resolved 
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QT^ga 


K 




Set I 


0.943 


1.421 


1.01 


2.09 


Setn 


0.907 


3.342 


1.04 


2.14 



Table B.l: Rescaled diquark normalizations and constants of photon-diquark couplings. 



ones, and F^^"^ in eqs. (B.5) and (B.6). Since the bare vertices satisfy the 
Ward-Takahashi identities, and since current conservation is maintained in the 
calculation of the electromagnetic form factors, the correct charges of both nucle- 
ons are guaranteed to result from the contributions to G£;(0) obtained with these 
bare vertices, T^^ and T^^f^ of eqs. (3.30,3.31). In order to reproduce these correct 
contributions, we then adjust the values for the diquark couplings, Qs and Qa, to be 
used in connection with the resolved vertices of eqs. (B.5) and (B.6). This yields 
couplings gl^^'^ and gl^^'^, slightly rescaled (by a factor of the order of one, cf. ta- 
ble B.l). We see that gauge invariance in form of the differential Ward identity 
provides an off- shell constraint on the diquark normalization where we have aver- 
aged over diquark off-shell states by performing the four-dimensional integration 
for the current matrix element. Once the "off-shell couplings" gl^^^ and g^^^ are 
fixed we can continue and calculate the contributions to the magnetic moment of 
the proton that arise from the resolved axialvector and transition couplings, F^j,"^ 
and Tf^f, respectively. These contributions determine the values of the constants 
K and Ksa for the couplings in eqs. (3.31) and (3.32,3.33). The results are given in 
table B.l. As can be seen, the values obtained for k, and by this procedure are 
insensitive to the parameter sets for the nucleon amplitudes. In the calculations of 
observables we use k,— 1.0 and Kga — 2.1. 

B.1.2 Pseudoscalar and pseudovector vertices 

The pion and the pseudovector current do not couple to the scalar diquark. There- 
fore, in both cases only those two contributions have to be computed which are ob- 
tained from the middle and lower diagrams in figure 3.1 with replacing the photon- 
quark vertex by the pion-quark vertex of eq. (4.14), and by the pseudovector-quark 
vertex of eq. (4.16), respectively. 

First we evaluate the flavor trace implicitly given in the diagrams of figure 3.1. 
As we employ only the third component (in isospace) of the pseudoscalar and the 
pseudovector currents to obtain gT^-NN and gA, cf. eqs. (4.2,4.6), we just need the 
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following flavor factors, 

(f5,ljflavor = ^Tr (r^t^jr^ (t^t^) = oJ , (B.9) 

(f^Jflavor = -iTrTV(T„T2) = -5^o. (B.IO) 

The flavor factors for the pseudovector vertices are half of the above factors. 

In the quark loop calculations we use the rescaled couplings from appendix 
B.1.1 which proved to yield consistent results with respect to the Ward identity 
for the nucleon. For the Dirac part of the vertex describing the pion coupling to 
the axialvector diquark we obtain, 

2 

~ in 

A {gryV{q-Q/A)V{q + Q/A) 



I 



(27r)4 (p2 + ml) ipl + ml) {pi + ml 



and fixes its strength (at = 0) to ^ 4.5, see table B.2. 

For the effective pseudovector-axialvector diquark vertex in eq. (4.18) it is 
sufficient to consider the regular part, since its pion pole contribution is fully de- 
termined by eq. (B.l 1) already. The regular part reads, 

J (27r)4 {pl + ml){pl + ml){pl + ml) ^ ' ^ 

[-Amyf" (pi +P2+ p^)" - Tr757"/^27^^i7''i?^3] • 



Although after the g-integration the terms in brackets yield the four independent 
Lorentz structures discussed in the paragraph above eq. (4.18), only the first term 
contributes to ^(^(O) (with ^1+^2+^3 = (3/4) (p^ + kd) + q). 

The scalar-axialvector transition induced by the pion is described by the vertex 

= j^jr9rV{q-Q/m{q + Q/^) X (B.13) 

(P2-P3)Pl - {P3-Pl)P2 + {PrP2)p§ 

(pI + ml) {pI + m^) (pI + m2) 

and the reverse (axialvector-scalar) transition is obtained by substituting Q — > —Q 
(or pi <-> P2) in (B.14). The corresponding vertex for the pseudovector current 
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Set I 4.53 4.41 3.97 1.97 
Set II 4.55 4.47 3.84 2.13 



Table B.2: Strengths for pion- and pseudovector-diquark couplings. 



reads 



[S'''^{ml ~ pi-p2 - P2-P3 - PrPi) + (B.14) 

{pm}f + {pm}f - {P2P3Y-^' . 

The short-hand notation for a(n) (anti) symmetric product used herein is defined 
as {piP2\± — P1P2 P1P2 • The reverse transition is obtained from Q —Q 
together with an overall sign change in (B.14). As already mentioned in the main 
text, the term proportional to provides 99 % of the value for gA as obtained 
with the full vertex. It therefore clearly represents the dominant tensor structure. 

As explained for the electromagnetic couplings of diquarks, we use these re- 
solved vertices in connection with the rescaled couplings gf^^ and g^^^^ to compute 
g-KNN and gA in the limit Q ^ 0. In this way the otherwise unknown constants 
that occur in the (pointlike) vertices of eqs. (4.17-4.20) are determined. 

As seen from the results in table B.2, the values obtained for these effective 
coupling constants are only slightly dependent on the parameter set (the only ex- 
ception being k^^^ where the two values differ by 8 %). For the numerical calcu- 
lations presented in section 4.3 we employ — 4.5, = 4.4, k,^^ — 3.9 and 



B.2 Calculation of the impulse approximation dia- 
grams 

In this appendix we discuss the difficulties in the formal transition from the Min- 
kowski to the Euclidean metric. These are encountered in the connection between 
Bethe-Salpeter wave function ^ and vertex functions $ in a general (boosted) 
frame of reference of the nucleon bound state. 
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k^=(l-^}Pf-k =(l-T])P,-p 



Figure B.l: The quark coupling diagram of the impulse approximation. 

As the generic example for this discussion, we have chosen the second dia- 
gram in figure 3.2 which describes the impulse-approximated contribution arising 
from the coupling of the photon to the quark within the nucleon. Furthermore 
it is sufficient to restrict ourselves to scalar diquark correlations, retaining only 

and The principle of the calculation will nevertheless become clear. We 
have singled out this diagram and, having added the appropriate momentum def- 
initions, it is shown in figure B.l. We evaluate the diagram in the Breit frame, cf. 
the momentum definitions in eqs. (3.65). 

Let us consider Mandelstam's formalism, eq. (3.26), in Minkowski space where 
it has been derived originally [63]. Here, the matrix elements between bound 
states are related to the corresponding Bethe-Salpeter wave functions in Minkowski 
space which for the moment are our nucleon Bethe-Salpeter wave functions ^m-^ 
Upon the transition to the Euclidean metric, the corresponding contribution to the 
observable, here to the nucleon form factors, is determined by the "Euclidean" 
Bethe-Salpeter wave function ^ e- In the rest frame of the nucleon bound state 
this transfer from M — > £■ of the Bethe-Salpeter wave functions commutes with 
the replacement of the wave by the vertex functions; that is, unique results are ob- 
tained from the Euclidean contributions based on either employing the Minkowski 
space wave functions or the vertex functions which are related by the truncation 

of the propagators of the constituent legs, here $m = (Gq"*^*^)"^ 



versa, = Gq 



^M> or, vice 

M 



M 

2 



At finite momentum transfer one needs to employ Bethe-Salpeter wave 

functions in a more general frame of reference, here we use the Breit frame in 
which neither the incoming nor the outgoing nucleon are at rest. As described in 
section 3.2.1, the "Euclidean wave function" ^'^ in this frame is obtained from 
the solution to the Bethe-Salpeter equation in the rest frame by analytic continua- 



'In the following the subscript m stands for definitions in Minkowski space and e for the 
corresponding ones in Euclidean space. 
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Wick rotation 



M 



M 



E 



+ residue terms 



Figure B.2: Interrelation of matrix elements in Minkowski and Euclidean space. The 
integral sign is shorthand for the four-dimensional integration over the relative momentum 
p, see equation (B.15). 



tion, in particular, by inserting complex values for the argument of the Chebyshev 
polynomials, see eqs. (3.66,3.67). This corresponds to the transition from left to 
right indicated by the arrow of the upper line in figure B.2. 

In the analogous transition on the other hand, when the truncated Bethe-Sal- 
peter amplitudes are employed, the possible presence of singularities in the legs 
has to be taken into account explicitly. In the present example, these are the single 
particle poles of the propagators of the constituent quark and diquark that might 
be encircled by the closed path in the p'^ -integration. The corresponding residues 
have to be included in the transition to the Euclidean metric in this case, which is 
indicated in the lower line of figure B.2. 

The conclusion is therefore that the naive relation between Bethe-Salpeter ver- 
tex and wave functions cannot be maintained in the Chebyshev expansion of the 
Euclidean spherical momentum coordinates when singularities are encountered in 
the truncation of the legs. Resorting to the Minkowski space definitions of vertex 
vs. wave functions, however, unique results are obtained from either employing 
the domain of holomorphy of the Bethe-Salpeter wave functions in the continua- 
tion to the Euclidean metric (with complex momenta) or, alternatively and techni- 
cally more involved, from keeping track of the singularities that can occur in the 
Wick rotation when the truncated amplitudes and explicit constituent propagators 
are employed. 

The rest of this section is concerned with the description of how to account for 
these singularities which, for our present calculations are affected by constituent 
poles for quark and diquark, give rise to residue terms as indicated in the lower 
right comer of figure B.2. 

To this end consider the quark contribution to the matrix elements of the elec- 
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^ ^ ^ Contour 

Imp" 




Figure B.3: Location of the relevant singularities in the impulse-approximate quark 
contribution to the form factors. The relative momentum p is the integration variable in 
the loop diagram corresponding to equation (B.15). 

tromagnetic current which is given by (cf. also eq. (3.39)) 

Kr-^^ = J^-0^,^Hp+{l-v)Q:Pf)D{k,)S{k,)T^^S{p,)^'{p,R^. 

(B.15) 

We are interested in the location of the propagator poles herein (in Minkowski 
space). For these poles, solving the corresponding quadratic equation for the ze- 
roth component of the relative momentum p° yields 

<^ (B.16) 

-TjUQ ± W{m„{p-rj/2Qf) 
44> (B.17) 

-VOOQ ± W{m„{p-{v/2-l)Q)^) 
^ (B.18) 

(1 -v)ooq± W{m.,,, {p - {r]/2 - 1/2) Qf) 

with W{m, p^) = \/p'^ + m'^ — ie. We used the momentum definitions given in 
figure B.l and in eqs. (3.65). 

For Q = 0, i.e. in the rest frame of the nucleon in which the Bethe-Salpeter 
equation was solved, the naive Wick rotation is justified for 1 — ^ < r] < 



/^pole,l 



kg 



le — 



rP 

Ppo\e,2 
; 2 2 

/^polc.S 
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since there is always a finite gap between the cuts contained in the hypersur- 
face Rep° = of the Rep° - Imp" - p space. As Q increases, these cuts are 
shifted along both, the and the p°-axis, as sketched in Figure B.3. This even- 
tually amounts to the effect that one of the two cuts arising from each propagator 
crosses the Imp°-axis. As indicated in the figure, the Wick rotation p° — > ip^ is 
no longer possible for arbitrary values of p^ without encircling singularities. The 
corresponding residues thus lead to 

{JT"" - J^-0^,^'{p+{l-v)Q:Pf)D{ka)S{p,)r^^S{k,)^%p,P,) 
+ ^ J-^s(^P^'ipUe,i,P+i^-v)Q,Pf)D{k,)S{k,) X 

+ analogous terms for ^'(/Cg) and (B.19) 

upon transforming equation (B.15) to the Euclidean metric. Here, the residue 
integral is evaluated at the position of the pole in the incoming quark propagator 
S{pq) on the Euclidean p^-axis 

pUi = -^VVW+Qy4. + iW{m,, {p - rj/2 Qf) , (B.20) 
where Res(5'(A;g)) denotes the corresponding residue, and the abbreviation 

ep = e{i^ujQ- W{m,, {p - Qf)) 

was adopted to determine the integration domain for which the encircled singu- 
larities of figure B.3 contribute. 

Analogous integrals over the spatial components of the relative momentum p 
arise from the residues corresponding to the poles in the outgoing quark propaga- 
tor S{kq) and the diquark propagator D{kd) as given in eqs. (B.17,B.18). 

One verifies that these cuts (as represented by the shaded areas in Fig. B.3) 
never overlap. Pinching of the deformed contour does not occur, since there are 
no anomalous thresholds for spacelike momentum transfer in these diagrams. 

We have numerically checked the procedure described here using the parame- 
ter set from the scalar diquark sector employing diquark-quark vertices of dipole 
shape, cf. table 2.5. We have calculated the contributions to the proton electric 
form factor from the (quark and diquark) impulse approximation diagrams for 
two different values of 77, 771 and r]2. Hereby 771 = 1/3 is close to the left limit of 
the ?7 range in which the (Euclidean) Bethe-Salpeter equation can be solved, cf. 
eq. (2. 154). The second value, r]2 = 0.4, lies safely in the middle of this range. For 
each 77 value, we calculated the form factor in Euclidean space using (i) boosted 
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Proton Electric Form Factor 

myM„=0.62, nnJM„=0.7, V(p'): dipole 

2.3 I ' , ' , ' , ' , ' , ' — 

vertex function: ri=1/3 

2.1 - wave function : Ti=1/3 

vertex function: r|=0.4 
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Figure B.4: The impulse-approximated contribution, corresponding to the two left di- 
agrams of figure 3.2, to the electric form factor of the proton. We have employed the 
dipole-y parameter set from the scalar diquark sector, cf. table 2.5. Results of the Bethe- 

Salpeter wave function calculations for the t] values 1/3 and 0.4 are compared to the 
respective vertex function plus residue calculations. The results are normahzed to the 
latter that uses rj = 1/3. 

wave functions and (ii) boosted vertex functions along witii tlie residue 
prescription form eq. (B.19). The results are found in figure B.4. For rji the 
Chebyshev expansion of the Bethe-Salpeter wave function to 9 orders still turns 
out insufficient to provide for stable numerical results. This is due to being too 
close to the limit of the allowed range in rj. The considerably weaker suppression 
of higher orders in the Chebyshev expansion of the wave function as compared 
to the expansion of the vertex function enhances the residual r^-dependence of the 
observables obtained from the former expansion at a given order, in particular, 
when it has to reproduce close-by pole contributions in the constituent propaga- 
tors. The impulse approximation contributions to Ge deviate substantially from 
those employing the vertex function and residue calculations in this case. On 
the other hand, for ri2 that has been employed to give the other results of figure 
B.4, unique results are obtained from both procedures. Both the Bethe-Salpeter 
wave function and vertex function calculation are in perfect agreement for values 
of the momentum partitioning that are closer to the middle of the range allowed 
to T]. Furthermore the difference between the vertex function calculations for r]i 
and r]2 show a good agreement and only deviate for larger momentum transfer 
(Q^ > 3 GeV^) by a few per cent. Thus, the independence of the form factor on 
the momentum distribution between quark and diquark holds in actual numerical 
calculations. 
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